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ABSTRACT: An investigation of the optic spectral 
componenis beating at the photodetector to 
recover the RF signal at the output of a sub- 
octave microwave photonic link is presented, 
by a theoretical and experimental approach. It 
is demonstrated the best efficiency is achieved 
when the carrier to sideband ratio (CSR) is O 
dB at a low bias voltage condition on Mach- 
Zehnder (MZM) intensity modulator. A RF power 
link gain improvement of 9.7 dB is demonstrated 
for a same photodetector incident optic power, 
compared with a link operating at quadrature bias 
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PHOTONIC LINK 


voltage condition on the MZM. 
KEYWORDS: Fiber optic links; low bias photonic 
links, microwave photonics. 


OPTIMIZAÇÃO DA POTÊNCIA ÓPTICA 

NO FOTODETECTOR PARA AUMENTO 

DO GANHO NO ENLACE ANALÓGICO A 
FIBRA ÓPTICA DE SUB-OITAVA 


RESUMO: A análise do batimento entre 
componentes ópticas espectrais no fotodetector 
para recuperação do sinal RF na saída de um 
enlace analógico a fibra óptica de sub-oitava é 
apresentada através de abordagens teóricas 
e experimentais. Demonstra-se que a melhor 
eficiência é atingida quando a diferença entre 
amplitudes das componentes da portadora 
óptica e das bandas laterais é O dB, condição 
de baixa tensão de polarização no modulador de 
intensidade Mach-Zehnder. Uma melhora de 9,7 
dB no ganho do enlace é demonstrada neste ponto 
de operação em comparação com o modulador 
Mach-Zehnder polarizado em quadratura, com 
mesma potência óptica incidente no fotodetector. 
PALAVRAS-CHAVE: Enlace analógico a fibra 
óptica; enlace fotônico sob baixa polarização; RF 
em fotônica. 


11 INTRODUCTION 

Microwave photonics (MWP) has been 
considered a promising technology to implement 
RF functionalities at high frequency and high 
bandwidth (BW) long distance transmission lines 
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and signal processing systems, where conventional electronics circuits based on coaxial 
cables, waveguides and printed circuits boards (PCB) face limitations, due to high propagation 
losses and circuitry complexities [2]. The electronic circuit behavior is strongly dependent to 
the wavelength, and it brings severe design limitations for high bandwidth systems, where 
the start frequency is far ffom the stop frequency. In contrast, microwave photonics principle 
basically consists in modulating an optical carrier with the RF signal, to transmits and processes 
it at the optic spectrum and finally recover it at the system output by photodetection. The RF 
signal modulates an optic carrier at frequencies about 200 THz, producing a very low change 
on the wavelength among all the system bandwidth. Besides, fiber optic also offers other 
advantages like electromagnetic interference immunity (EMI), low weight and size. Those 
characteristics have been attracted attention to modern microwave systems implementation 
on civilian and military application, such as 5G, radio over fiber, electronic warfare systems, 
radar and remote antennas [7-9]. The most spread architecture of MWP systems uses Mach- 
Zehnder interferometric modulator (MZM) to intensity modulate (IM) the optic carrier and 
direct detection (DD) to recover the electric signal, commonly referred as MZM IM-DD analog 
optical link [2]. 

Just like any other RF system, the main figures of merit to evaluate an MWP system 
are the power gain (G), the noise figure (NF) and the spurious free dynamic range (SFDR). 
On MZM IM-DD systems, the gain and the NF are strongly dependent of the incident 
photodetector optic power (POP). The higher the POP, the higher the gain. But, unfortunately, 
the photodetector output noise power increases as the POP grows up. Also, in order to avoid 
harmonic distortion, the MZM is commontly biased at the quadrature point of its electrooptic 
transfer function, where the second-order harmonics are completely suppressed [2]. 

The MZM IM-DD link gain depends on the squared optic power and there is no 
theoretical superior limit to it. However, in practice, it is limited by the photodetector optic 
saturation power, limited on few milliwatis for high speed commercial devices. On quadrature 
biased MZM photonic links, the optic carrier to sideband ratio (CSR) is high [8], approximately 
20 dB for low signal condition. So, most of the power relies on the optical carrier, which brings 
the photodetector to the saturation region with no significant contribution to recover the RF 
signal, since it results from the beating between the carrier and the optic sidebands [5]. 

To improve the RF recovering efficiency at the photodetector, the CSR reduction is 
mandatory. Some approaches use nonlinear optic effects to transfer power from the carrier 
to the sideband, such as four-wave-mixing [7] and parametric amplification of the sidebands 
[9]. Other approaches move the MZM bias away from the quadrature point to get the carrier 
power down, known as low biasing technique [5][6][10]. Those approaches have shown good 
improvement on the system gain, NF and SFDR, comparing to quadrature biased MZM links. 
As consequence, the system is limited to a sub-octave bandwidth because of the increase on 
the harmonic distortion. Previous works mentioned above explore those techniques focused 
on the system performance, not considering the photodetector optic power as a reference 
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variable. This work investigates the low biasing technique to obtain gain improvement, 
considering the POP as a constant parameter, compared with the quadrature condition. 

Section 2 shows the analytical analysis of the system. Section 3 shows simulations 
using the expressions deduced in the previous section, with some approximations and 
considerations. In section 4, experimental results are presented as a proof of concept. 
Conclusions are presented in section 5. 


21 THEORETICAL FOUNDATION 


The schematic diagram of an MZM IM-DD photonic link is shown on Fig. 1. Itis 
composed basically by a laser as the optical source feeding a Mach-Zehnder intensity 
modulator, followed by a photodetector connected to the MZM by a single-mode optical 
fiber. To get a general formulation, it is considered a dual drive configuration (DD-MZM), 
where two independent RF signal are injected into the modulator together with a DC bias 
voltage to stablish the operation point at the modulator electrooptic transfer function. 

Itis assumed the lightwave electric field can be expressed by: 


e(t) = Re[E,eilvot+bolg] (1) 

where E, is the electrical field amplitude, w is the optic angular frequency, q, is the 

initial phase and € is the electric field polarization vector. For a sake of clarity, it will be 

assumed the field polarization does not change through the optic system and can be negleci 

on the formulation, which will consider only phasor representation ffom now on and will be 
referred as optic signal. 

The intensity modulation process at the DD-MZM will produce at the output an optic 
spectrum composed by an optic carrier surrounded by a series of sidebands spaced by 
multiples of the RF modulating frequency. Neglecting all the losses in the modulator and in the 
fiber optic path, it can be expressed by [3]: 


Ea: 
En(t) — Rida 


n=+00 


a da (m) cos(Ppa Rm NPma) e lwo+nwm)tAnAms] (2) 


n=-00 
where q, is the absolute phase shift accumulated through the system, P = (Vit Vico)/2 Vis 


the common phase produce by the bias voltages V.. and V.., 


applied to the MZM arms 1 and 
2, respectively, V is the half wave voltage of the MZM, J (m) represents the coefficients of the 
Bessel functions of first kind of order n, m = 7Vn/V, is the RF modulation index, V. Is the 
modulating signal amplitude, Q = Vi Vsco)/2V, Is phase difference between the MZM arms 
produced bias voltages, ? =(b,-4.,)/2and 4 =(P  +$,,)/2 are optic phases produced by 
the two RF modulating signals, where 9 | and À , are its initial phases respectively and w 


is its angular frequency. 
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The total power of the spectrum represented by (2) is the sum of the individual powers 
of each optic component [1]. The optic power at the photodetector, P., can be expressed by: 


n=+00 


Pa =PiKpo D Un(m).cos(Ppa + nbma)|” (3) 


n=-00 


where P, is the laser optic power and K is the link coefficient loss, including the 
insertion loss of the DD-MZM. 

The instantaneous envelope power of the modulated signal is proportional to the 
product of (2) by its complex conjugated and it is related to its optic intensity [3]. It can be 
understood as the average power within one period of the carrier lightwave. The RF signal at 
the photodetector output is proportional to this envelope power [2]. From (2), the photodetector 
envelope power, P,, can be expressed by: 


p=+00 


1 
Pi) = Pikpo 15 Jp“ Gm[cos(2bpa + 2pbma) + 1] 


pis uo 


k=+00 p=+00 


1 
EA E a >. Jo (M)Jp-e mM |cos(2Ppa + 2pbma 


k=1 p=-00 


Fa kPma) ei(kwmt+kbms) 


+ (—Dcos(2Ppa — 20Pma + kPma e” 'Komt+kbms)] 


k=+00 p=+00 


1 
+5 DD JoMp-emcos(kbma) [08 Omttóms) 
k 


=1 p=-00 


TE e ttenito 
(4) 


where p represents the orders of the optical spectrum components, and kthe orders 
of the RF spectrum components, in other words, the electric harmonics frequencies. 
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Fig.1 - Microwave photonic link with a DD-MZM configuration and direct detection. 
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The small signal RF power gain can be calculated taking the fundamental harmonic 
frequency for m << 1. Considering the input and output impedances matched, it can be 
expressed by [4]: 


1 TP TAViN | 
G = 7 Kpo Na? cos ( V ) ZmtL (3) 
IA TM 


where n, is the photodetector responsivity, AV. is the DD-MZM voltage bias drift far ftom 





the quadrature bias point, expressed by AVac = [(Vac — Vacz) — Vr/2], Z. is the DD-MZM input 
impedance and Z is the load output impedance. 


31 SIMULATIONS 


It was considered the DD-MZM in a push-pull operation, where V.. = -V.,and the 
phase difference between q | and P ,ism. The DD-MZM V adopted for the simulations is 
3.5 V. For simulation, the phases É and à, were considered equal to O to simplify (2) and (4). 

lt was considered V = 0.25 V to keep the modulation index on small signal condition, 


resulting in m = 0.224. 
a. Bias voltage at the quadrature point 


Based on the previous considerations, the electric field amplitude spectrum from (2) 
is expressed by: 


En(n) = - Un(m) cos ( + n5)| 


which has the shape shown in Fig. 2, per unit of electric field. 


(6) 


Using (4) to plot Fig. 3, it is possible to realize that the second harmonic does not 
appear on the expected RF amplitude spectrum. 
Appling the quadrature condition in (3) results in: 


n=+o00 


P;=P; >. Un(mcos (+ n5)] (7) 


n=—00 


which reveals that, together (6), itis possible to decrease the carrier, n=0, and increase the 
sidebands, n= +1, by growing up the cosine bias argument Do toward to m1/2. Fig. 4 shows 
that condition. Equation (2) and the red lines in Fig.4 reveal that there is a high CSR penalty 
to eliminate the second-order harmonics, 
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Amplitude Spectrum of Modulated Signal in Quadrature 
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Fig.2 - Optic amplitude spectrum of the modulated signal, at quadrature bias condition and m=0.244. 
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Fig.3 - Optic power envelope spectrum per unit of P , at quadrature bias condition and Ko unitary. 


.e, the cosine amplitude factor in (6) and (7) is the same for the carrier and sidebands 


components. For low modulation indexes, the carrier power is much bigger than the 
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sidebands and does not contribute to recover the RF output signal efficiently. This effect 
goes down with the modulation index increasing, but it keeps still high for practical condition. 
Consequently, increasing the laser power will bring the photodetector to saturation without 
produce link RF power gain efficiently. 


b. Shifting bias voltage far from the quadrature point 


By inspection of (2) and starting from the quadrature condition, the best bias voltage 
to optimize the beating between the carrier and the sidebands can be calculated for a given 
modulation index. The best beating efficiency will occur when the amplitude of the two 
components are equals. From (2) this bias condition can be determined, resulting in: 


Ouitopm = tan (e) 


For the same modulation condition adopted for the quadrature case, m = 0.224, p, 


(8) 


d.opm 
is 1.46 and the amplitude spectrum of (2) becomes: 


Es T 
En(n) = 7 Un (m) cos (1.46 — n5)| (9) 
which are shown in Fig. 5. 

The blue dotted lines in Fig. 4 also shows, together (9), that the Bessel functions 
coefficients continues with a great difference between then, but the cosine amplitude factor 
attenuates the carrier and increases the sidebands, bringing the optics components to the 
same amplitude to achieve the highest beating efficiency. 


The total power incident at the photodetector is calculated from (3), expressed by: 


n=+00 


P.=P, >. Un(m).cos (1.46 + no) (10) 


n=-00 


The photodetector output RF signal spectrum is shown in Fig. 6, where it can be 
seeing the presence of the second harmonic. 
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Fig.4 - Trigonometric circumference plotting of the cosine argument (ppd+ n m/2). 


Red lines represent the quadrature condition. Dotted blue lines represent the low bias condition. 
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Fig.5 - Optic spectrum amplitude of the modulated signal, at optimum MZM bias condition. 
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Envelope Optic Power Spectrum in Low-Bias 


D 
o 
o 
Co 


Amplitude per P. unit 


1 2 S E o o) f 5) 9 
N.vs 
m 


Fig. 6 - Optic power envelope spectrum per unit of P and Ko unitary, at optimum MZM bias condition. 


c. Link RF power gain 


Fig. 7 shows a comparative plot of the RF power gain of the link, considering the two 
conditions assumed above and equation 5. It was considered Z =Z = so), and Ka and n, 
unitary. As can be seen in Fig. 7, the gain for optimal bias condition is 9.7 dB above the gain 


in quadrature, for the same incident optic power at the photodetector. 


41 EXPERIMENT 


Fig. 8 shows the schematic diagram of the experiment setup. Similarly to simulations, 
the experiment was carried out with the MZM operating upon a push-pull configuration. The 
optic carrier is ffom a DFB laser, tuneded at the wavelengh of 1552.54 nm, with a maximal 
power of 100 mW when drived by a current of 454.72 maA. 

A laser drive controls the temperature and the current in the DFB laser. The optical 


carrier is coupled into a Mach-Zehnder intensity modulator with 10. 
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Low-Bias Point vs Quadrature Bias Point 
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Fig.7 - Comparative plot of the RF link power gain. Quadrature condition in purple line and low bias 
condition in blue line. 


GHz of bandwidth, ffom EOspace Inc, with a mesured insertion loss of5 dB, a RF V. 
of 6.54 V, a bias V of 6.3 V and input impedance of 50 . A symetric DC power supply of +/-10 
V feeding a potenciometer was used to regulate the bias voltage into the MZM, measured 
by a digital voltimetre. An RF signal generator feeds the MZM with the modulating signal. 

To measure the optical spectrum, it was inserted before the photodetector a 90/10 
fiber coupler, taking 10% of the modulated optic signal to an optical spectrum analyzer (OSA) 
and injecting 90% into a high speed photodiode, ffom New Focus Inc, with responsivity of 
0.6 A/W, saturation power of 8 mW, output impedance of 50 O and bandwidth of 25 GHz. A 
multimeter was used to measure the incident power on photodiode, using a low frequency 
transimpedance amplifier available in the photodetector module, with 1 kO of gain. The RF 
photodetector output signal is injected into an RF amplifier, ffom Miteg Inc, to feed an electric 
spectrum analyzer (ESA). The RF power gain of the amplifier is 40 dB at 7 GHz, with an input 
and output impedances of 50 O), powered by a DC power supply of 15 V. The fiber path is a 
single mode fiber and the total losses are 3.3 dB. 
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Fig.8 - Schematic diagram of the experiment setup. 


To achieve an RF signal amplitude of 0.25 V, the same value used for the 
simulations, the signal generator power was adjusted to 3.44 dBm at a frequency of 7 
GHz. The first gain measurement was taken for the low bias condition. Measuring the 
optic signal on the OSA, the MZM bias voltage was adjusted to geta CSR of O dB as 
shown in Fig. 9 (left). 

So, the laser power was set to its maximum, P = 100 mW. The transimpedance 
amplifier output marked 48 mV at this point. The fundamental RF frequency output 
signal measure was -12.98 dBm, as shown in Fig. 9 (right). Taking out the RF amplifier 
gain, the intrinsic link gain was -52.98 dB. 

Then, the MZM bias voltage was adjusted to the quadrature point, where the 
second harmonic is canceled. The laser power was decreased to obtain the same value 
at the photodetector for the low bias condition, that is, 48 mV at the transimpedance 
output. On that condition, the laser power was 2.57 mW. The optic spectrum was 
measured, and it is shown in Fig. 10 (left), with a measured CSR of 16.12 dB. 

The fundamental RF output signal is -20.61 dBm, and its spectrum is shown in 
Fig. 10 (right), with no second harmonic observed. The intrinsic link gain is -60.61 dB, 
that is, 7.63 dB below the low bias condition gain. Due to the MZM electrooptic transfer 
function instability and the absence of an automatic bias control circuit, there was an 
uncertainty on the measurements and was considered as about 2 dB, after several 


sequential observations. 
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Fig.9 — Left: optic modulated signal spectrum at low bias condition. Right: RF output signal spectrum at 
low bias condition. 
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Fig. 10 — Left: Optic modulated signal spectrum at quadrature bias condition. Right: RF signal output 
spectrum at quadrature bias condition. 


5 | CONCLUSION 


The beating of the optic spectral components at the photodetector was investigated 
to optimize the power efficiency to recover the RF signal at the sub-octave link output, 
considering the photodetector saturation power as the bottleneck of the system. It was 
demonstrated the best condition is when the CSR is O dB, and there is a RF link power gain 
improvement of 9.7 dB for the same photodetector optic power. 
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ABSTRACT: A flare system in oil platform is a 
combustion stack used to burn off excess gases 
that cannot be processed and gases that have 
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to be eliminated in emergency shutdown to avoid 
the risk of explosion. Gas flow measurement in 
flares is still considered challenging because 
the measurement has to attend many specific 
demands, quite different ffom any other flow 
measurement applications. For these types of 
measurements many sensor technologies are 
available, however they are all very expensive, 
and only a few ones can attend all demands. 
The objective of this letter is to present a flow 
rate sensor based on Fiber Bragg Grating (FBG), 
cross-correlation and heatwave travel time 
techniques. The system was developed with 
inexpensive components and is litile-intrusive, 
capable of attaining the rangeability of the flare 
demand and independent of gas composition, 
pressure and temperature. 

KEYWORDS: Cross-correlation, Fiber Bragg 
grating (FBG), Flare systems, Flow rate, Optical 
fiber sensor. 


SENSOR OPTOELETRÔNICO APLICADO 
A MEDIDAS DE VAZÃO NA INDÚSTRIA 
DE ÓLEO E GÁS 
RESUMO: Um sistema de alívio (Flare) na 
plataforma de petróleo é uma chaminé de 
combustão usada para queimar gases em 
excesso que não podem ser processados e gases 
que devem ser eliminados em desligamentos 
de emergência para evitar o risco de explosão. 
A medição do fluxo de gás em lares ainda é 
considerada um desafio porque a medição 
deve atender a muitas demandas específicas, 
totalmente diferente de qualquer outra aplicação 
de medição de fluxo. Para esses tipos de 
medição, muitas tecnologias de sensores estão 
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disponíveis, porém todas são muito caras e apenas algumas podem atender a todas as 
demandas. O objetivo deste trabalho é apresentar um sensor de vazão baseado em Redes 
de Bragg (FBG), técnicas de correlação-cruzada e tempo de viagem por ondas de calor. 
O sistema foi desenvolvido com componentes de baixo custo e pouco intrusivo, capaz de 
atingir a rangeabilidade da demanda do flare e independe da composição do gás, pressão e 
temperatura. 

PALAVRAS-CHAVE: Correlação-cruzada, Rede de fibra de Bragg (FBG), Sistemas de alívio, 
Vazão, Sensor à Fibra Óptica. 


11 INTRODUCTION 


The oil and gas industry is responsible for a significant share of the greenhouse gas 
emissions. In the exploration and production of natural gas and oil, the gas which is burned in 
torches, known as flare or relief system, is a huge source of CO, emissions. 

A typical flare system is a combustion stack used to burn off excess gases that cannot 
be processed. It is also used as a safety measure in emergency shutdown when gases in the 
process plant have to be eliminated to avoid the risk of explosion. 

In addition to contributing to global warming and climate change, the burning of natural 
gas is considered a waste of valuable, nonrenewable energy resource. 

In this context, the need to quantify the waste of gas volumes correctly and accurate 
is evident. Due to the importance of this fact, frequent publications of regulatory directives are 
issued, aiming at improving the measurements of gas flow in flare systems. 

Nowadays, the gas flow measurement in flares is performed with very expensive 
equipment and it is still considered challenging because the measurement is quite different 
from other flow measurement applications. 

The following challenges are still to be overcome when measuring flow in flares 
(Shannon, 2017): 


* | Unpredictable nature of gas flaring, since the petroleum gas is a varied mixture 
of gaseous hydrocarbons whose main component is always methane CH. 


* | NVelocity range and accuracy. During 98% of total operating time the range may 
vary from a meager 0.03 m/s up to 5 m/s, which is the normal plant operation 
(known as pilot or purge). Accuracy demands are around 10%; 


* | Gas temperature may vary from -70ºC to above 150ºC; 
* Size of flare systems. The burner line diameter may vary from 8” to 40”; 


* Being the flare an important safety equipment, flowmeters should not impinge 
any impedance in the flow so that conventional flowmeters based on turbine, 
orifice plate, Coriolis and Vortex cannot be applied on flares. 


Several technologies are available for flow velocity measurement. The most used 
technologies are: 
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* | Ultrasonic Flow Meter. It measures the difference in transit time of pulses that 
travel from a downstream transducer to the upstream transducer and vice ver- 
sa. Itis highly accurate but presents a high cost per installation, ranging from 
$50,000 - $100,000 (Sage Metering, 2020); 


*- Pitot Tube. This is the oldest way to measure flow velocity, it is based in a diffe- 
rential pressure to determine the flow velocity. It presents low flow sensitivity and 
its response depends on the gas composition; 


* | Thermal mass flow meter. Their working principle is based on the measurement 
of the heat convection from a heated surface to the flowing fluid as the sensor 
temperature is proportional to the mass flow rate. This the most commontly me- 
thod to measure the mass flow of clean and unmixed gases or gases mixtures 
when the gas composition is consistent and known. Since this technique is ba- 
sed on the thermal properties of the fluid it is dependent upon fluid composition 
(Steinberg, 2013); 


- Optical flow meters. They are divided into two classes, Laser Doppler veloci- 
meters (LDV) and Doppler-Laser-Two-Focus Velocimetry (L2F). LDV is based 
in the Doppler shift of a laser beam to measure the velocity where two beams 
of laser cross the flow being measured generating fringes. L2F method, also 
known as Laser Transit Anemometry, measures the time of flight (ToF) of parti- 
cles crossing the two laser beams (Ruck, 1994). These two optical methods are 
independent of pressure, temperature and gas composition but are based on 
the presence of sub-micrometer size particles in the gas that shall be inserted 
by the system (Parker, 2007). 

In conclusion from the challenges and technologies listed above, all flowmeters 
present good performance in some aspects but bad performance in others. Normally, the 
choice goes to the most expensive one, the ultrasonic flowmeter. 

Another important feature that is preferable in equipment applied to the oil and gas 
industry is the capability of working without embarked electricity. This is because electrically 
passive equipment does not depend on the expensive security technology for protection 
against fire and explosions. In this sense, fiber optic sensors tend to be cheaper and less 
complex to install and maintain. 

The objective of this paper is to present a flow rate sensor based on Fiber Bragg 
Grating (FBG), cross-correlation and heatwave travel time techniques. The system is 
inexpensive and minimum intrusive, capable of attaining the rangeability of the flare demand 
and independent of gas composition, pressure and temperature. 


21 CROSS-CORRELATION TECHNIQUE 


The cross-correlation (CC) is frequently applied when one needs to recognize a 
short-duration signal inserted in a longer signal. CC presents many applications, such as in 


pattern recognition, for instance in biometric recognition of retina or fingerprint, or tracking a 
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person through a camcorder through face recognition (Wu, 2010). Other application of CC 
is presented in the work of (Wenbiao, 2009), who studied cross-correlation for application in 
velocimetry using electrostatic sensors. The result of the CC algorithm is a measurement of 
similarity between two signals against the displacement of one of them along the time. CC is 
used in signal processing and is defined as the correlation of a series against another series, 
shifted by a particular number of samples. 

The CC algorithm is easily understood when we exam the algorithm applied to produce 


Suppose two time series (X,, Yo having M samples each by Equations (1) and (2). 


[xy] = [Xo X1, X2, es XM—1 | (1) 


[17x] = [0,71 Y2, Es RR (2) 


Then, the cross-correlation function, Q, (t) , Is defined as Equation (3). 


M-—1 


Pyx(T) ni >. XrYt+r (3) 


t=0 
Where + is the time delay. 


31 FIBER BRAGG GRATING 


Fiber Bragg grating is a modulation of the refractive index of the core of an optical 
fibers, forming a grating. This grating reflecis a narrow band of the light guided in the fiber, 
centred at the Bragg wavelength, according to the Bragg law (Werneck, 2017) of Equation 


(4). 
Ap = 2Nery A (4) 


Where A, is the Bragg wavelength, n. is the effective refractive index of the core and 


ff 
A is the periodicity of the grating. Essentially, any external agent capable of changing A or 
n egg Wi displace the reflected spectrum centered at À,, either a longitudinal deformation or a 
temperature variation. The sensitivity of the Bragg wavelength with temperature is given by 


Equation (5). 


Ap 
Ei (arBç +N)AT 
B (5) 
Where AA, is the Bragg wavelength shift, AT is the temperature variation, ca; is the 


silica thermal expansion coefficient aa, = 0.55 X 10% “CC and n is thermo-optic coefficient 
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for a Ge-doped silica optical fiber n = 8.6 X 10%/“C. Thus, the sensitivity of the grating to 
temperature at the wavelength range of 1550 nm is AA,/AT= 14.18 pm /“C. 


41 MATERIALS AND METHODOLOGIES 


The general idea of the proposed system it to apply a current signal on a heat source 
upstream inside the duct. The heat source is initially comprised by a small resistor, which will 
be changed in the future by a laser source injecting light into an optical fiber with its tip coated 
with an absorbing material to generate heat. In this work we tested only the resistor as a heat 
source; the methods to produce the heat with laser are under study. As a result of the heat 
pulse, a sub-degree temperature variation will flow inside the duct at the same speed as the 
gas. When this temperature pulse reaches the sensors, they will respond accordingly with a 
time delay. But due to the fact that the heat produced will be vanishing along the tube, the two 
sensors will detect different signal intensities. In order to detect the time delay, a CC algorithm 
is employed and the flow velocity is calculated by the ratio between the sensors distance and 
the time delay. 

This technique was mentioned by the first time by (Xu, 1996) that used thermistors to 
measure flow velocity. They tested at flow speeds up to 30 mm/s, clearly far below from what 
is needed for flare flow measurement. 

lt is good to notice that (Ashauer, 1993) demonstrated a ToF technique using 
thermocouples, but the measured velocities were in the range of millimeters per second, not 
obtaining a linear relationship for flow velocities above 16 mm/s. 

Additionally, (Fernandes, 2010) applied the ToF and CC techniques to measure 
newborn air breathing velocity using two lasers beams perpendicular to the respiratory flow 
and two position sensor detectors (PSDs) as sensors. They demonstrated a good linearity in 
the flow velocity measurement up to 1 m/s. 

The schematic diagram of the proposed system in this paper is shown in 
Figure 1, where dis the distance between the two temperature sensors, FBG, and FBG,. 

The CC detecis the few-milliseconds delay between the two pulses and the flow 
velocity can be calculated by the ratio between the sensors distance d and the time delay T. 

The airflow is obtained by a turbine that controls the airflow velocity in a wind tunnel 
(see Figure 1), by the motor speed and a Pitot tube is used for calibration. 

The FBG interrogator is a Micron Optics model sif55 with 5 kHz interrogation frequency 
and +2 pm resolution. 
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Figure 1 - Schematic diagram of the system inside the wind tunnel. 





The Figure 2 shows a picture of the experimental setup, the wind tunnel and the 


sensors. 
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Figure 2 - Picture of the wind tunnel and the sensors location. 


5 | RESULTS AND DISCUSSION 


The Figure 3 shows the signals collected from the FBGs for a flow rate 
of 5 m/s. The blue line represents the output from the proximal FBG, while the orange line 
represents the output from the distal FBG. The signal-to-noise ratio is low because the 
temperature pulses collected by the FBGs are too close to the accuracy of the interrogation 
equipment, +2 pm, recalling that 1ºC in an FBG represents about 14 pm displacement of its 


central wavelength. 
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Figure 3 - The output signals produced by the two FBGs at a flow rate of 5 m/s. The signal-to-noise ratio 
is too small to allow any temperature variation distinction. Blue line: FBG,; 
Orange line: FBG,. 


As expected, no direct cross-correlation is possible with this low signal-to-noise ratio, 
as the noise produced by the interrogation equipment modifies differently each signal, making 
them absolutely uncorrelated. In order to clear the signals from the noise, a low-pass, first- 
order Butterworth filter was applied using Matlab software. The command used is [b,a] = 
butter (n,W.) which returns the transfer function coefficients of an n” order lowpass digital 
Butterworth filter with normalized cutoff frequency W.. 

The command returns b and a which are row vectors of length n+1, representing the 
coefficients of the filter's transfer function. We decided by a second-order filter (n = 2) and the 
cut-off frequency to be W = 850 mHz because the heat pulse frequency was 500 mHz. Since 
n =2, vectors a and b will contain 3 numbers each, that are used on the transfer function of 
the digital filter, at Equation (6). 


ie B(z)  b(1) + b(2)z7* + b(3)z7? 
= 40)" a) ra raG (6) 


However, since filters change signal phase and since the signals shown in 
Figure 3 are different ffom each other, even using the same filter one can change differently 
the phase in each signal, disturbing the cross-correlation results. In order to circumvent this 
effect, the MatLab function filtfiltwas used. The syntax is y = filtfilt (D,a,x) that performs a zero- 
phase digital filtering by processing the input data, x, which contains 125,000 points, in both 
forward and reverse directions. 

The result of this operation is a zero-phase distortion with a transfer function equals 
to the squared magnitude of the original filter and an order that is double the order specified 
by band a. 

The Figure 4 shows the result of these two operations applied in real time to the signals 
shown in Figure 83. 
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Figure 4 - Temperature signals captured by the two FBGs as a result of the signal processing 
techniques applied to the signals shown in Figure 3. The upper line (blue) represents the temperature 
sensed by FBG1 and the bottom line (black) the temperature sensed by FBG2. Horizontal axis snows 

time (h:m:s). 


Notice that the delay between the signals is now visible and that, while the 
proximal sensor measured a pulse of about 1ºC, the distal sensor measured only 
about 0.2ºC. This temperature, according to the FBG sensitivity to temperature 
shown above, represenis a Bragg wavelength shift of only 5 pm, too close to 
+2 pm, the FBG interrogator accuracy. 

Before applying the CC procedure, itis necessary to perform another signal processing, 
because the two signals shown in Figure 3, although with a noticeable phase delay, are 
completely different ffom each other, both in shape and in amplitude. Any attempt to cross- 
correlate these signals would result in a triangular-shape output with zero delay, meaning no 
correlation at all. To circumvent this, we have to normalize both signals by the average and 
the standard deviation using the Matlab command: T = (a-mean(a))/std(a) where a stands for 
each series. 

Finally, Figure 5 shows the result of the CC algorithm performed in Matlab on the 
signals shown in Figure 4 after normalization. The maximum probability occurred at a delay 
of 135 ms. 
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Figure 5 - The cross-correlation results. Y-axis plots the probability of similarity whereas 
X-axis plots the time delay. 
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From this calculated delay, corresponding to the lag time that the heat pulse takes to 
travel 77 cm from the upstream sensor to the downstream sensor, itis possible to calculate an 
airflow velocity of approximately 5.70 m/s. 

This process was repeated twenty times each, for different velocities resulting in the 
plot shown in Figure 6. 
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Figura 6 - Plot of the measured velocity by the FBG system versus calibrated velocity measured by the 
Pitot tube. 


Notice that, while for flow velocities up to 6 m/s the standard deviations are smaller 
than 5%, for higher velocities, the error bars increase substantially reaching up to 10% for 
speeds higher than 9 m/s. The reason for this is attributed to two limitations: The first one is 
the sampling rate provided by the commercial interrogator. It produces output values at regular 
frequency, rather from being continuous, which causes phase errors in the CC calculations. 
The second reason is its +22 pm Bragg wavelength uncertainty, leading to limitations in 
detecting sub-degree-centigrade temperatures. 

High-precision ToF measurement is the main request to flow rate measurement of 
ultrasonic gas flowmeters and also for the system presented in this work. However, the 
measurement accuracy of CC depends on high-speed of an analogical-to-digital converter 
(ADC), which is not the case of commercial interrogators. 

In order to improve the accuracy and resolution of ToF measuremenis, 
(Sun, 2019) proposed the cross-correlation calculation performed after digital filtering and 
interpolation. Such technique worked well in our case but only up to 10 m/s. For higher speeds 
the conclusion is that, either the sampling rate needs to be implemented or an interrogation 
system with continuous output signal be applied. 

To circumvent these limitations, we propose the use of a bench-top interrogation 
technique such as those employed in (Dante, 2019) consisting of wavelength-division 
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multiplexing (DWDM) or Fabry-Perot edge filters, respectively, which provide a higher output 
sensitivity to FBG center wavelength displacement, together with a continuous, real-time 
output. 


61 CONCLUSIONS 


A new methodology to measure volumetric flow rate inside flare ducts based 
in heatwave travel time, FBG and cross-correlation techniques was presented. Results 
demonstrate that the system attend the main demands for a flare flowmeter showing 
independence of gas composition, good linearity and accuracy and low cost as compared 
to conventional flowmeters. On the top of that, the system can be built so as to not require 
electricity for its operation, promising to be a good candidate for an all-fiber explosion-proof 
system. 

The system presented a linear behavior for flow velocities measurement up to a velocity 
of 10 m/s. Flow velocities above that present an up to 10% precision, due to limitations of the 
commercial interrogation system used. 

It is suggested an edge filter interrogation technique to be applied in future tests to 
interrogate the FBGs in order to improve signal-noise-ratio, allowing the measurement of 
higher velocities. 
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RESUMO: Foram produzidos, caracterizados 
e mantidos em operação dois sensores óticos, 
consistindo de um sensor FBG (rede de Bragg em 
fibra Ótica) e de um sensor RFBG (rede de Bragg 
regenerada em fibra ótica). Ambos os sensores 
foram produzidos em fibra ótica monomodo 
padrão G-652, dos quais o sensor RFBG foi 
produzido em fibra hidrogenada e submetido a 
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POR 16 SEMANAS 


tratamento térmico para regeneração. Os dois 
sensores foram encapsulados com tubo capilar 
em aço inoxidável (AISI 304), calibrados em 
temperatura na faixa entre 5 “C e 60 ºC. Os 
respectivos sensores foram fixados em chapa 
metálica (parede), para medir temperaturas 
nessa faixa de operação, com o intuito de 
comparar durabilidade (envelhecimento), bem 
como precisão e exatidão em relação a um 
termopar tipo K. Ao final de 16 semanas, não 
foram observadas degradações significativas no 
sinal de cada sensor e o erro de medição foi < 
AO 

PALAVRAS-CHAVE: FBG; RFBG; Durabilidade 
de rede de Bragg; Encapsulamento. 


ANALYSIS OF AGING, PRECISION 
AND ACCURACY IN FBG AND RFBG 
OPTICAL SENSORS THAT MEASURE 

TEMPERATURES BETWEEN 5 ºC AND 60 
ºC FOR 16 WEEKS 


ABSTRACT: Two optical sensors were produced, 
characterized and maintained in operation, 
consisting of an FBG sensor (Fiber Bragg Grating) 
and an RFBG sensor (Regenerated Fiber Bragg 
Grating). Both sensors were produced in Standard 
Single-Mode Fiber G-652, of which the RFBG 
sensor was produced in hydrogenated fiber and 
subjected to heat treatment for regeneration. The 
two sensors were encapsulated with a stainless 
steel capillary tube (AlSI 304), temperature 
calibrated in the range between 5 “C and 60 ºC. 
The respective sensors were fixed in a metal 
sheet (wall), to measure temperatures in this 
range of operation, in order to compare durability 


Capítulo 3 


(aging), as well as precision and accuracy in relation to a K-type thermocouple. At the end 
of 16 weeks, no significant degradations were observed in the signal of each sensor and the 
measurement error was < 2 ºC. 

KEYWORDS: FBG; RFBG; Bragg grating durability; Encapsulation. 


11 INTRODUÇÃO 


As redes de Bragg em fibra ótica (FBG) foram descobertas em 1978 por K. O. Hill 
(HILL et al., 1978). Desde sua descoberta muitos trabalhos foram desenvolvidos em torno 
desses sensores óticos (OLIVEIRA, 2012; COSTA et al., 2015; MOURA et al., 2017; LU et 
al., 2018; PENG et al., 2019; FENG et al., 2020), devido as suas propriedades intrínsecas 
de medição de temperatura e deformação, ao mesmo tempo que são imunes a interferência 
eletromagnética (FRAZÃO, et al., 1999). Algumas aplicações desse tipo de sensor para 
medição de temperatura podem ser vistas na indústria petroquímica, pois tais sensores são 
aptos a operar em áreas classificadas já que não utilizam corrente elétrica (TERADA, 2010). 
Na biomecânica, os sensores baseados em redes de Bragg em fibra ótica utilizam-se da 
deformação para mensurar, por exemplo, força de mordida, aquisição de dados na aplicação 
de aparelhos ortodônticos e monitoração de problemas no maxilar (FIORIN et al., 2019). 
Também há estudos sobre o envelhecimento de sensores baseados em FBG, para avaliação 
de durabilidade e desempenho, por exemplo, no monitoramento de umidade em um ambiente 
na presença de gases residuais (ALWIS et al., 2017) e no monitoramento de vibração quando 
tais sensores foram encapsulados em metal (GUO; XIONG; LIU, 2019). 

O método de inscrição sob máscara de fase por uma fonte laser ultravioleta é 
amplamente utilizado para gravação de redes de Bragg em fibra Ótica (HILL et al., 1993), 
porque o método apresenta maiores estabilidade e repetibilidade da gravação quando 
comparado a outros métodos interferométricos (FOKINE, 2002; PAL et al., 2003; MIHAILOV 
et al., 2004). Para a gravação de redes sementes saturadas, para posterior processo de 
regeneração (OLIVEIRA, 2012), é necessária a etapa de hidrogenação da fibra ótica seguida 
da gravação da FBG. A hidrogenação aumenta a fotossensibilidade na fibra (PATRICK et al., 
1995; EGAN et al., 1996) e uma maior fluência laser na banda do ultravioleta produz redes 
sementes saturadas, as quais são caracterizadas por espectro alargado. A regeneração pode 
ser obtida através de tratamento térmico adequado da fibra ótica, ao final do qual é produzida 
a rede de Bragg regenerada em fibra ótica (RFBG). A RFBG é caracterizada por elevada 
durabilidade em altas temperaturas (CANNING et al., 2009). 

Os objetivos do presente trabalho são comparar a durabilidade entre dois sensores 
por rede de Bragg, um do tipo RFBG e um outro do tipo FBG padrão, durante um período de 
16 semanas em temperatura ambiente e também avaliar a precisão e a exatidão dos referidos 
sensores em relação a um padrão fornecido por termopar tipo K. Ao longo do período, foram 
realizadas medições semanais da temperatura obtida no termopar e em ambos os sensores 


através da aquisição de espectros em reflexão. 
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21 METODOLOGIA 


Os sensores térmicos foram produzidos em segmentos de fibra Ótica, monomodo 
padrão (SSMF G-652 Draka8), com -33 cm de comprimento. O sensor 1 foi gravado em 
fibra Ótica não hidrogenada, seguido de encapsulamento, calibração em temperatura e 
aplicação em campo. O sensor 2 teve sua rede semente gravada em fibra ótica hidrogenada, 
com posterior tratamento térmico para regeneração, também seguido de encapsulamento, 
calibração em temperatura e aplicação em campo. 

2.1 HIDROGENAÇÃO 


A câmara de hidrogenação é formada por tubo inoxidável, de 5 mm de diâmetro 
interno e 1 m de comprimento, comportando várias dezenas de fibras. Assim, um trecho de 
fibra Ótica monomodo padrão (G-652), com 1 m de comprimento, foi inserido na câmara de 
hidrogenação sob pressão de 80 bar, temperatura ambiente. A fibra Ótica permaneceu na 
câmara por um período de sete dias e, então, foi retirada para se obter um segmento de fibra 
necessário para a posterior gravação da rede semente do sensor 2. 

2.2 GRAVAÇÃO 


Ambos os sensores foram gravados com 3 mm de comprimento, utilizando o método 
de gravação sob máscara de fase. O laser utilizado nas gravações foi o Coherent? Xantos 
XS, com comprimento de onda de emissão em 193 nm. 

O sensor 1 (FBG), possui os seguintes parâmetros de gravação: energia de 1,80 mJ/ 
pulso, frequência de 250 Hz e tempo de exposição de 20 minutos. A máscara de fase utilizada 
foi a lbsenê com pitch de 1055,2 nm. 

O sensor 2 (RFBG), possui os seguintes parâmetros de gravação: energia de 1,00 mJ/ 
pulso, frequência de 250 Hz e tempo de exposição de 10 minutos. A máscara de fase utilizada 
foi a lbsen€ com pitch de 1058,9 nm. Como a fibra ótica utilizada para a gravação do sensor 
2 estava hidrogenada, foi possível utilizar menor fluência /aser neste caso. A rede produzida 
resultou em eficiente regeneração após o tratamento térmico específico. 

Os diferentes pitchs das máscaras de fase resultam em distintos comprimentos de 
onda de Bragg em cada sensor, porém não produzem efeitos apreciáveis na durabilidade e 
na precisão dos mesmos. 

Na Figura 1, pode ser visto um esquemático do sistema de gravação utilizado no 
Laboratório Multiusuário de Fotônica, na Universidade Tecnológica Federal do Paraná, 
campus Curitiba. Para a monitoração das gravações foi utilizado o interrogador ótico (Micron 
Opticsê SM 125) com resolução de 0,001 nm. 
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Figura 1 — Esquemático do sistema de gravação de FBG. 


Fonte: Adaptado de Oliveira (2012). 


2.3 TRATAMENTO TÉRMICO PARA REGENERAÇÃO 


A rede semente, que produzira o sensor 2, foi submetida a um tratamento térmico para 
obtenção da rede de Bragg regenerada. O referido sensor foi inserido em um forno elétrico 
tubular, que possui controle de temperatura até 1000 “C e foi desenvolvido no próprio grupo 
de pesquisa. O forno foi programado, para seguir uma rampa de aquecimento, iniciando 
em 22 ºC e atingindo 900 “C em 30 minutos. O forno permaneceu em 900 “C por mais 80 
minutos, até concluir a regeneração térmica da RFBG, resultando no sensor 2, seguido de 
resfriamento natural pela troca de calor com o ambiente. 

O setup de monitoração durante o tratamento térmico do sensor 2 constituiu-se de 
interrogador ótico, forno elétrico tubular e controlador de temperatura, conforme a Figura 2. 











1. Notebook 4. Interrogador VA 7. Forno 
2. Cabo de rede 5. RFBG 8. Termopar tipo K 
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Figura 2 — Setup para monitoração do tratamento térmico do sensor 2. 


O forno elétrico tubular é formado por uma resistência de níquel-cromo, 


posicionada a 15 cm da entrada, envolvendo um tubo em alumina com 12 mm de 
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diâmetro. A potência é fornecida por um relé eletrônico comandado pelo controlador Novus 
1200, conforme mostra a Figura 3. A rede de Bragg foi posicionada no centro da área de 
aquecimento, que possui comprimento de -3 cm. Para o monitoramento da temperatura 
interna do forno utilizou-se termopar tipo K ligado ao controlador, obtendo a temperatura na 
região próxima à FBG. 


Forno 





Tubo de 
alumina 









sensor 
sem encapsulamento 







[ aos E 
Scm Área de 
aquecimento 


Termopar K 


Controlador 
Figura 3 — Detalhes da parte interna do forno. 


Fonte: Adaptado de Oliveira (2012). 


2.4 ENCAPSULAMENTO 


Os sensores 1 e 2 foram encapsulados utilizando um tubo capilar em aço inoxidável 
(AISI 304), com diâmetro interno de 0,3 mm e espessura de parede de 0,6 mm. Esta proteção 
é necessária para aplicações em campo, pois as fibras Óticas são razoavelmente sensíveis 
a impacto e curvatura. 

A fibra Ótica possui uma proteção primária em acrilato, e outras camadas de 
materiais plásticos, que foram removidas para fabricação desses sensores e substituídas 
pelo encapsulamento em aço inox, conforme a Figura 4. A fixação entre o tubo capilar e a 
capa plástica externa do cordão ótico foi obtida utilizando tubo termorretrátil do mesmo tipo 
utilizado em isolamento de cabo elétrico. 
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Figura 4 — Esquemático do encapsulamento da FBG para medição em reflexão. 


Fonte: Adaptado de Oliveira (2012). 


2.5 CALIBRAÇÃO 


Os sensores 1 e 2 foram calibrados para obtenção da curva de sensibilidade térmica. 
Para isso, utilizou-se uma pastilha termoelétrica Peltier (TEC) e controlador específico, 
com as temperaturas ajustadas para 5 “ºC, 10 “ºC, 20 “C, 30 ºC, 40 ºC, 50 “Ce 60 ºC, pois 
os referidos sensores serão operados em temperatura ambiente. Para cada temperatura 
selecionada foi aguardado um tempo de 5 minutos para adequada estabilização e adquiridas 
as medidas, de 1 em 1 minuto, totalizando 11 aquisições para cada temperatura. Para cada 
nível de temperatura foram realizadas as médias para os comprimentos de onda de Bragg, 
resultando na obtenção da curva de sensibilidade térmica. 


2.6 APLICAÇÃO EM CAMPO 


Os sensores foram fixados em chapa metálica (parede), próxima a uma janela do 
laboratório, com o auxílio de grampos metálicos parafusados. O laboratório é condicionado em 
22 ºC +2ºC, porém, como a chapa metálica tem contado direto com o exterior, as variações 
de temperatura são apreciáveis. O propósito para os sensores foi medir a temperatura na 
chapa metálica, correlacionar com as medições de um termopar tipo K e também observar 
possíveis degradações na amplitude espectral (envelhecimento). 

O esquema que representa as ligações entre os equipamentos, os componentes 
elétricos e os sensores pode ser visto na Figura 5. Para evitar a fadiga ou deficiência de 
acoplamento nos conectores óticos, optou-se por seleção via chave ótica do sensor a ser 
medido em cada operação. A chave ótica (JDS UniphaseO 1:2 channel) é ativada por tensão 
continua de 5 V, que altera o caminho ótico entre C e porta 2 (leitura do sensor 2), para C e 
porta 1 (leitura do sensor 1). Como referência de temperatura, utilizou-se multímetro (ICELO 
MD-6111, com resolução de 1 ºC) em conjunto com termopar tipo K. Os dois sensores e o 


termopar foram acoplados termicamente à chapa metálica com o auxílio de pasta térmica. 
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Figura 5 — Equipamentos, componentes elétricos e sensores para medição da temperatura ambiente. 


Fonte: Adaptado de Oliveira (2012). 





31 RESULTADOS 


Conforme a Figura 6, para o sensor 2 foram obtidos o espectro inicial (rede 
semente), em linha contínua, e os espectros finais após a regeneração, em linha 
tracejada representando o espectro da RFBG sem encapsulamento e em linha traço-ponto 
representando o espectro da mesma RFBG com encapsulamento. Os espectros foram 
obtidos em temperatura de 22 ºC. É possível verificar a perda da saturação da rede de 
Bragg após a regeneração e seu deslocamento para um comprimento de onda menor. O 
encapsulamento não provocou alteração espectral significativa, indicando que a rede foi 
mantida sem deformação. 
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Figura 6 — Espectros em reflexão da FBG inicial (rede semente) comparada aos espectros regenerados 
após o tratamento térmico (RFBG com e sem encapsulamento). 


Durante o tratamento térmico para obter o sensor 2, foram adquiridos os espectros, 
com posterior leitura de suas amplitudes, utilizando software comercial, resultando no 
gráfico mostrado na Figura 7. A temperatura (triângulos) foi elevada de 22 “C a 900 “C em 
30 minutos e foi mantida em 900 ºC por mais 80 minutos, totalizando aproximadamente 
110 minutos. A amplitude do sinal refletido (círculos) permaneceu próxima de 26,4 dB até a 
temperatura atingir -850 “ºC, porém a rede perdeu a saturação nessa fase. O processo foi 
seguido de redução gradativa na amplitude do sinal, atingindo um mínimo próximo de 8,5 
dB aos 47 minutos. Por último, ocorreu a fase da regeneração (aumento gradativo do sinal) 
estabilizando a amplitude em aproximadamente 20 dB. 
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Figura 7 — Amplitude do sinal refletido comparada à temperatura do sensor 2. 
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Para a calibração dos sensores 1 e 2, foi utilizado uma pastilha termoelétrica Peltier 
(TEC). O sensor 1 apresentou sensibilidade térmica de 12,23 nm/ºC e o sensor 2 apresentou 


sensibilidade térmica de 11,80 nm/ºC, obtidas através de ajuste linear (Figura 8). 


Sensor 1 Sensor 2 
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Figura 8 — Curva de sensibilidade térmica (a) do sensor 1 e (b) do sensor 2. 


Durante 16 semanas, ambos os sensores monitoraram a temperatura ambiente (Figura 
10). Na Figura 10.a pode-se observar, nos espectros do sensor 1, insignificantes variações 
nas intensidades dos sinais refletidos e os pequenos desvios no comprimento de onda são 
devido às diferentes temperaturas durante tais aquisições. Na semana 7, foi observado 
um deslocamento na linha de base, o que sugere ocorrência de flutuação na fonte ótica do 
interrogador. Da mesma forma, isso ocorre com o sensor 2, conforme mostra na Figura 10.b. 
Ainda na Figura 10.b, observa-se razoável correlação com os espectros mostrados na Figura 
10.a, semana a semana, sugerindo alguma flutuação no sistema de interrogação/leitura e de 


temperatura entre as medições. 
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Figura 10 — Espectros em reflexão nas semanas 0, 7 e 15 (a) para o sensor 1 e (b) para o sensor 2. 
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A Figura 12 mostra as amplitudes dos sinais refletidos por ambos os sensores, ao 
longo das 16 semanas. Observa-se um padrão nas flutuações das amplitudes dos sinais 
refletidos pelos sensores, corroborando com uma provável flutuação na potência da fonte 
ótica do interrogador, como também pode ser observado nos espectros correspondentes 
as semanas 0, 7 e 15, mostrados na Figura 10. Aparentemente não ocorreu degradação 


apreciável em ambos os sensores durante o período observado. 
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Figura 12 — Variações nas amplitudes dos sinais refletidos, correspondentes ao sensor 1 e ao sensor 2, 
ao longo das 16 semanas. 


A Tabela 1 correlaciona as temperaturas obtidas pelos sensores 1 e 2 com a 
temperatura de referência, indicada no multímetro (termopar tipo K), ao longo do experimento. 
Para obter a temperatura através dos sensores 1 e 2, comparou-se o comprimento de onda 
de Bragg nos espectros adquiridos a cada semana com as respectivas curvas de calibração. 
Não foi observada divergência superior a 2 ºC em nenhuma das medições. 
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Temperatura Temperatura Erro absoluto Temperatura Erro absoluto 


Semana no multímetro no sensor 1  dosensor 1 no sensor 2 do sensor 2 





[ºC] [C] [ºC] [ºC] [ºC] 

0 27 27 0 27 0 
1 19 21 -2 21 -2 
2 27 Rr 0 27 0 
3 21 22 1 23 -2 
4 26 26 0 26 

5 27 27 0 27 

6 20 21 1 21 4 
7 25 25 0 26 1 
8 24 24 0 Bo 1 
9 24 25 1 26 -2 
10 23 23 0 24 1 
11 23 2 0 24 1 
tz 22 22 0 23 1 
13 23 23 0 24 1 
14 28 rá 1 29 1 
15 33 Er 1 33 0 





Tabela 1 — Temperaturas do multímetro comparadas às temperaturas dos sensores ao longo das 16 
semanas. 


41 CONCLUSÃO 


As leituras obtidas pelos sensores 1 e 2 referenciadas ao multímetro, apresentaram 
diferenças máximas de 2 ºC, eventualmente estando dentro da faixa de erro dos instrumentos. 
O experimento relatou 16 semanas e nesse tempo não foi possível mensurar perda de 
intensidade do sinal refletido pelos sensores óticos, até pelo fato de ter ocorrido alguma 
flutuação no próprio sistema interrogador, havendo uma variação de -2 dB entre os dois 
sensores, o que não influencia a análise dos resultados para os objetivos propostos. 

Devido a razoavelmente baixa temperatura de operação, que foi inferior a 40 ºC, a 
durabilidade e exatidão observada nos dois sensores óticos são razoavelmente próximas, 
o que eventualmente poderá mudar com o transcorrer do experimento, devido ao sensor 
2 ter sido produzido através de regeneração em FBG, normalmente mais estável. Ambos 
os sensores se mostraram adequados para aplicação em campo e apresentaram razoável 
precisão e exatidão em relação ao medidor de temperatura de referência. O encapsulamento 
proposto revelou-se eficiente, tanto pela robustez, quanto pela baixa influência sobre as 
características de uma FBG sem cápsula. 

Os dois sensores apresentados neste trabalho continuam sob análise na mesma 
instalação, possibilitando novas conclusões em médio/longo prazo. 
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RESUMEN: A inícios de la expansión de los 
Sistemas Eléctricos de Potencia se publicaron 
varias teorias respecto a las condiciones óptimas 
para transmitir la energia eléctrica. Una de ellas 
es la enunciada por Lord Kelvin, en Inglaterra 
en el afo 1881, teoria conocida como la Ley de 
Kelvin. En la actualidad los niveles de tensiones 
de transmisión y distribución están definidas por 
normas internacionales y regionales, en este 
estudio se analizan los paradigmas relacionados 
a que valor de estas tensiones normalizadas 
sería el óptimo para el caso de estudio y para ello 
se aplica el fundamento de la Ley de Kelvin, con 
el fin de evaluar inversiones para la expansion 
del sistema eléctrico de distribución. Se debe 
mencionar que en la práctica este concepto tiene 
limitaciones, por lo que en este trabajo se disefia 
una metodologia utilizando algoritmos en donde 
se integran los valores anuales de las variables 
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que intervienen en el fundamento teórico, tambien 
se definien las restricciones de tal forma a que 
los resultados optenidos sean aplicables. Así, 
gracias a la metodologia desarrollada, se hallan 
las curvas de costos anuales de construccion, 
depreciacion, intereses para el financiamiento 
de las líneas y perdidas de energia eléctrica, 
donde la función objetivo será minimizar la suma 
de los mismos. Para ello se utilizan funciones de 
comandos de Matlab, obteniendo de esta forma 
la sección óptima del conductor como asi también 
la tensión de distribución óptima asociado a este 
para una proyección de 10 afos. 

PALABRAS CLAVE: Ley de Kelvin. Cálculo 
de Pérdidas de Energia, Costos Anuales 
de construcción de líneas de distribución. 
Optimización estática. 


AVALIAÇÃO DE INVESTIMENTO 
COM OTIMIZAÇÃO DE SECÇÃO 
DO CONDUTOR E TENSÃO DE 
DISTRIBUIÇÃO. APLICAÇÃO DOS 
ALGORITMOS DA LEI DE KELVIN 


RESUMO: No início da expansão dos Sistemas 
Elétricos de Potência, várias teorias foram 
publicadas sobre as condições ótimas de 
transmissão de energia elétrica. Uma delas é a 
enunciada por William Thomson, na Inglaterra 
em 1881, teoria conhecida como Lei de Kelvin. 
Atualmente os níveis de tensões de transmissão 
e distribuição são definidos por normas 
internacionais e regionais, neste estudo são 
analisados os paradigmas relacionados a qual 
valor dessas tensões normalizadas seria o Ótimo 
para o estudo de caso e para isso são analisados 
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os fundamentos da Lei de Kelvin, em a fim de avaliar investimentos para a expansão do 
sistema de distribuição de energia elétrica. De referir que na prática este conceito apresenta 
limitações, pelo que neste trabalho é desenhada uma metodologia utilizando algoritmos onde 
se integram os valores anuais das variáveis que intervêm na fundamentação teórica, sendo 
as restrições também definidas de forma a que os resultados obtidos são aplicáveis. Assim, 
graças à metodologia desenvolvida, encontram-se as curvas de custos anuais de construção, 
depreciação, juros de financiamento das linhas e perdas de energia elétrica, onde a função 
objetivo será minimizar a soma delas utilizando funções de comando do MATLABO; obtendo 
assim a seção ótima do condutor, bem como a tensão de distribuição Ótima associada a ele 
para uma projeção de 10 anos. 

PALAVRAS-CHAVE: Lei de Kelvin. Cálculo de perdas de energia. Custos anuais de 
construção de linhas de distribuição. Otimização estática. 


EVALUATION FOR INVESTMENT WITH OPTIMIZATION OF CONDUCTOR 
SECTION AND DISTRIBUTION VOLTAGE. APPLICATION OF KELVIN'S LAW 
ALGORITHMS 


ABSTRACT: At the beginning of the expansion of Electric Power Systems, several theories 
were published about the optimal conditions for electric energy transmission. One of them 
is enunciated by William Thomson, in England in 1881, a theory known as Kelvin's Law. 
Currently, the transmission and distribution voltage levels are defined by international and 
regional standards, in this study the paradigms related to which value of these normalized 
voltages would be the optimal for the case study are analyzed, and for this, the fundamentals 
of Kelvin's Law are analyzed, in order to evaluate investmentis for the expansion of the electric 
energy distribution system. It should be noted that in practice this concept has limitations, 
so in this work a methodology is designed using algorithms where the annual values of the 
variables that intervene in the theoretical foundation are integrated, and the restrictions 
are also defined so that the results obtained are applicable. Thus, thanks to the developed 
methodology, the annual construction cost curves, depreciation, interest on financing lines 
and electricity losses are found, where the objective function will be to minimize their sum 
using MATLABO command functions; thus obtaining the optimum conductor section as well 
as the optimum distribution voltage associated with it for a 10-year projection. 

KEYWORKS: Kelvin's Law. Calculation of energy losses. Annual costs, construction of 
distribution lines. Static optimization. 


11 INTRODUCCIÓN 


En Paraguay existen varias restricciones para la expansión del sistema de distribución 
de energia eléctrica, principalmente en áreas urbanas densamente pobladas. Una de ellas 
representa a los problemas estructurales para la edificación de nuevas obras con el fin 
atender la demanda de suministro de energia eléctrica mediante la construcción de nuevas 
Estaciones de Transformación, como así también ampliaciones de la red de media tensión 
de distribución. 

Para el caso de la expansión del sistema de distribución, el medio físico donde se 
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deben hacer las ampliaciones de los locales de transformación, generalmente están muy 
limitados en cuanto al espacio para la construcción e instalación de ductos y salidas de nuevos 
alimentadores, esto hace que inevitablemente se tenga aglomeración de los conductores en 
las canaletas o ductos, lo que hace que la capacidad de conducción de corriente real de 
los mismos sea considerablemente menor, esto conlleva pérdidas considerables de energia 
por efecto Joule y en el peor de los casos, la ocurrencia de fallas en los conductores con la 
consecuencia de un corte de energia eléctrica, valorado a partir del Costo de Energia no 
Suministrada (CENS). 

Otra problemática desde el punto de vista estructural se observa en el sistema 
de transmisión, en donde se tiene la dificultad de emprender la construcción de nuevas 
Estaciones de Transformación en áreas urbanas densamente pobladas, debido a la escasez 
de terrenos de las dimensiones adecuadas para instalación de las mismas. 

Estas limitaciones, son en la actualidad algunas de las principales condicionantes 
para la expansión del Sistema Eléctrico de Distribución en Paraguay. 

En base al análisis de la problemática enunciada, a fin de optimizar las inversiones y 
de obtener mayor capacidad de transporte de energia eléctrica en el sistema de distribución, 
en este estudio se aplica la Ley de Kelvin. De esta forma se encuentra el momento y la 
sección de conductor económicamente conveniente de manera a incrementar la línea troncal 
del alimentador estudiado y seleccionado como muestra. 

Con este mismo fundamento, se realiza un análisis similar con otros valores de tensión 
nveces superior a la tensión de distribución actual, considerando que la relación de capacidad 
de transporte en Vatios a una tensión superior es n?, ya que con ello se tendria la posibilidad 
de transportar n?veces más potencia en las mismas condiciones de infraestructura en cuanto 
a espacio físico disponible. 

Para ello, en este trabajo se investigan los criterios actuales para definición de niveles 
de tensión de distribución y se desarrolla una metodologia para hallar la tensión óptima para 
centros de carga con las características del alimentador seleccionado como muestra. 

Todas las variables que intervienen en la aplicación de Ley de Kelvin son integradas 
al algoritmo desarrollado en este estudio. 

Se debe hacer mención que, en la práctica, el fundamento de la Ley de Kelvin tiene 
ciertas limitaciones para su aplicación, la cuales se citan a continuación: 

| No es fácil estimar la pérdida de energia en la línea sin curvas de carga 
reales, que no están disponibles en el tiempo de estimación. 


Il. El costo anual a cuenta de intereses y la depreciación en el desembolso de 
capital no tienen una forma lineal. 


Il. Esta ley no tiene en cuenta varios factores como la densidad de corriente 
segura, resistencia mecánica, pérdida de corona, etc. 


IV. El tamafio del conductor determinado por esta ley puede no siempre puede 
ser practicable, puede ser demasiado pequehio para el transporte seguro de 
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la corriente necesaria. 


V. El interés y la depreciación en el desembolso de capital no pueden ser 
determinados con precisión. [1]. 

Con la metodologia a implementar se lograrã gracias a la programación y 
manipulación de datos masivos, estimar las curvas de cargas reales determinando 
posteriormente las perdidas por efecto Joule, utilizando el Software CYMDIST '! y teniendo 
en cuenta diferentes valores de proyección de la demanda de energia eléctrica. También 
se obtienen las curvas referentes a los costos de los materiales y la variación de estos en 
función a la sección y tensiones de servicio de los mismos. 

Para que el conductor a definir, utilizando el concepto de la Ley de Kelvin, sea 
seguro (referida a la densidad de corriente), se utilizará en el algoritmo el Costo de la 
Energia Eléctrica no Suministrada (CENS), de tal forma a que este valor sea la restricción 
a las soluciones posibles. 

De esta manera, una vez obtenidas las matrices de estas variables la función 
objetivo será minimizar la suma de la relación entre el costo anual de interés y depreciación 
referente a los costos de construcción y el área de la sección del conductor, como así 
también la relación entre el costo anual de energia desperdiciada y el área de sección del 
conductor, así se obtendrá la sección óptima económica para la proyección de la demanda 
especificada. 

La misma metodologia será aplicada a otros niveles de tensión normalizados a fin de 
verificar la variación de los costos y la conveniencia respecto a la construcción de nuevas 
Estaciones de Transformación. 


OBJETIVOS DEL TRABAJO 


Tomando como base lo expuesto en el apartado anterior, en el presente trabajo se 
propone como objetivo general “Hallar el valor óptimo de tensión de distribución, como así 
también la sección de conductor óptima para el alimentador seleccionado como muestra, 
desarrollando una metodologia basada en la Ley de Kelvin, con la cual se pueda evaluar 
opciones de inversión para un plazo de 10 afos”. 

Considerando la necesidad de que el sistema de distribución se adapte a la expansión 


del sistema eléctrico, además se instauran los siguientes objetivos específicos: 


* Identificar los criterios y normas a ser tenidos en cuenta para adoptar un deter- 
minado nivel de tensión en distribución. 


* | Desarrollar algoritmos a fin de observar el comportamiento de las variables y 
realizar los cálculos necesarios para la aplicación de la Ley de Kelvin para dis- 
tintos escenarios. 


* | Estimar los costos del financiamiento para la construcción y depreciación de 


1 Paquete base para el análisis del sistema de distribución del software CYME, utilizada en la Administración Nacional 
de Electricidad (ANDE) 
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nuevas líneas de distribución. 


* | Hallar los costos anuales de las pérdidas de energia eléctrica por efecto Jou- 
le en conductores a tensiones normalizadas para el alimentador seleccionado 
como muestra. 

Las preguntas guía para el presente trabajo son: i) Ya que en el segmento de 
generación y en muchos casos en transmisión, se encuentran cantidades de excedentes 
de potencia. ; Qué alternativa sería la óptima para distribuir dicha cantidad de excedente 
teniendo en cuenta las restricciones estructurales existentes? é; Cómo se pueden integrar 
las distintas variables que componen la determinación de un valor de tensión y sección de 
conductor óptimo en distribución? En relación con las preguntas se plantean las siguientes 
hipótesis: |) El valor óptimo de tensión y sección óptima del sistema de distribución en 
estudio es el que presenta menores costos anuales totales para el tiempo analizado. ii) El 
desarrollo de un algoritmo basado en la Ley de Kelvin permitirá evaluar los costos anuales 
totales para diferentes niveles de tensión y secciones de conductores. 


REVISIÓN DE LA LITERATURA 
Las Tensiones de Distribución Normalizadas en Paraguay se definen en el Pliego de 
Tarifas Nº21 [2]. 
Tensiones trifásicas (entre fases) 
Muy alta tensión: 220.000 Voltios 
Alta tensión: 66.000 Voltios 
Media tensión: 23.000 Voltios 
Baja tensión: 380/220 Voltios 
Tensiones monofásicas (entre fase y neutro) 
Media tensión: 23.000 / V3 Voltios 


Este pliego a su vez está basado en el Reglamento Para Instalaciones Eléctricas de 
Media Tensión, que data del afio 1975, así como también del Reglamento para instalaciones 
de Baja Tensión que data del aho 1971 [2]. 

En los últimos ahios se construyeron las líneas de transmisión ITAIPU-VILLA HAYES 
y AYOLAS-VILLA HAYES, en 500kV, por lo que este nivel de tensión pasa a formar parte de 
los niveles de tensiones utilizados en el Sistema Eléctrico Paraguayo. 

Los conductores normalizados y utilizados en el sistema de distribución eléctrica en 


Paraguay y a los que está enfocado el presente estudio son: [5]. 
Conductor Desnudo de Aleación de Aluminio 


* Conductor Desnudo de Aleación de Aluminio, de 35 mm? 


* Conductor Desnudo de Aleación de Aluminio, de 50 mm? 
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* Conductor Desnudo de Aleación de Aluminio, de 70 mm? 

* Conductor Desnudo de Aleación de Aluminio, de 95 mm? 

* Conductor Desnudo de Aleación de Aluminio, de 150 mm? 

* (Conductor Desnudo de Aleación de Aluminio, de 240 mm? 

* (Conductor Desnudo de Aleación de Aluminio, de 300 mm? 
Normas internacionales 


Definición de los niveles de tensión según la Norma IEC 60038 


Tensión máxima Tensión nominal del Tensión máxima Tensión nominal del 
para el equipo (kV) sistema (kV) para el equipo (kV) sistema (kV) 





Tabla 1- Tensiones nominales entre 1 kV y que no excedan de 36 kV. [6]. 


Ley de Kelvin 


El área de conductor más económica es aquella para la cual el costo total anual de la 
linea de transmisión es mínimo. Esto se conoce como Ley de Kelvin. Fue Lord Kelvin quien 
lo declaró por primera vez en 1881. “El costo anual total de la línea de transmisión se puede 
dividir ampliamente en dos partes: Tasa anual de desembolso de capital y costo anual de 
energia desperdiciada en el conductor” [7]. Si bien la definición de la Ley de Kelvin se refiere 
a sistemas de transmisión, para este estudio, este concepto será aplicado a sistemas de 
distribución *. 

Este costo anual se puede ilustrar gráficamente trazando el costo anual contra el 
àrea 'a' del conductor como se muestra en la Ilustración 1. En el diagrama, la línea recta (1) 
muestra la relación entre los costos anuales (es decir, P1 + P2a) y el área de la sección 'a' 


2 Elaborado por el autor 
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del conductor. De manera similar, la hipérbola rectangular (2) proporciona la relación entre el 
costo anual de energia desperdiciada y el área de sección 'a'. Al sumar las ordenadas de las 
curvas (1) y (2), se obtiene la curva (3). Esta última curva muestra la relación entre el costo 
anual total (P1 + P2a + P3/a) de la línea y el área de la sección 'a”. El punto más bajo de la 


curva (es decir, el punto P) representa el área más económica de la sección. [1]. 





mosrimcr 


E) 





Costos anuales 


Ilustración 1 - Representación gráfica de la Ley de Kelvin [1]. 


a: Sección del conductor. 

P|: Costo de la línea de transmisión o distribución que no depende de la sección del 
conductor. 

P,: Costo de la línea de transmisión o distribución que depende de la sección del 
conductor (Interés para el financiamiento de la construcción de la línea y depreciación del 
material). 

P,: Costo de las pérdidas por efecto Joule. 

Para la elección del voltaje óptimo, se obtiene una curva para el costo total de 
transmisión contra el voltaje, como se muestra en la Ilustración 2. El punto más bajo (P) en la 
curva proporciona la tensión de transmisión económica. Por lo tanto, en el presente caso, OA 


es la tensión de transmisión óptima [1]. 
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Costos de capital 





O A 


Voltaje de Transmisión 


Ilustración 2 - Voltaje de transmisión óptimo [1]. 


Por otro lado, la tendencia actual es seguir ciertas fórmulas empíricas para encontrar 
la tensión de transmisión económica. Así, de acuerdo con la práctica estadounidense, el 
voltaje económico entre las líneas de un sistema trifásico AC es: [7]. 

V=5,5v0,62+ 2 
150 (1) 
|: longitud de la línea de transmisión o distribución. 

P: Potencia a transmitir. 

V: Tensión económica 


Proyección de la Demanda 


En la Tabla 2, se presenta una comparación de las tasas de crecimiento de la 
demanda máxima consideradas por la Administración Nacional de Electricidad (ANDE) en 
su plan de expansión y de las resultantes del análisis de demanda Ilevada a cabo en el 


marco del estudio de tarifas eléctricas. 


= , : Caso Caso 


Tabla 2 - Comparación de Demanda máxima - tasas anuales acumulativas proyectadas 2012-2022 [8]. 





3 Estudios Energéticos Consultores 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 4 ES 


Costo Marginal de Generación 


En función a estudios referenciados, se definió como Costo Marginal de Largo Plazo 
de la Energia (CMLPE) de 35,66 USD/MWh. Que es el costo a utilizar para hallar los costos 
por perdidas de energia [8]. 


Determinación del Costo de la Falla 


Los resultados de estudios referenciados oscilan en el rango [2,94 USD/kWh a 4,09 
USD/kWhl] los cuales resultan consistentes con las experiencias internacionales [8]. 


MARCO METODOLÓGICO 


Universo de estudio 


Los análisis se basaron en la Estación Transformadora de 220/66/23 kV Ilamada 
Estación San Lorenzo, en adelante ES-SLO, y como muestra se ha seleccionado el 
alimentador de 23 KV SLO-18. 


Técnicas de recolección de datos 


La selección de nivel de tensión de distribución para zonas urbanas será hecha con 
la técnica del análisis documental. Así, las tensiones normalizadas serân las expuestas en la 
Tabla 1 del presente estudio. 

Con respecto a la definición de los Costos de construcción de líneas de distribución 
para diferentes secciones de conductores normalizados, será hecha con la técnica del 
análisis documental. De esta manera se tienen los costos modulares para la construcción 
de los alimentadores aéreos de Aleación de Aluminio Desnudo de Media Tensión (MT). En la 


Tabla 3 son representados estos valores en USD/km. 


Materiales 
Construcción de Línea de MT Mano de Total 
Estructuras/ obra 
Conductor 
postes 


Línea Trifásica de 3x35mm? 9.199 3.565 13.248 
Línea Trifásica de 3x5Omm? 10.119 3.921 14.573 


Línea Trifásica de 3x150mm? 15.981 6.906 23.728 
Línea Trifásica de 3x240mm? 1.093 20.775 7.485 29.354 
7 Línea Trifásica de 3x300mm? 1.421 27.008 8.113 36.542 


Tabla 3 - Costos modulares de materiales y mano de obra (USD/km), para construcción de 
alimentadores MT [5]. 





Línea Trifásica de 3x95mm? 15.319 6.371 22.496 
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Interés y depreciación de líneas de distribución 


Los valores hallados correspondientes a la Depreciación Anual para los diferentes 


tipos de líneas adoptadas en este estudio fueron: 


ce23= [68073 74880 87198 113357 118259 153736 199857] 
Ccc23= [3583 3941 4589 5966 6224 8091 10519] 
p23= ce23 + cc23 


23—-vd 
d23 = PT (2) 


vu 
ce23 = costos de postación (estructuras, aisladores, postes), tensión 23kV 
cc23 = costo de conductores, tensión 23kV 
cil = (cons23 — d+ a) xl (3) 
a = amortización 
vu = Vida Útil 
vd = Valor de desecho (se considera 0) 


Los intereses para el financiamiento de las obras para los conductores se hallaron de 


la siguiente forma: 


cons 23 = cc 23 + ce 23 + cm 23 = costo de la obra (USD) (4) 

cí 23 = cís + cíl = costo anual de interés del financiamiento (USD) (5) 
a = cons 23/n (6) 

cís = (cons 23 — 2,” a) * ís (7) 


n = número de periodos de pagos 

is = Interés fijo 

| = Interés variable (tasa Libor [10]) 

ce23 = costos de postación (estructuras, aisladores, postes), tensión 23kV 

cc23 = costo de conductores, tensión 23kV 

cm283 = costo de mano de obra, tensión 23kV 

Una vez obtenidos los valores anualizados de los intereses y depreciación se definen 
los mismos, a fin de tener la curva 1 (Ilustración 1) a ser analizados con el fundamento de 
la ley de Kelvin. La matriz de resultados se muestra en la Tabla 4 con su correspondiente 


representación en la Ilustración 3. 
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DT Simm 
me/ ss [80 [10 [o [10 [20 | 0. 


Tabla 4 - Costos anuales (USD) de interés y depreciación - tensión 23kV. 





Cargo Anual de Intereses y Depreciación 






Referencia Sección |... 
(eex) . 














1,8 ee Do — OO do 
A s-—|2=S0mmê 
O rad-e|jo = FUMO 
=) 124 nte? 
o it 5 = 150mm? 
TD 08 6 = 240mm? 
O qual e 7 = 300mm? 
Õ E ane Í 


Afios Sección 


Ilustración 3 - Costos Anuales de Interés y Depreciación, tensión 23kV. 


Los resultados obtenidos indicados en la Tabla 4 son mostrados en función del tiempo 
y para cada sección de conductor estudiado. Se puede observar que para una sección 
menor se tiene menores costos anuales referentes a la depreciación y gastos en intereses 


en relación a los conductores de secciones superiores que requieren mayor inversión inicial. 
Pérdidas de energia eléctrica, para líneas MT de 23kV 


El Costo de la Energia no Suministrada (CENS), fue hecha con la técnica del análisis 
documental, se ha utilizado el valor de 4.09 USD/kWh, que es el resultado obtenido por [8]. 

Las pérdidas y costos por efecto Joule en las líneas se hallaron mediante simulaciones 
hechas en el software CYMDIST *. 

En la Tabla 5 son representados los valores de las pérdidas en los conductores para 
distintos factores de carga, calculadas en el software. 


4 CYMDIST: paquete base para el análisis del sistema de distribución del software CYME, utilizada en ANDE. 
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fm | me | tó | tome | tio | toe | ate | tó 
Carga | E 
(KW) ao érdidas (KW) 
2000 | 5565 | ao | ao [| 2260 | 1670 | 1000 | 66 | 55 — 


Tabla 5 - Pérdidas en los conductores por efecto Joule, tensión 23kV (Simulación efectuada en el 
software CYMDIST). 





A fin de tener cuantificadas y anualizadas las pérdidas se elabora una curva de 
pérdidas característica (ilustración 4) a partir de la Tabla 5 para cada sección de conductor 
y para diferentes valores cargas horarias (CAH), que fueron datos suministrados por el 
Centro de Operaciones de ANDE. 


Ecuaciones obtenidas a partir de la matriz de pérdidas 


p35 (i) = 0,00129*CAH (i) 2, 0208 (8) 
p50 (i) = 0,009 *CAH (i) 2, 0229 (9) 
p70 (i) = 0,0071*CAH (i) 2, 0035 (10) 
p95 (i) = 0,0053*CAH (i) 2, 0008 (11) 
p150 (i) = 0,0031*CAH (1)2,0122 (12) 
p240 (i) = 0,002*CAH (i) 2, 0177 (13) 
p300 (i) = 0,0018*CAH (i) 1,9982 (14) 


p (i) = Energia pérdida (kWh), para cada tipo de conductor. 
CAH = Carga máxima horaria del alimentador (A) SLO — 18. 
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Curvas de perdidas en conductores para el alimentador SLO-18 
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Ilustración 4 - Curva de pérdidas por sección de conductor. 


Principal del alimentador SLO-18. 


Cálculo del costo de pérdidas anuales en USD, por cada tipo de conductor 


Desarrollado el algoritmo, se aplica una tasa de crecimiento de la demanda de 9,3%, 
se obtiene la matriz C, ..o,10' AU contendrá la demanda horaria proyectada para el tiempo de 
estudio. Si bien es poco probable que la configuración del alimentador se mantenga invariable 
en el tiempo, para este estudio se considera constante. 

El algoritmo compara cada uno de estos valores con la capacidad de conducción de 
cada uno de los tipos de conductores principales estudiados, definiendo como restricción el 
límite de conducción de estos, para así aplicar el criterio del CENS. 

Queda así definida la matriz cp(i,j) , que es la matriz del costo anualizado de las 
pérdidas para cada tipo de conductor para un periodo de 10 afios, a ser analizadas con el 
fundamento de la ley de Kelvin y representados en la Ilustración 5 y en la Tabla 6 en forma 


matricial. 


DIO mm 
nu so | 70 | 95 | 150 | 240 | 30 | 
ER 
ES 
ES 


10.172.649 3.752.970 636.745 79.223 217.928 18.538 15.083 
13.770.329 5.995.614 1.018.798 201.971 35.412 22.181 18.016 


17.842.561 9.167.516 3.189.729 538.055 52.977 26.541 21.519 
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Tabla 6 - Costos anual de energia desperdiciada en conductores, CENS más Joule (USD), tensión 
23kV. 


Cargo anual de Pérdidas de Energia 
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Ilustración 5 - Costos anuales de energia desperdiciada en conductores. 


CENS más Joule (USD), tensión 23kV. 


Los resultados obtenidos indicados en la Tabla 6 son mostrados en función del tiempo 


y para cada sección de conductor estudiado. Se puede observar que para una sección menor 


se tiene mayores costos anuales, esto es debido principalmente a la aplicación del costo de 


la energia eléctrica no suministrada. 


RESULTADOS OBTENIDOS 


La metodologia desarrollada permite el cálculo de los costos anuales de líneas de 


distribución. Estos cálculos pueden ajustarse a los distintos escenarios a que se ve sometido 


la expansión del sistema eléctrico, como ser las variables referidas a los intereses para el 


financiamiento de las obras, la demanda de energia eléctrica asociada a la variación del 


PIB, la variación del costo de los materiales, etc. En este estudio se hicieron los cálculos 


con las condiciones más exigentes para el sistema, y los resultados son representados en la 


siguiente animación en 3D. 
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21 CONCLUSIONES 


Se identificaron las tensiones normalizadas según norma IEC 60038, esto es muy 
importante ya que la mayoria de los fabricantes de equipos que son utilizados en sistemas de 
distribución, los fabrican para estos valores. Por lo que los cálculos y estimaciones también 
fueron hechos para estos. 

Se definieron los costos de construcción de líineas de media tensión utilizando como 
fuente la base de datos de ANDE para líneas con conductores de hasta 150mm?, para los 
conductores de secciones superiores, se estimaron con factores proporcionales. Para los 
costos de construcción de líneas de media tensión a tensiones superiores se utilizaron 
ecuaciones y normas para la determinación de las distancias horizontales y verticales 
respectivamente, estimando la diferencia de costos de los materiales con respecto a las 
lineas de 23kV para la construcción de líneas de hasta 44kV. Para este estudio se consideró 
que la línea no tiene otra de menor tensión en el mismo recorrido, ya que si fuese así, las 
distancias eléctricas verticales serían mayores y necesariamente de deberian aumentar la 
altura de los postes. Para esta hipótesis, según norma, no fue necesario aumentar las alturas 
de los postes. 

Se logró, a través de entrevistas y revisión documental, conocer y tener en cuenta 
las condiciones contractuales para el financiamiento de las obras de expansión del sistema 
de distribución, así también, se estudió la metodologia y los valores para el cálculo de la 
depreciación de los bienes, necesarios para la aplicación de la ley Kelvin. 

Mediante el programa CYMDIST, se logró calcular las pérdidas para diferentes 
secciones de conductores para línea de media tensión de 23kV y obtener las ecuaciones 
de pérdidas en función de la carga, como así también para tensiones superiores (34,5kV), 
según norma IEC 60038, también se ha simulado para una tensión de 44kV, a pesar de 
no corresponder a la referida norma, debido a que es también una tensión de distribución 
utilizada en la región. Una vez obtenidas las pérdidas anuales para una proyección de la 
demanda de 10 ahos, se obtuvieron las matrices de pérdidas para cada sección de conductor 
y tensión de distribución. 

Halladas las pérdidas por efecto Joule en las líneas de distribución, aplicando el 
algoritmo desarrollado, se obtuvieron las matrices de costos de pérdidas anuales, para cada 
sección de conductor y tensiones de distribución. 

Se halló el conductor óptimo con el fundamento de la Ley de Kelvin, se obtuvieron los 
resultados para el escenario propuesto. Los cálculos se efectuaron en el escenario donde se 
adoptaron los valores más altos correspondientes a tasas de intereses y tasa de crecimiento 
de la demanda. Se obtuvieron las matrices de costos anuales para cada sección de conductor 
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y tensión de distribución, llegando a la conclusión de que no es conveniente aumentar el 
nivel de tensión de distribución aplicado a la muestra en estudio, debido a que, para estos 
casos, los costos globales a tensiones superiores son mayores a los costos para la tensión 
de distribución a 23KV. 

Con todos los objetivos específicos logrados, y superadas las limitaciones de la Ley 
de Kelvin, se concluye que se logró el objetivo general de este estudio, el cual es definir el 
valor óptimo de tensión de distribución para la muestra en estudio, siendo así la tensión de 
23kV, el valor óptimo para este caso. También se determinó la sección de conductor óptima, 
comprobando que con el algoritmo y la metodologia utilizada se puede definir el momento 
más conveniente de ampliación de la red de distribución estudiada, o aumento de sección del 


mismo teniendo en cuenta todas las variables que intervienen en la Ley de Kelvin. 
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ABSTRACT: The energy extracted from 
petroleum-derived fuels is becoming more 
expensive and of difficult availability. When these 
materials are transformed into electrical energy, 
they radiate heat energy and as a consequence 
of this combustion, pollutes the environment. On 
the other hand, having water resources is not 
easy when they are not available and cannot be 
processed in their natural form. It should be noted 
that the Electricity Generation from the Nuclear 
reaction of radioactive elements is very expensive. 
It is true that with a small percentage of atomic 
fuel, heat energy can be generated more than 
1000 times than the usual fuels. Consequently, 
the use of this kind of energy requires greater care 
in handling it. The radioactive waste generated by 
this atomic reaction cannot be destroyed, it can 
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SYSTEMS STABILITY 


only be stored away from natural life. There are 
many methods of creating Clean Electric Energy, 
from Renewable Resources such as water, wind, 
biomass, Solar Light, static energy stores, etc. 
These methods make the production of clean 
electrical energy moderately expensive. Unless 
there is a specific and mandatory regulation of 
a Central Government for the use of this type of 
Generation and interconnected to the National 
Electricity Network. This Electrical Interconnection 
should be done in an interesting way in terms 
of economics and the Stability of the Electrical 
Power System. In this research, a portion of an 
Electric Power Transmission Subsystem is used 
for the study; electrically located towards the 
North-West of the main “Electric Center” of the 
National Electricity Network. A photovoltaic plant 
of 5OMWPp of capacity in alternating current, and a 
Battery Energy Storage System of 15MW in direci 
current were introduced, where there are very 
long Transmission Lines, Electrical Substations 
with very radial loads. In this part of the electrical 
network, the Diesel Generation is usually a 
necessity, both to supply feeder loads as well as 
to control Voltage at times of maximum demand. 
This power system is normally connected to the 
National Electricity Network. 

KEYWORDS: Fuzzy logic, voltage control, 
artificial intelligence. 


Capítulo 5 


ESTRATÉGIAS DE CONTROLE DE TENSÃO BASEADAS EM ALGORITMOS 
DE LÓGICA FUZZY COM GERAÇÃO ALTERNATIVA, LIMPA E RENOVÁVEL 
OPERANDO COM OUTRA GERAÇÃO ELÉTRICA CONVENCIONAL COM 
CARGAS RADIAIS NA ESTABILIDADE DE SISTEMAS DE ENERGIA 


RESUMO: A energia extraída dos combustíveis derivados do petróleo estão-se tornando 
mais cara e de difícil disponibilidade. Quando esses materiais são transformados em energia 
elétrica, eles irradiam energia térmica e, como consequência dessa combustão, poluem 
o meio ambiente. Por outro lado, ter recursos hídricos não é fácil quando eles não estão 
disponíveis e não possam ser processados em sua forma natural. Deve-se notar que a 
Geração da Eletricidade a partir da reação nuclear de elementos radioativos é muito custosa. 
É verdade que com uma pequena porcentagem de combustível atômico, a energia térmica 
pode ser gerada mais de 1000 vezes do que os combustíveis usuais. Consequentemente, o 
uso desse tipo de energia requer maior cuidado em sua manobra. O lixo radioativo gerado 
por esta reação atômica não pode ser destruído, só pode ser armazenado longe da vida 
natural. Existem muitos métodos de geração de energia elétrica limpa, a partir de recursos 
renováveis, como água, vento, biomassa, luz solar, armazenamento de energia estática, 
etc. Esses métodos tornam a produção da energia elétrica limpa moderadamente custosa. 
Exceto regulamentação específica e obrigatória de Governo Central para a utilização desta 
modalidade de Geração e interligada à Rede Elétrica Nacional. Esta Interligação Elétrica deve 
ser feita de forma interessante em termos de Economia e Estabilidade do Sistema Elétrico 
de Potência. Nesta pesquisa, uma parte do Subsistema de Transmissão de Energia Elétrica 
é usada para o estudo; localizado eletricamente a Noroeste do principal “Centro Elétrico” da 
Rede Elétrica Nacional. Foi introduzida uma central fotovoltaica de 5OMWPp de capacidade 
em corrente alternada e um Sistema de Armazenamento de Energia a Baterias de 15MW em 
corrente contínua, onde existem Linhas de Transmissão muito longas, Subestações Elétricas 
com cargas muito radiais. Nesta parte da rede elétrica, a Geração Diesel costuma ser uma 
necessidade, tanto para suprir cargas dos alimentadores quanto para controlar a Tensão nos 
momentos de máxima demanda. Este sistema de energia é normalmente conectado à Rede 
Nacional de Eletricidade. 

PALAVRAS-CHAVE: Fuzzy Logic, controle da Tensão, Inteligência Artificial. 


1 | GENERAL INTRODUCTION 


In this part of the research, based on the previous studies, we will include the 
Bidirectionality of the equipment, the flywheels of inertia and we will analyze the results 
to advance towards the investigations of renewable energy (specifically Solar Light) with 
PVP. The main objective of this study is to determine the best management of the reactive 
power for optimal voltage control at strategic points where the stresses are being felt. The 
secondary objective is to use the resources available in the power system under study. 
All the signals that will go to the control sensors (primary and secondary voltage controls, 
generation prices) through a powerful and versatile artificial intelligence (Al) tool; such as 
complex fuzzy logic (FL) decision algorithms. 
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Electromechanical Characteristics 





220kV 
SE-LPA bus 1 
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Figure 1. PES in studies with injection NES, DG, PVP and BESS. 


ELECTRICAL SUBSTATION 


It remains practically the same, except for the 50MWp (AC) and 15MW (DC) PVP and 
BES [23] System respectively, connected directly to the SE-FIL 66kV bus load center. 


TRANSMISSION LINES 


In general terms, the Transmission Lines did not suffer great variations. 


GENERATION 


As for the Generation, it was fully adapted with GD. Electric Substation has its own DG, 
requesting to operate without interconnection with the NEN, to operate in Islands or to operate 
in fully autonomous and independent Sub-Islands. The Electric Generators (Diesel) and (Gas) 
are distributed in 66kV and 23kV, very close to the Load. One PVP and BES [23] System will 
be connected directly to the bus of SE-FIL 66kV. 


POWER TRANSFORMERS 


In general terms, the Power Transformers of system in study did not suffer great 
variations. In the SE-FIL, the PT that will serve as the source for PVP and BESS [23] [24] will 
be two (2) of 66 / 23kV of 30MVA each. 


COMPENSATION 


A flywheel is added to the SE-FIL and SE-MES to accelerate or brake under- 
frequency or over-frequency respectively and of an electrodynamic nature. 
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LOADS 


In general terms, the Loads of system in study did not suffer great variations [22]. The 
loads and DG were divided into parts equal or proportional to the loads of each SE in order 
to facilitate the practicality, performance and academic analysis of the systems under study. 


COMMERCIAL FEATURES 


In this part it was not investigated in depth contractual, but obviously the regulation 
of buying and / or selling the energy of the NEN, Diesel, Gas or Solar must be studied. This 
research focuses directly on the Voltage Stability of the SEP under study. 


CURRENT CONTROL ACTIONS FOR THE OPERATION OF THE ELECTRICITY 
NETWORK 

Currently the control actions are carried out by varying the voltages of the neighboring 
SE. Load shedding according to need or making a connection in parallel with a Generator 
Distributed (DG) and PVP/BES system during convenience operative or temporarily during 
times of high demand or bad weather. All this information would be processed thanks to the 
programming and intervention of the MATLABO FUZZY DESIGN [21] to carry out similar 
decisions. 

The actions of control of Voltage and Operation of the electromechanical equipment 
would be left to the decision of the FLPRC strategically located throughout the study PES. 


PHOTOVOLTAIC POWER PLANT AND BATTERY [24] ENERGY STORAGE 
SYSTEM (PREVISION) 

To the SEP studied, a PVP and BES [23] systems are added in the Electric Power 
Transmission Levels format, directly connected to the 66kV bus of the SE-FIL, as shown in 


figure 29. 


THE FOLLOWING DATA WERE EXTRACTED FROM SOLARGISO 


Dates in [15] 

Global tilted irradiation Yearly average: 5,418 kwh/m* per day. 

Air temperature Yearly average: 25ºC. 

Specific photovoltaic power output Yearly average: 4,284 kwh/m? per day. 
Performance ratio Yearly average: 79,1%. 

Address: Filadelfia, Boquerón, Paraguay. 

Geographical coordinates: -22,380451º,-060,03574º 

Time zone: UTC-04, America/Asuncion [PYT]. 
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Elevation: 139m. 

Land cover: Mosaic cropland (>50%) / natural vegetation. 

Population density: tinh./kmº?. 

Terrain azimuth: flat. 

Terrain slope: 0º. 

Below are pictures of project horizon and sun path. Also, they are detailed in the image 
of day length and solar zenith angle. 
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Figure 2. Project horizon and sun path. Adapted from [15]. 


In the same way, were obtained, the data of the Solar Energy Potential and 


Meteorology. 


Global horizontal irradiation (GHI): 5,139 kWh/mº per day. 
Direct normal irradiation (DNI): 4,707 kWh/m* per day. 
Diffuse horizontal irradiation (DIF): 2,010 kWh/mº per day. 
Ratio of diffuse to global irradiation (D2G): 0,391. 

Air temperature: 24,7ºC. 

Cooling degree days (CDD): 2.527 degree days. 

Heating degree days (HDD): 70 degree days. 
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Figure 3. Day length and solar zenith angle. Adapted from [15]. 


On the other hand are shown, but from the same information source [15], the 
Photovoltaic Electricity data and its performance. 

Specific photovoltaic power output (PVOUT specific): 4,284 kWh/kWp per day. 

Total photovoltaic power output (PVOUT total): 78.189.816,6 kWh 

Global tilted irradiation (GTI theoretical): 5,418 kWh/m2 per day. 

Long-Term PV Power Potential Data (25-Year Average) from the same information 
source [15]. 

Specific photovoltaic power output: 4,005 kWh/kWp per day. 

Total photovoltaic power output: 73.084.033,5 kWh. 

Performance ratio: 73,9 % 

Capacity factor: 16,7 % 


GLOSSARY 


*- | Cooling degree days: Quantifies energy demand needed to cool a building. 
“Cooling degree days” are a measure of how much (in degrees), and for how 
long (in days), outside air temperature was higher than a specific base daily 
average temperature (18ºC). Yearly and monthly values are aggregated from 
daily values 


- | Capacity factor: The ratio of an actual electrical energy output over a year to 
the maximum possible electrical energy output over a year expressed in %. 
The maximum possible power production is the AC installed capacity times the 
number of hours in a year, while the actual production is the amount of electricity 
delivered annually from the project. 


- Ratio of diffuse to global irradiation: Ratio of diffuse horizontal irradiation and 
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global horizontal irradiation (DIF/GHI). 


-  Diffuse horizontal irradiation: Average yearly, monthly or daily sum of diffuse 
horizontal irradiation. 


- Direct normal irradiation: Average yearly, monthly or daily sum of direct nor- 
mal irradiation. 


- Global horizontal irradiation: Average annual, monthly or daily sum of global 
horizontal irradiation. 


- Global tilted irradiation: Average annual, monthly or daily sum of global tilted 
irradiation. 


* | Global tilted irradiation (theoretical): Average annual, monthly or daily sum of 
global tilted irradiation without consideration of terrain shading. 


- Heating degree days: Quantifies energy demand needed to heat a building. 
“Heating degree days” are a measure of how much (in degrees), and for how 
long (in days), outside air temperature was lower than a specific base daily 
average temperature (18ºC). Yearly and monthly values are aggregated from 
daily values. 


- | Performance ratio: Ratio between specific AC electricity output of a PV system 
and global tilted irradiation received by the surface of a PV array (PVOUTspe- 
cific/GTI). 


- Specific photovoltaic power output: Yearly and monthly average values of 
photovoltaic electricity (AC) delivered by a PV system and normalized to 1 kWp 
of installed capacity. 


- | Total photovoltaic power output: Yearly and monthly average values of pho- 
tovoltaic electricity (AC) delivered by the total installed capacity of a PV system. 


*- Air temperature: Average yearly, monthly and daily air temperature at 2 m abo- 
ve ground. 


THEORETICAL DATA OF PHOTOVOLTAIC PANELS DISCLOSED BY THE 
SAME SOURCE OF INFORMATION [15] 

System size: Installed capacity: 5SOMWp 

PV module type: c-Si - crystalline silicon (mono or polycrystalline). 

Geometry of PV modules: Azimuth: 0º « Tilt: 21º. 

Relative row spacing: 2,5 

Inverter type: Centralized high-efficiency inverter [97.8% Euro efficiency]. 

Transformer type: High efficiency transformer [0.9% loss]. 

Snow and soiling losses at PV modules: Monthly soiling losses up to 3.5 % * Monthly 


snow losses up to 0.0 %. 
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Cabling losses: DC cabling 2 % * DC mismatch 0.3 % * AC cabling 0.5 Yo. 

System availability: 99.5 %.BES System. 

BESS [24] equivalent to 30% of the total installed capacity of the PVP, or 15MW. 
(tentative and theoretical values). [23] Local ESS is a local provider of storage capacity for 
storing and delivering energy. This equipment is very necessary when sunlight is not present. 
This storage equipment must be capable of delivering a certain amount of energy to the 
Electric Grid (EG). 

[23] The parameters of the Li-ion battery are defined via look-up tables based on 
experimental data. 

[24] The lifetime of each type of PV array is 20 years. Itis assumed that the replacement 
time of the BESS is four years. The minimum and maximum rated power values of each PV 
array are 200kW and 20MW, respectively. For the BESS, the minimum and maximum rated 
capacity values of each battery are 20kWh and 1MWh, and the minimum and maximum rated 
power values of each battery are 20kW and 1MW. 

Parameters of Batteries (Bullich-Massagué et al; 2020; Mariaud et al., 2017). 


Battery 1 Battery 3 
Description Lead-acid Nas 
Charging/discharging efficiency (%) 91,0 88,0 
Cycle number (cycle) 600,0 2500,0 
soc”! 0,4 0,2 
soc" 0,9 0,8 
Unit power capacity cost (S-kw"*) 400,0 500,0 
Unit energy capacity cost (S-kwh?) 600,0 1200,0 
O &M coefficient (%) 1,0 1,0 





Table 1. Parameter of Batteries [adapted from 24]. 


PHOTOVOLTAIC POWER PLANT AND BATTERY ENERGY STORAGE SYSTEM 
(IN OPERATION) 

Data of Operation [16]. 

Address: Del Salto Department, Uruguay. 

Geographical coordinates: -31.262800º, -57.863060º 

Name: PVP “EL NARANJAL”. 

Start-up by “INGENER GROUP” and managed by “COLIDIM S.A”. 

Installed Capacity (MW): 5OMWp 

Solar paneis and respective trackers: 200.000 panels. 

Area destined for PV Generation: 145Ha. Or 1.449.000m?. 

Photovoltaic Power Plant and Battery Energy Storage System For a proposed case 


study. 
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[23] [24] For this case, a PVP/BESS Complex is of 5OMWp (AC) and 15MW (DC) are 
proposed, with double solar panels (200.000), with solar flow followers or trackers [17], 
connected directly to the 66kV Buss of SE-FIL, as shown in the figure for an area of 150Ha. 


up to 2.40 m 
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up to 2.38 mt 





























Figure 4. Tracker-Dueitto: all in one. Adapted from [17]. 


Appearance of the PVP and BES System to be used in the SE under study. 


50 MWPp(AC) 
15MW (DC) 


BES 
(30%PV) 







2 x30MVA 


FILG6KkV 






Figure 5. PVP and BES System (own creation). 


Historical maximum demand (MW) in the SEP under study hour by hour. The different 
voltages of the buses involved are also displayed: 220, 66 and 23kV respectively. Simulation 
results are compared in RP (Power flows) vs. Values registered on February 23, 2021, the 
day that the historical maximum demand of 3,777MW was registered. The Load registered 
correspond to SE-LPA 
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Table 2 and Figure 6. Historical Maximum demand (MW), 02/03/2021. (own creation). 


RESULTS AND DISCUSSIONS. GENERAL EQUATION OF TRANSMISSION 
LOSSES 


Loss = 100% — (era) 6 
t 


(6) Where L are the Loads and G, are the total Generation. 


LatLasla 
Loss = 100% — (rmfttztita =) 7 
nao 00% DG, +DG>+---DGn +NEN 


(7) Where L are the Loads, DG, are Distributed Generation and NEN is the National 
Electricity Network. 


Loss = 100% — (DE Ra CEVETDESS) 8 

(1) Where L are the Loads, DG, are Distributed Generation, NEN is the National 
Electricity Network, PVP is Photovoltaic Power Plant and BESS is Battery Energy Storage 
System. 


The losses are calculated below for each notable case presented. 


Chapter Load Generation Loss 


2 99,6 69,4 
3 29,3 99,7 
4 26,3 26,3 





Table 3. Calculation of Transmission Losses for each case. 
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RECORD OF HISTORICAL VALUES VS PROPOSALS OF PVP AND BESS 


The following shows the registered values of Voltage Values, after the Intervention of PVP [15] [16], GD 
[12], NEN and BESS [23] [24] for generation in the record of maximum demand, Values in MW, 


HISTORICAL VOLTAGE IN PU. 2021/02/23 - SE-LPA 
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Figure 7. Voltage in different Buses of SEP in studied. 


It is observed that the Voltage values in 66 as 23kV respectively remain constant, 
thus maintaining a fine control of them. TAP positions are gained and Transmission losses 


decrease drastically. 


HISTORICAL MAXIMUN DEMAND in MW. 2021/02/23 - SE-LPA 
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Figure 8. Generation NEN, DG, PV and BESS as a result of joint operation. 


The historical Demand versus the generations NEN, DG [12] , PV [16] , BES [23] 
in terms of Active Power (MW) are plotted below. The great intervention of the additional 
generations is observed, maintaining the system with a lot of reliability, performance, 
sturdiness and with a wide Voltage control. Losses for Transaction decrease dramastically. The 
partial or total feeding of the feeders becomes very easier. 

The green curve corresponds to the Active Power registered in 220kV buses of the 
SE-LPA. The blue curve corresponds to the (simulated) Active Power in 220kV buses of the 
SE-LPA after the intervention of the other Electric Generations. The black curve corresponds 
to the generation (simulated) of the PV in 66kV buses of the SE-FIL. The yellow curve 
corresponds to the Generation (simulated) of the BESS in 66kV buses also of the SE-FIL. 
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And the red curve corresponds to (simulated) the entire DG in each Electrical Substation in 
the vicinity of the load (LPA, FIL and MES). 


PROGRAMMING IN MATLABO [21] THE NEW PROPOSAL OF PV AND BESS 


Programming for the new proposal of PV Injection and BES [23] system related to 
other generation sources to control Voltage and supply feeder service using MATLABO Fuzzy 
Logic Design [21] Tool. 


PV and BES System 
Fuzzy IN (%vsh) 


50% 


30% 





Figure 9. Behavior during the day regarding the injection of Active Power from 


PVP/BESS COMPLEX TO SEP AND THE SAME FOR THE NIGHT 
(THEORETICAL FORM AND OWN CREATION) 

Behavior during the day regarding the injection of Active Power from PVP to SEP and 
the same for the night (Theoretical form). These signals would be an input value for the Fuzzy 
Brain Together with another Voltage Signal, since this Power Electronics equipment does not 
have physical inertia, so it is impossible for it to control transients of an electrodynamic nature 
and it must have other equipment capable of facing such situations. For protection in this case 


was considered a Normal Voltage of Operation between 0,950 and 1,030 pu. 
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Figure 10. Voltage input in 66kV Buses for PVP/BESS Complex. 


Figure 38 shows another of the FB input signals to associate with the triangular 
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decision output. 


PV and BES System 
Fuzzy IN (%vsh) 


PV and BES System 
Fuzzy OUT (% vs MW) 





Figure 11. Inputs, processing and decision making. 


Process by which the inputs are processed by the fuzzy brain based on established 


rules to associate them with the triangular decision output. 


PV and BES System 
Fuzzy OUT (% vs MW) 





0 6,25 12,5 18,75 25 31,25 375 43,75 50 


Figure 12. Output with triangular decision after going through Fuzzy processing 


The next step would be to carry out the programming including the PVP and the BESS 


[24], for later conclusions. 


FIRST PROGRAMMING PROPOSALS IN MATLABG 


First programming proposals in MATLAB, taking into account the operational limits of 
Voltage and availability in the PV / BESS generation as input, and the Generation relation in 
MW for the establishment of the rules. 

After carrying out some programming in MATLABO, satisfactory results were obtained 
regarding the Participation of the PV Generation and BESS from 00:00 to 24:00 on a hot 
summer day and in full sun. The PVP / BESS Complex operating range has also been related 
in terms of Voltage limits from 0,950 to 1,030 pu; that is, ffom 62,7kV to 68kV. 

Topology of the construction and relationship of input, processor and output data, 
through the CF to the different FLRPC installed throughout the PES. 


PRELIMINARY RESULTS AFTER PROGRAMMING IN MATLABO 


MATLAB?º graphics defines a surface by the Z-coordinates of points above a rectangular 
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grid in the x-y plane. The plot is formed by joining adjacent points with straight lines. Surface 
plots are useful for visualizing matrices that are too large to display in numerical form and for 
graphing functions of two variables [21]. 


The following graph (Surface Graphics) shows the axes in three dimensions: 





GR iroo Os 


Participation% woltagePU 


Figure 13. Surface Graphics. Results of the programming of the PVP / BESS complex. 


In figure 13, Axis “X”, which indicates the Active Power Injection ffom the PVP /BESS 
complex (5 to 5SOMW) and that its operation is subject to the value of the 66kV bus voltage 
where the Installation is interconnected, participation of Solar Generation and storage (00:00 
to 23:00 of the day). 

As indicated above, the “Y” axis corresponds to the participation of Solar Generation 
and storage in percentage terms (see figure 37), reported hour by hour (00:00 to 24:00 a 
day). This corresponds directly to one of the inputs to the FB and it is strongly related to the 
FLRPC3 on the 66kV bus. 

As for the “Z” axis, itis also a voltage information in pu. input to the FB directly linked 
to the FLRPC3 of the 66kV bus where the equipment is connected. To limit the operability of 
the PVP / BESS complex, values of 0,950 < V <1,030 pu were established. These voltage 
values are established between approximately 62,7 /68kV, since this equipment lacks physical 
inertia, so for electrodynamic events, it is convenient to disconnect it below or above these 
values. FL Programming is found in Appendix 2. 


CONCLUSION OF THIS CHARTER 


With the DG installing directly on the load bus, a significant reduction in transmission 
losses has been achieved. Voltage control has become more efficient and practical. The 
possibility of autonomy in the form of islands is clearly observed when the loads are not fed 
from the NEN. 

The PVP with the BESS fully integrated to the 66kV bus of the SE-FIL makes the 
performance of the operation optimal, since it increases the reliability of the very radial power 
system in the area. 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 5 


The autonomy of the loads becomes a practicality, given the ease of operation between 
their own networks when the sources of the National Electricity Network are interacting, 
Distributed Generation (Diesel and Gas with lower Carbon emissions), Photovoltaic 
Generation (Generation clean and renewable) and the storage of Electric Energy through 
Battery Banks (BESS). 

For futures works, it would be interesting to investigate the dynamic behavior of all 
these sources interacting all at the same time, always looking for the operational and economic 
convenience of the Company. 


THE SATISFACTORY RESULTS, IN SUMMARY, ARE AS FOLLOWS: 


The programming in the Fuzzy Logic interface is very flexible in terms of the 
information to be processed, since different units of measurement can be handled and they 
are finally independent of the result to be obtained. 

The use of the AIFL is very powerful in terms of the complex decisions that the FB 
must make when the information conditions are minimal. 

A clear performance of the Al method used was demonstrated thanks to its 
practicality and flexibility, as long as the input and desired output information are correctly 
processed. 

Finally, the use of this AIFL demonstrates a strong versatility for any type of 


programming and complex decisions. 
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RESUMO: Este trabalho técnico apresenta 
um resumo dos estudos de energização no 
“Regime Transitório e Dinâmico” considerando 
o comissionamento da Linha de Transmissão 
de 500kV entre as Subestações AYOLAS e 
VILLA HAYES com 350km de extensão (sem 
ter reatores no final do extremo de VHA). Com 
o comissionamento desta linha de transmissão 
de 500kV, pretende-se interligar blocos de carga 
previamente energizados em níveis de 220kV, 
alimentados a partir da Usina Hidrelétrica ITAIPU 
(compartilhada entre o Paraguai e a Brasil). 
PALAVRAS-CHAVE: Estudos de transitórios 
eletromagnéticos  (ETEmg), Estudos de 
transitórios eletromecânicos (ETEmc), efeito 
Ferranti, Sobretensões temporárias (TOV), LT 
500kV AYO-VHA, Central hidrelétrica Itaipu (CH- 
IPU), régime permanente (RP). 


“STUDIES OF ELECTROMAGNETIC AND 
ELECTROMECHANICAL TRANSIENTS” 
OF THE TL 500KV AYOLAS-VILLA HAYES; 
ENERGIZATION WITHOUT REACTOR 
SINCE THE ITAIPU HYDROELECTRIC 
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POWER PLANT 


ABSTRACT: This technical work presents a 
summary of the energization studies in the 
“Transitional and Dynamic Regime” considering 
the commissioning of the 500kV Transmission 
Line between AYOLAS and VILLA HAYES 
Substations with 350km in length (without having 
reactors at the end of the VHA end). With the 
commissioning of this 500kV transmission 
line, the intention is to interconnect load blocks 
previously energized at 220kV levels, fed from 
the ITAIPU Hydroelectric Power Plant (shared 
between Paraguay and Brazil). 

KEYWORDS: Electromagnetic transient studies 
(ETSmg), Electromechanical transient studies 
(ETSmoc), Ferranti effect, Temporary overvoltages 
(TOV), LT 500kV AYO-VHA, Itaipu hydroelectric 
plant (CH-IPU), permanent regime (RP). 


11 INTRODUÇÃO 


Quando uma linha de transmissão (LT) é 
projetada, o comissionamento desta, em termos 
de energização, operação, proteção e aspectos 
energéticos ou económicos, devem justificar 
o alcance da mesma. Os estudos (inerentes 
a energização e operação) de RP, transitórios 
eletromagnéticos e eletromecânicos, devem 
prever o comportamento e o desempenho da 
operatividade da LT. Em alguns casos devido 
a assuntos políticos, económicos, carência de 
tempo ou situações fora do contexto técnico, não 
permitem completar todas as obras associadas 
ao comissionamento final da obra principal. A 
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LT 500kV AYO-VHA e seus elementos associados surgiu como um nexo para a integração 
dos sistemas elétricos de Argentina, Brasil e Paraguai. Estudos de interligação regional 
se encontram em desenvolvimento, além disso, a LT pode ser operada em sua condição 
transitória, para obter o maior proveito da mesma no sistema de potência de Paraguai. Este 
estudo aporta informação importante e relevante para a operação da LT dentro do sistema 
elétrico da ANDE. 

Finalmente, neste trabalho técnico se apresentam os resultados dos estudos 
da energização e operação (desde a CH-IPU), da LT de 500kV de AYO-VHA, com um 
comprimento aproximado de 350km, injetando em seu extremo terminal uma potência reativa 


de aproximadamente 350MVAr, devido principalmente ao efeito Ferranti. 


21 ABORDAGEM TEÓRICA. EFEITO FERRANTI E TRANSITÓRIOS 
ELETROMAGNETICOS 


2.1 LT DE GRANDES COMPRIMENTOS E EFEITO FERRANTI 


Em situações donde uma LT está operando em condição radial ou em períodos de 
carga leve, a tensão no terminal receptor poderia ser muito mais elevada que do emissor. 
Este fenómeno é denominado efeito Ferranti e converte a uma LT em uma “geradora de 


potência reativa” não lineal em função da tensão [1]. 





A equação (1) descreve que a potência reativa capacitiva (MVAr) é diretamente 
proporcional ao quadrado da tensão (Volts) e inversamente proporcional à reatância 
capacitiva (ohm) da LT. 

O efeito “Ferranti” descreve o estranho fenômeno que baixo certas condições de 
frequência e comprimento de linha um aumento na tensão pode ser observado na mesma, 
quando o extremo se encontra aberto, em relação a uma tensão de entrada senoidal. O efeito 
foi descoberto no final do século XIX durante a instalação de um sistema de distribuição 
baseado em CA na Grã-Bretanha. Na verdade, foi durante a guerra técnica e comercial que 
atingiu os Estados Unidos entre a Westinghouse e Nikola Tesla como principais propulsores 
e Thomas Edison. Os primeiros propuseram um sistema CA em quanto Edison propunha 
comercialmente um sistema de CD. No Reino Unido, foi “Sebastian Ziani de Ferranti” quem, 
como fervoroso defensor dos sistemas CA, instalou um sistema de distribuição de corrente 
alternada com níveis de média tensão e reduzindo os transformadores. Em um sistema de 
transmissão CA, Ferranti observou que, adicionando seções de distribuição adicionais, isto 
é, aumentando o comprimento total da linha de transmissão, a tensão na linha aumentava 
localmente (Figura 1). De fato, o fenômeno foi observado pela primeira vez na Linha de 
transmissão Deptford - London, que em algumas lâmpadas de fibra de carbono aumentou 
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a luminosidade quando seções foram adicionadas na distribuição. Deve notar-se que o 
sistema tinha uma carga de apenas um par de lâmpadas de baixo consumo, enquanto a 
potência do gerador excedia ligeiramente 935 kW. Assim, Ferranti determinou que o efeito 
foi acentuado com o LT aberto na sua extremidade [2]. 





Figura 1. LT entre os nós 'A'e 'B', separadas por uma distancia elétrica 'd”. 


2.2 TRANSITÓRIOS ELETROMAGNÉTICOS 


A análise de transitórios (o simplesmente transitórios) dos circuitos elétricos é tão 
importante como os analise de estado estável. Quando ocorrem estes transitórios, as 
correntes e tensões em algumas partes do circuito podem superar varias vezes os valores 
de comportamento normal e em consequência pode deteriorar o equipamento. O estado 
transitório pode ser distinguido de seu estado estável, mesmo que durante os transitórios, 
magnitudes como correntes, tensões, potência e energia, se modificam no tempo, enquanto 
que em estado estacionário permanecem invariáveis, é dizer, constante (em operação 
de CD) ou periódica (em operação de CA) que possui ângulos de fase e amplitudes 
constantes. A causa dos transitórios pode ser qualquer tipo de mudança nos parâmetros 
do circuito e/ou na configuração, que geralmente ocorre devido à comutação, curto circuito 
e/ou circuito aberto, mudanças de operação das fontes, etc. As mudanças das correntes 
e tensões durante os transitórios não são instantâneas e têm um tempo de retardo, além 
de isso que são extremadamente rápidos com uma duração de milissegundos ou mesmo 
microssegundos. Porém, estes mudanças rápidas nem são instantâneas (ou abruptos) porque 
os processos transitórios são logrados a través do intercambio de energia, que geralmente 
são armazenadas nos campo magnéticos das indutâncias e/ou nos campos elétricos das 
capacitâncias. Qualquer variação na energia não pode ser abrupta (porque a potencia é uma 
derivada da energia, p = dw /dt), e que está em contraste com a realidade física. Todas as 
variações transitórias, que também são chamadas respostas transitórios (o simplesmente 
respostas), se desvanecem e, logo de sua desaparição, se estabelece uma nova operação 
em estado estável. Então é possível dizer que o transitório descreve o comportamento do 
circuito entro dois estados estáveis: um que anterior às mudanças, e outro novo que surge 
depois das mudanças [3]. 
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31 ETEMG DA LT 500KV AYO-VHA DESDE LA CH-IPU 


Nestes tipos de estudos são também abordados os estudos de reconexão e rejeição 
de carga e devem se desarrolhar dentro das atividades do desenho básico dos seguintes 
objetivos principais [4]: 


* Avaliar as sobretensões transitórias máximas que vai impor nas subestações e 
os terminais das linhas de transmissão. 


* Avaliar a dissipação de energia dos descarregadores de linha, considerando 
o dimensionamento dessa equipe desde o ponto de vista da capacidade de 
absorção de energia (kJ/kV). 


* Verificar a adequação da coordenação de isolamento das estruturas de trans- 
missão de surtos de manobras, com a integração do projeto básico. 


3.1 CONSIDERAÇÕES PRINCIPAIS 


As manobras de energização da LT 500kV AYO-VHA, foram analisadas nas barras 
das subestações AYO e VHA para a verificação de possíveis sobretensões temporárias e/ 
ou permanentes das equipes involucrados da LT, sempre com fonte da CH-IPU, como se 
observa na figura 2. É abordado que os fenómenos analisados cumpram com os requisitos 
básicos em termos de energização transitória de linhas de transmissão elétrica. Os casos 
de energização da LT foram feitas a traves de um programa computacional de transitórios 
eletromagnéticos (EMTP - Electromagnetic Transient Program). 


CH-IPU 







LT 500kV AYO-VHA 350 km 


: TOPOLOGIA 







80MVAr ES-AYO . 
2x 375MVA 


VuSas 


' 
' 
[ Í f “a 


Bloco de Carga A Bloco de CargaB 


Figura 2. Circuito equivalente utilizado para os ETEmMg da LT 500kV AYO-VHA. 









Na figura 2 se observa o circuito equivalente a ser utilizados nos ETEmg. Serão 
mostrados os resultados obtidos da energização da LT desde cada extremo terminal a 
través dos interruptores 'S e 'S,. 
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3.2 CRITÉRIOS DE TENSÃO ADOTADOS 


* | Em barras de 500kV, para sobretensões transitórios de manobra (SIL), se ad- 
mite um nível de isolamento de 1.175kV com um margem do 15%, é dizer, os 
valores máximos de surtos de tensão obtidos nas simulações de manobras não 
devem superar os 999kV. 


* | Em barras de 220kV, se admite um SIL de 1.050 kV com um margem do 15%, 
é dizer, os valores máximos de surtos de tensão não devem superar os 892kV. 


* As tensões máximas de RP devem ser inferiores a 1,100 p.u. em barras próxi- 
mas à geração. 


* | Nos casos de estudo de energização da LT foram abordadas energizadas sis- 
temáticas, estadísticas e finalmente determinísticas para obter as maiores soli- 
citações de tensão nas barras involucradas. 


* | Desempenho dos descarregadores de 500kV de E < 3276kdJ. 


3.3 SEQUENCIA DE ENERGIZAÇÃO DOS ETEmg DOS DIFERENTES 
CENÁRIOS 


* Às energizações foram realizadas de forma sequencial e separadas. 

* | Os casos de estado estacionário são análogos a os de regimen transitório. 
* | Osistema é alimentado desde a CH-IPU (Figura 2). 

*— ALT 500kV AYO-VHA é energizada desde o interruptor 'S 

*— ALT 500kV AYO-VHA é energizada desde o interruptor 'S,. 


* | Os blocos de carga 'A'e 'B' são interligados a través dos interruptores S. e S,. 


3.4 CONDIÇÕES INICIAIS DOS CASOS DE FLUXO DE POTENCIA 


*— ACH-IPU alimenta ao HDVC de FURNAS e ao Sistema Interligado Nacional de 
Paraguai antes da energização da LT 500kV AYO-VHA. 


*  Atensão de pre-energização em barras da ES-VHA está na ordem dos 0,900 
p.u. (450kV). 


* | Os blocos de carga “A e B” estão previamente energizados com fonte da CH-I- 
PU. 


* | O extremo terminal da LT de 500kV a ser energizada conta como reator de 
neutro de 80MVAr e está conectada à barra de 500kV dessa subestação. 


* || Todos os elementos de controle de tensão estão operando em condições nor- 
mais e com o sistema de potência preparado para a energização com fonte da 
CH-IPU. 
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3.5 RESULTADOS OBTIDOS DOS ETEmg. 


A continuação, são apresentados os resultados obtidos dos ETEmg das energizações 
da LT desde o terminal de AYO e desde o terminal de VHA. OS resultados são apresentados 
em termos de tensão nas barras de 500kV das EESS VHA e AYO, e a energia dissipada pelos 


descarregadores de sobretensão das barras análogas a VHA e AYO. 


3.5.1 ENERGIZAÇÃO DA LT 500KV AYO-VHA DESDE A ES-VHA 
(interruptor “S.*) 


ES-VHA - 500kV (kV / pu) - Lado da CH-IPU 


Antes de 
energizar a LT 





Maior surto de |Post-energización 


dalT 





Energização desde VHA 








tensão 
kv pu | kv pu | kv pu 






ES-AYO - 500kV (kV / pu) - Terminal da LT o inicia soMvar ESAVO 1, 
feio 600 MVA 2x375 MVA qo 
Antes de Maior surto de |Post-energización T 


energizar a LT tensão da LT 220kV 


kv pu | kv pu | ky pu Becodcanaa [——— ————] Bocodecamas, 


0 0,000 892 1,785 538 1,076 








Figura 3. Energização da LT desde VHA e aberto no terminal de AYO. 


A figura 3 mostra as tensões das barras de VHA e AYO, quando a LT é energizada 
desde o interruptor “S,”. Não são observadas sobretensões temporais que fiquem fora dos 
critérios normais adotados. Não se registram tensões de post-energização que superem os 


limites estabelecidos dos 1,100 p.u. em RP. 


ES-VHA 500kV - Lado da CHHPU ESAYO 500kV - Terminal daLT 





Figura 4. Tensões em barras de VHA e AYO 500kV. 
Da “figura 4” se extraem as seguintes informações: 


* Máxima sobretensão transitória na “ES-AYO”, 892kV (f-f) = 728kV (f-n) 


<< 999kV (f-n) 


RMS 


PEAK PEAK' 


* | Tensão de estabelecimento a 1 s na “ES-AYO”, 310kV (f-n) ,s << 318kV (f-n) 
é dizer, 1,076 p.u. < 1,100 p.u. 


RMS” 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 6 


Descarregadores de Sobretensão — 500kV 
ES-VHA - Lado da CHAPU 


Descarregadores de Sobretensão — 500kV 


ESAYO - Termina da LT 


0.2 0.4 0.6 0.8 (9 1 





Figura 5. Energia dissipada pelos descarregadores de sobretensão de VHA e AYO 500kV. 


Com respeito ao rendimento dos descarregadores de sobretensão ZNO de 500kV, a 


“figura 5” mostra a seguinte observação: 


* Não se observa excitação dos descarregadores de sobretensão ZNO de 500kV 
da ES-VHA (barra mais próxima à CH-IPU. 


- É registrada a excitação dos descarregadores de ZNO da ES-AYO, com uma 
dissipação de energia na ordem dos 16kJ, muito inferior aos 3.276kJ estabele- 
cidos como critério. 


3.5.2 ENERGIZAÇÃO DA LT 500KV AYO-VHA DESDE A ES-AYO 
(interruptor “S,”) 


SE-VHA - 500kV (kV / pu) - Terminal da LT 








Antes de Maior surto de |Post-energización 
energizar a LT tensão daLT ME 
CH-IPU Energização desde AYO 
kv pu | kv pu | kV pu 


350 km 


0 0,000 EE. 1,946 733 1,467 





SE-AYO - 500kV (kV / pu) - Lado da CH-IPU Asa ESVHA somvar ESAYO . 
6: coomva 2x375MVA H% 
Post-energización top? *ep 


da LT 


Maior surto de 
tensão 


Antes de 
energizar a LT 





















kv pu | kv pu | kv pu 





Figura 6. Energização da LT desde VHA e aberto no terminal de VHA. 


A figura 6 mostra as tensões das barras de VHA e AYO, quando a LT é 
energizada desde o interruptor “S,”. Não são observadas tensões temporárias 
que ultrapassem os valores abordados nos critérios em quanto a surtos de tensão 
transitórias, porém, as tensões de post-energização registram valores muito por 
encima das tensões permitidas de regimen permanente como V < 1,100 p.u. 
(efeito Ferranti). 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 6 


ES-VHA 500kV - Terminal da LT ES-AYO 500kV — Lado da CHHPU 





485kV 
0. 01 


Figura 7. Tensões em barras de VHA e AYO 500kV. 


Da figura “figura 7” se extraem as seguintes informações: 


* Máxima sobretensão transitória na “ES-VHA”, 973kV (f-f) 
< 999kV (f-n) praxe 


= 794kV (f-n) 


RMS 


PEAK 


* | Atensão de estabelecimento a 1 s na “ES-VHA”, 423kV (f-n) ,us>> 318kV (f-n) 
nus: é dizer, 1,467 p.u. >> 1,100 p.u., 3/% mais do estabelecido nos critérios de 
tensão de RP. 


Descarregadores de 
Sobretensão — 500kV 
ESVHA- Terminal da LT 





Figura 8. Energia dissipada pelos descarregadores de sobretensão de VHA e AYO 500kV. 


Com respeito ao rendimento dos descarregadores de sobretensão ZNO de 500kV, a 


“figura 5” mostra a seguinte observação: 


* | São observados uns valores dissipação de energia dos descarregadores da 
ES-VHA de comportamento crescente que poderia ocasionar a perda de isola- 
mento dos elementos do sistema, principalmente para os casos de tensão de 
post-energização (RP). Este fenómeno se intensificaria com a perda e esgota- 
mento dos controles de tensão do sistema de potência como Centro Estático de 
Reativos, Banco de capacitores, TAP e outros. 


- É registrada a excitação dos descarregadores de ZNO da ES-AYO, com uma 
dissipação de energia na ordem dos 20kJ, muito inferior aos 3.276kJ estabele- 
cidos como critério. 
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41 ABORDAGEM TEÓRICA. TRANSITÓRIOS ELETROMECÂNICOS 


Nesta parte do trabalho técnico se abordará mis fenómenos eletromecânicos 
associados à energização da LT de 500kV AYO-VHA, e la explicação dos resultados obtidos 
em termos dinâmicos do sistema de potência. O objetivo destes estudos é o de complementar 
e respaldar os resultados obtidos previamente nos análises de regimen permanente e de 
transitórios eletromagnéticos. Por enquanto, se pretende fazer um barrido de todas as 
variantes possíveis dos fenómenos que poderiam acarrear como o resultado da energização 
da LT de 500kV. Os estudos de estabilidade eletromecânica dos sistemas elétricos estão 
relacionados com o análise do comportamento de “este” despois das perturbações. O tipo de 
perturbação do modelado que se deve utilizar na representação do sistema elétrico. Como 
resultado destas perturbações, que geralmente se devem a mudanças estruturais repentinos 
na rede elétrica, o sistema deixa o ponto de funcionamento estável nele que se encontrava 
e tende a estabelecer em outro ponto de operação. As unidades geradoras estão sujeitas a 
acelerações e desacelerações de tal intensidade que certas unidades ou grupos de unidades 
podem perder sincronismo entre elas o com o sistema. Dependendo da natureza e duração da 
perturbação, o comportamento dos controles eletromecânicos das unidades geradoras pode 
ser amortecido ou não, terminando em um novo ponto de funcionamento estável o colapso 
do sistema [6]. Possivelmente, a classificação mais importante do fenômeno dinâmico é o 
rango do tempo de resposta natural. A classificação típica é mostrada na figura 9. O conceito 
similar foi apresentado em [7]. 


DE escargas 
aipia 
Surtos | Surtosdetensão — | Surtosdetensão — 
Transitórios no 
E stator e ressonância 
sub-síncrona 





E Pa 
ramo 


C a de Rr a de 
Rr de velocidade e carga 


Caldeiras e dinâmicas de 
longo prazo 


107 10º 10º 0,1 10 10º 10º 
(segundos) 


Figura 9. Rango de tempos do fenómeno Dinâmico. Adaptado de [7]. 
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5 | ETEMC DA LT 500KV AYO-VHA DESDE LA CH-IPU 


Na seguinte seção se analisaram as duas condições topológicas quando a LT 500kV 
AYO-VHA é energizada desde a CH-IPU como o objetivo de determinar o comportamento 
dinâmico dos sistema de potência. Foram feitos todos os modelos dinâmicos do sistema de 
potencia estudado como central elétrica de geração, compensadores estáticos de potência 
reativa e outros parâmetros elétricos (transformadores, LT, cargas, capacitores, etc.). 


5.1 CONSIDERAÇÕES PRINCIPAIS 


* | Os blocos de carga estão previamente energizadas em níveis de 220kV com 
fonte da CH-IPU. 


* | A energização da LT é feita a partir dos interruptores “S ” e “S,” e de forma 
separada. 


* | ALT é energizada desde o interruptor “S ”. Logo os blocos de carga “A” e “B” 
são interligados a través do interruptor “S,” (figura 3). 


* ALT é energizada a través do interruptor “S,”. Logo os blocos de carga “A” e 
“B” são interligados a través do interruptor “S,” (figura 6). 


* | Todos os casos dinâmicos são análogos aos casos de RP e de transitórios ele- 
tromagnéticos. 


* Através do uso de um programa computacional de ETEmc, foram feitas as aná- 
lises de sobretensões temporais a frequência industrial fundamental (TOV). As- 
Sim, OS principais resultados obtidos a frequência industrial são apresentados a 
seguir e respaldados por os estudos de RP e de ETEmMO. 


5.2 RESULTADOS OBTIDOS DOS ETEmc. 


5.2.1 Energização da LT de 500kV AYO-VHA desde o terminal de VHA 


As seguintes análises correspondem à energização a través do interruptor 'S * com 
fonte da CH-IPU. 





Sobretensões temporárias e/ou permanentes Sobretensões temporárias e/ou permanentes 
de energização da LT. ES-VHA — Lado da CH-IPU. de energização da LT. ES-AYO - Terminal da LT. 
a) Tensão dinâmica | ão dinâmi 
e, 4 ita Tensão dinâmica 


ia 


| I 
| Fechamento do anel de 500kV 
entre a ES-VHA e ES-AYO 







Fechamento do anel de 500kV 
entre a ES-VHA e ES-AYO 





Tensão de pre-energização da LT | 


| 
es + 
el ACT. 


Tensão de pre-energização da LT 

















tempo -segundos Tempo - vegundos 


Figura 10. TOV em barras da ES-VHA e ES-AYO com posterior fechamento de anel. 
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Para as condições de operação estudadas, energizações da LT 500kV AYO-VHA 
desde o terminal de AYO com fonte de IPU, se observa tensões dinâmicas e de post- 
energização dentro dos critérios estabelecidos como são mostradas na "figura 10". Porém, 
é importante ter em conta que o valor da tensão de RP na barra da ES-AYO (extremo terminal 
aberto da LT) está próximo aos 1,100 p.u. (Limite de tensão de régime permanente em 


emergência). 




















Sobretensões temporais e/ou permanentes Sobretensões temporais e/ou permanentes 
de Energização da LT. Barra de 500kV de la CH-IPU e ES-VHA de energização da LT. Barras de 220kV 
né cá AYO 220kV | = 
k TRI220kV 
| Energização da LT =! W. VHA 220kV | 
il 500kV AVYO-VHA | Fechamento do anel de 500kV | | uu! RE Pç E 
| , u ' ' entre a ES-VHA e ES-AYO = Fa 


e) 7 SAN 220kV | e GUA 220kV | 


MD 220kV 








Tensão de pre- 
energização da LT | 











1111111 eee 





Tempo - segundos 


Figura 11. Á esquerda TOV nas barras de 500kV da CH-IPU. À direita TOV em barras de 220kV de 
varias subestações. Energização da LT desde o terminal de VHA. 


A “figura 11”, à esquerda, mostra a tensão em barras de 500kV da SE-MD. Se 
observam tensões na ordem dos 1,42 p.u. (521kV) de corta duração, não representando 
valores que estimulem anomalias ao controle dos geradores da CH-IPU. Na mesma figura 
a direita, são mostradas as tensões temporais e/ou permanentes em barras de 220kv de 
varias estações do sistema analisado. Observe-se que são registradas tensões sustentadas 
por debaixo dos 1,100 p.u. (242kV), inferior aos valores considerados nos limites. Aos 70 
segundos é fechado o anel de 500kV entre as barras de VHA e AYO, não observando tensões 
que violem os limites exigidos. Em consequência, a energização desde o terminal da ES-VHA 


faz que seja mais viável. 


5.2.2 Energização da LT de 500kV AYO-VHA desde o terminal de AYO 


As seguintes análises correspondem à energização a través do interruptor 'S,” com 
fonte da CH-IPU. 
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Sobretensões temporárias e/ou permanentes Sobretensões temporárias e/ou permanentes 
de energização da LT. ES-VHA - Terminal da LT. de energização da LT. ES-AY O — Lado da CH-IPU. 


A fis eee | q E | | 






Tensão dinâmica 


Tensão dinâmica | 


on! 1 Tensão de pre-energização da LT | Í 1,100 p.u. 











LL. eee eee e—-—s. eee 


Figura 12. TOV em barras da ES-VHA e ES-AYO. 


A “figura 12º” mostra que para a condição de energização analisada são observadas 
tensões por encima dos critérios de operação que foram vistos nos resultados de régime 
permanente e de ETEm9g, em consequência, (desde o ponto de viste deste análise) não é 
recomendável energizar a LT desde o terminal da ES-AYO ('S,). Também, como consequência 
de estes elevados níveis de tensão, existe altas possibilidades de que as equipes de controle 
de tensão fiquem fora de serviço empiorando a o fenómeno da sobretensão. 

















Sobretensões temporais e/ou permanentes Sobretensões temporais e/ou permanentes 
de Energização da LT. Barra de 500kV de la CH-IPU e ES-VHA de energização da LT. Barras de 220kV 
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Figura 13. Á esquerda TOV nas barras de 500kV da CH-IPU e ES-VHA. À direita TOV em barras de 
220kV de varias subestações. Energização da LT desde o terminal de AYO. 


A “figura 13”, à esquerda, mostra a Tensão em barras de 500kV da CH-IPU. São 
observadas tensões na ordem dos 1,020 p.u. (510kV), e não são observadas anomalias no 
controle de tensão dos geradores da CH-IPU. Na mesma figura à direita, são registradas 
tensões temporais e/ou nas barras de 220kV. Neste caso aparecem tensões sustentadas 
não permissíveis acima dos 2,300 p.u. (506kV), donde os controles de tensão se encontram 
esgotados e apeirando o isolamento dos componentes do sistema elétrico, em consequência, 
não é recomendável fazer a energização da LT de 500kV desde este terminal. 
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61 CONCLUSÕES FINAIS DA ANÁLISE REALIZADOS 


6.1 CONCLUSÕES ACERCA DOS RESULTADOS DE ETEmg. 


Finalmente, de acordo às resultados obtidos em simulações de análises de transitórios 
eletromagnéticos da energização da LT 500kV AYO-VHA sem reator, é concluído que quando 
é abordada a topologia estudada, não é recomendável energizar la LT desde o terminal 
“S, » devido às altas solicitações de tensão provocadas por a ausência de reatores nas 
extremidades da LT (efeito Ferranti) principalmente por altos valores de tensão de post- 
energização da LT. Por conseguinte, a LT de 500kV solo deve ser energizada e operada 
desde o interruptor “S, ”. 

Para realizar a energização do terminal de VHA é necessário manobrar previamente 
bancos de capacitores do sistema para mitigar os impactos produzidos pela energização da 
LT, principalmente em situações posteriores à energização. Al mesmo tempo, és necessário 
manter tensões na ordem dos 0,900 p.u (450kV) em barras de 500kV da ES-VHA no momento 
da energização da LT. Esta medida estaria destinada a mitigar os impactos resultantes da 
energização da LT. 


6.2 CONCLUSÕES ACERCA DOS RESULTADOS DE ETEmc. 


Com referência às análises obtidos em simulações de transitórios eletromecânicos 
da energização da LT 500kV AYO-VHA e posterior fechamento do anel entre a ES-VHA e 
ES-AYO é concluído que de acordo às condições de operação estudadas, são observadas 
tensões dinâmicas e de post-energização dentro dos critérios estabelecidos quando a LT é 
energizada desde o interruptor “S ” e posterior fechamento desde o interruptor “S,”. Além 
disso, é importante ter em conta que os valores de RP em barras de 500kV estão próximas 
aos 1,100 p.u. (550kV). 

Finalmente não é recomendável energizar a LT desde o terminal da ES-AYO devido à 
altas solicitudes de tensão principalmente as do post-energização (resultados já registrados 
nos ETEmog.), produzidas e intensificadas por a alta injeção de reativos na LT (efeito Ferranti) 
e com ausência de reatores nas extremidades da LT. Por conseguinte, a LT 500kV AYO- 
VHA solo pode ser energizada e operadas desde o interruptor “S,” da ES-VHA quando es 
alimentada com fonte da CH-IPU. 
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RESUMO: Transmissão de energia sem 
fio (Wireless Power Transfer, WPT) é uma 
opção de ganho em mobilidade e conveniência 
enquanto dispositivos elétricos são carregados. 
Metamateriais são usados para aumentar a 
eficiência de transmissão da energia por meio de 
acoplamento magnético ressonante. Um sistema 
WPT foi implementado em um software capaz 
de solucionar problemas eletromagnéticos em 
3D, onde três configurações foram simuladas: 
inicialmente sem Metamateriais, com Split 
Ring Resonators, e com uma linha espiral 
supercondutora (Metamateriais supercondutores) 
que foi desenvolvida neste trabalho. Uma 
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SUPERCONDUTORES 


investigação da potência e da eficiência destes 
sistemas foi realizada por meio de simulações. 
A distância entre as bobinas foi aumentada 
desde 4 cm até 10 cm, e com desalinhamento 
horizontal que variou até 3 cm. Os Metamateriais 
mostraram-se eficientes conforme se pode ver 
nos resultados. 

PALAVRAS-CHAVE: Transmissão de Energia 
Sem Fio; Acoplamento Magnético Ressonante 
Forte; Metamateriais; Supercondutividade. 


WIRELESS POWER TRANSFER 
THROUGH COUPLED MAGNETIC 
RESONANCE WITH CONVENTIONAL 
AND SUPERCONDUCTING 
METAMATERIALS 


ABSTRACT: Wireless Power Transfer (WPT) 
is an option to gain mobility and convenience 
while charging electrical devices. Metamaterials 
are used to increase the energy transmission 
efficiency by coupled magnetic resonance. AWPT 
system was implemented in a 3D electromagnetic 
solver, where three configurations were 
simulated: initially without Metamaterials, with 
Split Ring Resonators, and with a superconductor 
spiral line (superconducting Metamaterials) that 
was designed in this work. An investigation of 
the power and efficiency of these systems was 
carried out through simulations. The distance 
between the coils was increased from 4 until 10 
cm, and the horizontal misalignment varied up 
to 3 cm. The metamaterials showed themselves 
efficient as can be seen in the resulis. 
KEYWORDS: Wireless Power Transfer; Strongly 
Coupled Magnetic Resonance; Metamaterials; 
Superconductivity. 
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11 INTRODUCTION 


MMs are artificial structures effectively homogeneous, that is, their average cell size is 
much smaller than the wavelength. Because of it, the refractive phenomenon overlaps to the 
scattering and diffraction in the propagation of the wave [4]. 

MMs are also called Left-handed materials (LHM) because they do not obey the “right- 
hand rule”. Furthermore, these elements can present a negative refractive index (NRI), and 
for that, they can change the behavior of the wave according to the reversal of Snell's law. 
In order to obtain an NRI, a metamaterial must have at least the electric permittivity or the 
magnetic permeability negative in the same frequency range. 

The Wireless Power Transfer (WPT) is a technology able to transmit electromagnetic 
energy from a source to an electrical charge between a gap in the air. It can offer energy to 
small devices [1], biomedical implantis, portable devices, networks sensors, Internet of Things 
(loT), robots and electrical vehicles [2]. 

This technology is an option to replace the traditional methods of transmission by 
cables. It features advantages over security in low frequency, mobility, and convenience. 
Furthermore, it increases flexibility on devices whose battery replacement is expensive or 
dangerous [3]. 

There are many techniques for WPT, among them the strongly coupled magnetic 
resonance (SCMR) that was used in this work. This technique works over the magnetic field 
oscillation between two coupled coils that operate in resonance. 

However, the WPT has some problems regarding the cost of implementation and the 
power decay with distance. The efficiency decay comes from the increase of the distance 
between the source and load, misalignment between the transmitter and receiver, and 
inherent problems like radial propagation and reflection. 

In order to mitigate losses in WPT, and consequently increase the transmission 
distance and efficiency, the application of metamaterials (MMs) has been studied. MMs are 
made by natural elements and can present unnatural properties. When properly arranged and 
submitted to electromagnetic fields, they work like perfect absorbers. 

There are two types of MMs. The “conventional metamaterials”, such as the Split Ring 
Resonators (SRRs) and the thin wires (TWs), and currently “superconducting metamaterials” 
have been developed. In these two forms, the MMs have been employed to increase the 
electromagnetic waves absorptivity capturing the magnetic flux dispersed in the air. 

In this work, an investigation of the applicability of metamaterials is made. The behavior 
of conventional and superconducting MMs in wireless power transfer by strongly coupled 
magnetic resonance SCMR is evaluated when there are load variation and misalignment. A 
spiral line of superconducting MM was projected. The system is simulated without MM and 
with conventional MM. An investigation of the power and efficiency behavior according to 
the increase in the distance between the transmitter and receiver coils and to the horizontal 
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misalignment between them, for all system configurations, is also made. 

Section Il presents the theory of conventional and the called superconducting MMs. 
Section Ill presents the simulated systems and the design of the superconducting MM unit 
cell. Section IV presents the results and the Conclusion is in Section V. 


21 METAMATERIALS 


MMs are artificial structures effectively homogeneous, that is, their average cell size is 
much smaller than the wavelength. Because of it, the refractive phenomenon overlaps to the 
scattering and diffraction in the propagation of the wave [4]. 

MMs are also called Left-handed materials (LHM) because they do not obey the “right- 
hand rule”. Furthermore, these elements can present a negative refractive index (NRI), and 
for that, they can change the behavior of the wave according to the reversal of Snell's law. 
In order to obtain an NRlI, a metamaterial must have at least the electric permittivity or the 
magnetic permeability negative in the same frequency range. 


A. Conventional metamaterials 


Among the so-called conventional metamaterials, there are two very used models: Thin 
Wires and Split Ring Resonators. The TWs models have a negative electric permittivity, and 
the SRRs have a negative magnetic permeability. In this work, it is used SRR in simulations 
because of the SCMR technique, where the magnetic part of the evanescent fields prevails. 
The ring resonators have a cut to operate at a resonance frequency in which the 
wavelength is significantly greater than its diameter. 
B. Superconducting metamaterials 


The superconducting MMs are different in their constructive characteristics comparing 
to the conventional ones. To exhibit superconducting behavior, the material must be below the 
critical temperature (), when the MM can present the Meissner effect and have its resistivity 
considerably reduced and no magnetic flow inside the conductor. 

Low losses are achieved in superconducting MMs, higher quality factor, and smaller 
wavelength when compared with the conventional MMs [5]. Furthermore, superconducting 
MMs can be miniaturized, and then they can be applied in small and medium-sized devices. 

In this paper, the two types of metamaterials and their respective behaviors in the 
simulated system of WPT are investigated. 


31 SIMULATED SYSTEMS 


This section are presented the general data for the three system configurations 
simulated. In this work, a superconducting metamaterial unit cell was also projected. 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 7 HE 


A. General data 


The system was simulated on CSTÊ Studio software using the Finite Integration 
Technique (FIT) in time domain. The transmitter and receiver copper coils are flat and have 
7 loops, width of 2.49 mm, spacing of 1.0 mm, inner radius of 22 mm, outer radius of 47 mm. 
Furthermore, this system has an inductance of 4.73 uH, relative electrical permittivity () of 4.8, 
tangent loss (6) of 0.015, as shown in Fig. 1. These data and geometry were obtained in [6]. 
These coils are in an Flame Retardant 4 (FR4) substrate with dimensions of 100 mm x 100 
mm x 1.5 mm. 

In the transmitter coil was used a discret port of voltage, defined by a starting point 
and an endpoint. These two points were connected by a perfectly conducting wire (visualized 
by a thick blue line) and the respective port source (indicated by a red cone) in the center 
of this wire, as shown in Figures 1 and 2. This port type realizes an ideal voltage source, 
exciting with constant voltage amplitude. This discret port is used in all simulations, c.a., with 
a peak voltage of . In the receiver coil, a resistor of 100 O is used (in the without misalignment 
simulations between the coils). 





Fig. 1. (a) Transmitter coil (b) Receiver coil. 


The operating frequency in which the system was projected is 1 MHz. In this frequency, 
there are many applications in WPT such as self-resonant structures [7] and biomedical 
implants [8]. 

The percentage efficiency (n %) can be calculated through the relation between the 
load power (P,) and source power (P), both in watts: 


P, 


Ss 


B. SCMR without any Metamaterial 


Initially, the system was simulated with the transmitter and receiver coils with only air 
between them, as shown in Fig. 2. 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 7 





Fig. 2. SCMR system without metamaterial: (a) Complete system. Lateral view showing details of 
transmitter coil (b) and receiver coil (c). 


C. SCMR with conventional metamaterial 


The SRR model simulated in this work is based on the prototype of [4]. The metamaterial 
slab has 12 unit cells arranged in 3 columns and 4 lines, 300 cm? in total and it is placed at 
haliway between the coils, as shown in Fig. 3. 





Fig. 3. SCMR system without metamaterial: (a) Complete system. Lateral view showing details of 
transmitter coil (b) and receiver coil (Cc). 


Every unit cell of the slab has two copper rings presenting width of 1.63 mm, spacing 
of 1.47 mm, height of 0.02 mm, inner radius of 19.87 mm, with overall of 50 mm x 50 mm x 1.5 
mm and one capacitor of 47 nF in each ring. 


D. Design of the Superconducting Metamaterial Unit Cell 


A superconductor unit cell was designed in this work. Its dimensions were reduced 
when compared with conventional MM. The substrate of has dielectric constant of 25, tangent 
losses (0) of 0.0005, and critical current density of in 77 K. Furthermore, its dimensions were 
defined with an average size much smaller than the wavelength to respect the effective 
homogeneity principle [4]. The dimensions for every cell are 40 mm x 40 mm x 0.5 mm, and 
the area of the whole slab is 144 cm? (9 unit cells in 3 lines and 3 columns, 52% lower the 
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conventional slab simulated), as shown in Fig. 4. 

The superconductor metamaterial has a spiral line geometry of YBCO, a superconductor 
material that generally is used in a critical temperature of 90 K. After the definition of the unit 
cell, the dimensions of the spiral line were defined after computer simulations and analysis: 
an outer radius of 38 mm and an inner radius of 2 mm. 

The equivalent circuit of a superconducting MM is an RLC circuit and it is necessary to 
place capacitors to achieve the resonance. From [9] we have: 


C(pF) = 0.035D, (mm) + 0.06 (2) 
where C is the capacitance in pF, with a value of 1.39 pF in this case. 
After that, the geometric inductance (£) was defined from: 


obtaining L= 0.0182 H, where f is the resonance frequency. 
The number of loops ('N) of the spiral conductor is determined by: 





D 


avg Ho 
where the factor a is: 
(Do = D;) 


— NTO 7 5 
“(Do + Dj) = 


In (5) D, is the outer diameter, D is the inner diameter and the average diameter is 
calculated as: 


(Do — D;) 
Davg = 5 (6) 


It was obtained N = 43.37, which was considered N=43. The designed unit cell is shown in 
Fig. 4 (a) and the complete system in Fig. 4 (b). 
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Fig. 4. Unit Superconducting MM: (a) Unit cell (b) Slab of superconducting metamaterial (c) Complete 
System. 


41 RESULTS 


In this work, a transmission-line (TL) network periodically loaded using 

The results of the simulations and the comparison between them are presented in this 
Section. The system was simulated without MM, with Conventional and Superconducting 
MM. 


A. Comparison of the configurations without misalignment 


The distance between the coils varied from 4 cm to 10 cm. The system was simulated 
without MM, and with conventional and superconducting MM. The comparison can be seen 
in Fig. 5. 


System efficiency comparison 


————— Without MM 
Conventional MM 
oupercondutor MM 


Efficiency [%] 





Distance [cm] 


Fig. 5. Efficiency comparison. 


The best results were at 4 cm of distance between the transmitter and receiver coil. 
Without the MM, the power in the load resistor was 45.14 W, and the efficiency of 13.23%. 

With the conventional metamaterial, the power in the load was 69.59 W and efficiency 
of 23.46% in the transmission, as shown in Table |. There was an increase of 77.32% in the 
efficiency and 54.16% in the power after inserting the conventional MM. 
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In the presence of the superconducting metamaterial, the efficiency was 16.11% (an 
increase of 21.76%) and the power was 49.34 W (an increase of 9.30%) when compared to 
the simulations without metamaterial. 


The MMs were responsible to direct the disperse magnetic flux. Both types of MM were 


able to increase the power flux and transmission efficiency. 


one | without Conventional MM Superconduct. MM 


Distance (cm) | P(W) n (%) P(W) n (%) P(W) n (%) 


45.14 13.23 69.59 23.46 49.34 16.11 
2108 a00 | 1864 [| 2310 | 765 


ren 
NE 
[e [om [as [ais | neo [rom | om 
[e [ts [om [rm | ru [as | om 
[e om [os [em | ss [r06 | om 


Table 1 - Efficiency and power comparison. 





B. Comparison of the configurations with misalignment 


An investigation was carried out on the horizontal misalignment between the 
transmitter and receiver coils. It is important in cases where there is no absolute certainty 
about the position of the recipient (biomedical implants, for example). The coils were initially 
with their centers 4 cm apart and then there was a horizontal displacement from 1 to 3 cm, 
as shown in Fig. 6. The simulations, again, occurred without MM, with the conventional and 
superconducting MM and there was a comparison between them. 





Fig. 6. Displacement of the receiver unit. 
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The system efficiency shows an increase after the insertion of the metamaterial as can 
be seen in Fig. 7. 


Efficiency comparison with horizontal misalignment - O until 6 cm 


—— Without MM | 
Conventional MM | 
Superconducting MM I 


Efficiency [9%] 





Distance [cm] 


Fig. 7. Efficiency comparison with misalignment between coils. 


Even with 3 cm of displacement, the efficiency was increased with metamaterials: 
49.41%, from 7.67 to 11.46% with the conventional and 20.33%, from 7.67 to 9.23% with the 
superconductor, as shown in Table Il. 


Misalignment á ã a 
0 | 152 20.52 161 


1314 2045 1533 
108 1051 1260 


Table 2 - Efficiency and power comparison with misalignment. 





Although the system is very sensitive to the misalignment, it still presents a better 
efficiency compared to no MM setup. Thus, the MM was able to increase efficiency in wireless 


power transfer even with horizontal misalignment between coils. 


5 | CONCLUSION 


After inserting the two configurations of MMs and comparing their results with the 
simulations without MM, it was verified an increase of 77.32% and 21.76% in the efficiency 
with conventional and superconducting respectively. Even with a horizontal misalignment of 3 
cm between the transmitter and receiver coils the efficiency was also increased after inserting 
the MM. 
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RESUMO: A degradação induzida por potencial 
(PID) é uma problemática que nos últimos anos 
tem sido foco de pesquisas e estudos sobre 
o desempenho do módulo fotovoltaico (PV) 
em condições de campo, tendo em vista as 
consequências ocasionadas por esta degradação. 
Mesmo com um extenso material nessa área, a 
compreensão sobre o fenômeno PID ainda é 
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REVISÃO 


incompleta, mas deve-se levar em consideração 
que a diversidade tecnológica e ambiental são 
fatores que implicam nas técnicas de reversão. 
Este artigo visa realizar uma revisão crítica a fim 
de fornecer uma visão geral e ampla da literatura 
disponível para promover a compreensão do 
estado atual da pesquisa PID. O papel consiste 
em apresentar as definições do mecanismo PID 
embasadas por estudiosos e pesquisadores bem 
como a influência da temperatura, umidade e da 
tensão na progressão do PID e as metodologias 
de detecção e reversão e as medidas preventivas 
em módulos PV c-Si. 

PALAVRAS-CHAVE: Detecção, revisão literária, 
PID, fotovoltaico, reversão. 


POTENTIAL-INDUCED DEGRADATION 
(PID): REVIEW 


ABSTRACT: Potential induced degradation (PID) 
is a problem that in recent years has been the 
focus of research and studies on the performance 
of the photovoltaic module (PV) under field 
conditions, in view of the consequences caused 
by this degradation. Even with extensive material, 
the understanding of the PID phenomenon is still 
incomplete, but it must be taken into account 
that technological and environmental diversity 
are factors that imply in reversal techniques. 
This article aims to conduct a critical review in 
order to provide an overview and broad view of 
the literature available to promote understanding 
of the current state of PID research. The role is 
to present the definitions of the PID mechanism 
supported by scholars and researchers as well 
as the influence of temperature, humidity and 
tension on the progression of the PID and the 
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detection and reversal methodologies and preventive measures in PV c-Si modules. 
KEYWORDS: Detection, literary review, PID, photovoltaic, reversal. 


11 INTRODUCTION 


The relevance of photovoltaic generation (PV) worldwide, along with the growing 
energy demand, has led to greater participation of this source in the energy matrices of the 
countries. The wide use of photovoltaic systems has led to greater investment in research and 
innovation of this energy source, seeking to optimise cell efficiency, reduce production costs, 
the development and discovery of new materials for photovoltaic modules. Another aspect 
that has been the subject of studies is the degradation of photovoltaic modules by several 
factors, compromising performance, reliability, as well as service life. The degradation is 
related to several factors, being able to originate in the manufacturing process of the modules 
and also the exposure conditions, due to the incidence of the radiation itself, the ultraviolet 
radiation, temperature and humidity since they are exposed to different environmental and 
climatic conditions. 

The most common problems in modules with crystalline silicon solar cells are interrupted 
interconnections, broken cells, corrosion, delamination of the encapsulant, discoloration of 
the encapsulant, broken glass, bypass diode and weld tape failure, as well as degradation 
induced by electrical potential, hot spots and bubble formation on the backsheet, among 
others (NDIAYE, 2013). 

A degraded photovoltaic module can continue to perform its main function, which is to 
generate electricity from sunlight, even if its use is not as efficient as initially was. However, 
the degraded state of the module can be more problematic when the degradation exceeds a 
critical point (CHARKI,2013). Itis important to note that manufacturers consider a PV module 
to be degraded when its energy reaches a level below 80% of its initial energy. 

Among the many problems affecting the performance of the PV module, potential- 
induced degradation (PID) has attracted the most interest in recent years, more precisely from 
2010, this is due to the fact that PID occurs frequently in photovoltaic systems and solar parks, 
where the voltage of strings are higher. 

The PID is related to the difference in potential of the module in relation to earth 
combined with high humidity and high temperatures, causing loss in power due to unwanted 
leakage current to the ground, since this decreases the efficiency of the cells (FIGUEIREDO, 
2015). 

Although most studies on potential-induced degradation are based on p-type crystalline 
silicon photovoltaic modules, recent research seeks to include photovoltaic modules 
composed of other cellular materials, as perovskite, thin film solar cells of Indium Gallium 
Copper Disselenete (CIGS) or whose technology is PERC. Thus, this article aims to carry out 
a bibliographic review on potential-induced degradation in c-Si photovoltaic modules. 
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21 ELUCIDATION OF PID IN PV MODULES 


The low maintenance cost, the non-emission of pollutant gases and the lifetime of more 
than 25 years are fundamental requirements for the implementation of photovoltaic technology 
in power plants. However, over this time, some problems begin to arise, whether with the 
inverters installed, whose average operating time is 10 years, whether by other decomposition 
mechanisms that decrease the efficiency of solar generation or contribute to the destruction 
of solar cells. In 2010, potential induced degradation (PID) of crystalline silicon solar cells 
was considered one of the major degradation mechanisms (PINGEL, 2010), this is due to the 
fact that PID causes a serious initial energy loss, negatively impacting the performance and 
reliability of the PV module under field conditions (WANG et al, 2016) . Thus, several global 
PV research institutes have leaned towards PID-related research in photovoltaic modules, in 
which the term PID was coined by Pingel et al, basing a large amount of work in order to solve 
and verify its causes (HYLSKY et al, 2018). 

To achieve greater solar power generation, in the photovoltaic systems connected 
to the grid, PV modules are connected in series to generate a high output voltage, while 
the module frame is grounded for safety reasons (Fig.1) in the meantime solar panels may 
be exposed to high voltage between the solar cell and the module structure. This potential 
difference causes leakage currents to travel between the metallic structure of the module to 
the solar cells, or vice versa, depending on the position of the module in the string, through 
the front glass and the encapsulation. The effect is stronger the closer the module is to the 


negative pole of the photovoltaic system. 





Figure 1. A simplified schematic diagram of a PV system with a floating potential. 


Fonte: LUO, 2016. 


Studies on the degradation mechanisms caused by high voltages and leakage 
currents of photovoltaic modules were initiated in 1985 by the Jet Propulsion Laboratory 
(JPL), more precisely in crystalline silicon (c-Si) and amorphous silicon (a-Si) modules thin 
film, these were performed outdoors and involved electrochemical and galvanic corrosion, 
electromigration, and heating by hot spots. Due to the different environmental and climate 
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conditions and the diversity of solar panel specifications and technologies, there is also 
emphasis on research to mitigate the effecis and causes of PID and to clarify the relative 
sensitivities of various commercial absorbers today (HYLSKY et al, 2018). Since there is no 
constant pattern for the emergence of the effects of degradation on the module. 

The conditions necessary for the occurrence of PID involve (i) environmental factors 
as well as factors on (li) system, (iii) module and (iv) cell level (PINGEL et al, 2010). While 
the environment is set for each individual installation, it is possible to prevent the PID by 
properly controlling only one of the factors (ii), (iii) or (iv) (SCHUTZE et al, 2011). In this way, 
all these factors will be addressed below. 

Humidity and temperature negatively affect the performance of the photovoltaic 
system. In humid climate, the unwanted leakage current decreases the cell performance 
(SCHWARK et al, 2013), since the high humidity and high temperature facilitate the reaction 
and diffusion of the metal ions in the encapsulated glass, since this moisture forms a film on 
the glass, which becomes electrically conductive and therefore the PID effect occurs more 
easily. 

One of the most favorable circumstances for the occurrence of PID is the sunrise 
of sunny days, period of occurrence of dew, due to the formation of droplets forming a thin 
layer of water on the glass that causes the increase of humidity. This increase causes the 
resistivity to decrease and facilitates current leakage through the glass. Throughout the 
day the humidity decreases due to the temperature rise, thus causing a linearly decreasing 
leakage and the PID effect does not become relevant during the generation process. It 
should be noted that the moisture in the encapsulation causes a constant current for a while. 
At the end of the day this current reaches almost zero, but after this time, the temperature 
decreases and the PV module becomes moist again (HOFFMANN, 2014). 

Itis importantto note that just as high temperatures cause an increase in the degrading 
effects caused by PID, high temperatures are used during the regeneration processes of 
modules against PID. As it is not possible to change the environmental conditions in which 
the photovoltaic plants are located, the research is based on examining and addressing the 
PID through the system and the photovoltaic modules. 

Regarding the system level, a difference in potential and signal from the module 
are the aspects that provide the most significant impacis. These in turn depend both on 
the position of the module in the array, the topology of the grounding, the extent of the 
photovoltaic arrangement and the type of inverter connected to the system. 

The choice of the protective glass and the encapsulating material have been 
shown to influence the occurrence of PID, since they are in direct contact with the cells. 
The various types of encapsulation significantly affect PID, modules using EVA (ethylene 
vinyl acetate) are more prone to PID occurrence, since all different substitutes were able to 
prevent. Studies point out that acetic acid contained in EVA together with moisture may be 
responsible for dissolving metal ions at the glass interface (SCHWARK et al, 2013). 
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Base resistivity, emitter sheet strength and anti-reflective coating (ARC) properties 
are parameters that influence the sensitivity of the PID effect at the cell level. The anti- 
reflective coating has the function of increasing light capture, hence the module power 
conversion is higher. However, the properties used in the anti-reflective coating of the solar 
cell are factors that significantly affect the effect due to its dependence on the silicon nitride 
ARC (Sin) (NAGEL, 2011). 

Another factor contributing to the emergence of the PID effect is the accumulation of 
ions in the ARC, since this leads to the formation of an electric field in this region, causing 
an antipassivation effect, increasing the recombination near the surface of the photovoltaic 
cell, in addition to acting as an undesirable path to escape currents (LAUSCH et al, 2014). 


31 DETECTION AND REVERSAL TECHNIQUES 


In order to prevent or mitigate PID effects, it is critical that its detection is done in 
a rapidly and assertively manner, because the degradation mechanisms can result in high 
yield losses of 20%, or more in cells that aren't PID resistance. To prevent degradation, all 
the componenis of the module such as encapsulant and glass should be resistant to PID, 
the solar cell also must be resistant and to achieve this, a possibility is to optimize the anti- 
reflexive coating (KAPUR, 2015). 

The phenomenon occurs most frequently in photovoltaic modules that are closer to the 
negative pole, when C-Si type p cells are used and in high-voltage array systems (KAPUR, 
2015) and various research shows that the type of encapsulant used in the module can directly 
affect the effect (HOFFMANN, 2014). High yield losses may be an indication of degradation by 
the PID effect (HACKE et al, 2011). 

PID can be detected in a number of ways, with thermographic imaging (CÂMARA, 
2019), electroluminescence imaging, open circuit voltage measurement, plotting the |-V 
curves of the modules (ISLAM, 2018), shunt resistance measurement (SCHWARK et al, 
2013), and the use of the mass spectrometer to observe the sodium migration within the 
module (ISLAM, 2018). 

Electroluminescence imaging (Fig. 2) has gained a growing interest due to its quick 
resolved defect detection possibilities such as, detection of recombination, resistive and 
optical losses of the modules, low diffusion lengths, high series resistance and shuntis. It has 
been observed from the electroluminescence imaging of the on-site degraded photovoltaic 
module that a performance degradation happens due to different types of cell defects, such 
as, localized shunting, cracks and front contact grid interruptions (NAUMANN et al, 2013). 
To capture the electroluminescence images the module is operated as a light emitting diode, 
and due to recombination effects the emitted radiation can be detected with a camera that 
photographs the emitted photons. In the images the damaged areas have less or no brightness, 
while those without defects have a shine (SCHWARK et al, 2013). Modules affected by PID 
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develop a lower operating voltage compared to other modules in the array, and these voltage 
differences can cause hot spots. This way the thermographic imaging assists in the detection 
of PID as the hot spots indicate the temperature increase of the affected cells (HACKE et al, 
2011). As we can see in the the image (Fig. 3), the cells closer to the frame show hotter spots 
than the rest of the cell (CÂMARA, 2019). 


SDECORaE::: 
CARRRRRES 


Figure 2. The left image is of the module before degradation and the right image after degradation. 


Fonte: (HOFFMANN, 2014). 








Figure 3. Thermographic image of a module with PID traces. 


Fonte: (OH, 2015). 


When comparing electroluminescence images and electroluminescence images, there 
is a slight correlation between the regions with the highest temperature and those with the 
lowest brightness (PINGEL, 2012). 

The shunt resistance (Rsh) of the modules who has the PID resistance encapsulant 
proved to be approximately steady, while the ones that aren't resistant to PID had shown 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 8 HE 


a decrease in the shunt resistance (ISLAM, 2018). Such resistance is also affected by the 
different voltages of the degraded modules (PINGEL, 2012). 

As the migration of sodium from the anti-reflective coating to the solar cell has been 
observed after PID tests, the use of the mass spectrometer will observe its accumulation in 
both the encapsulant and the solar cell (NAUMANN et al, 2013). The module parameters as 
maximum power point, open circuit voltage, short circuit current, FF (fill factor), and leakage 
current are drastically reduced as both decrease when the module is affect by PID (HACKE 
et al, 2011). 

The |I-V curve (Fig.4) of the affected modules show that the maximum voltage is not 
the same when compared with the modules that did not suffer from the effect (PINGEL, 2012). 
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Figure 4. Curves of modules with different back sheets before and after 24h PID test. 


Fonte: (KOCH et al, 2012). 


Photovoltaic modules shall always be connected in series, taking into account the 
technical characteristics of the inverters used. The effect that this event brings to large 
projects is the high voltage created, causing problems detected at both the system level, the 
panel level and the cell level (HYLSKY et al, 2018). These effects are the major causes that 
can be minimized or even reversible. Climatic factors such as temperature are allies against 
PID, it has been shown that panels stored at temperature around 100 degrees Celsius for 
10 hours lead to a recovery of close to 100% (Fig. 5) (PINGEL, 2012). 

In 2010, the most important factor for the appearance of PID was characterized 
by the large potential difference between cell and soil at the level of the system. For 
this situation, a reverse potential must be acting in order to reduce this DDP (Potential 
Difference) (PINGEL, 2010), a system contained in the inverter referring to the string of 
solar plates called PID doctors can be used for this type of solution. This adaptation will 
reverse the polarity referring to the polarity of energy generation, causing the reduction 
of the generated potential, without this electric field, the sodium ions that were previously 
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deposited in the stacking failures will be removed (HYLSKY et al, 2018). Field research 


for module regeneration shows reversals performed in two distinct strings. In one string a 


grounding kit was used, in another, was used a SMA PV Offset Box that applies a positive 


potential in the panels. In both cases, significant recovery was demonstrated over time (Fig. 


6) (OH, 2015). 
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Figure 5. PID recovery using a temperature of -100ºC. 
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Fonte: (PINGEL, 2010). 
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Figure 6. Relative recovery of time degraded panels versus PVO Box and PV-Grounding. 


Fonte: (SUGIMURA, 1999). 


Reverse voltage system techniques have PV recovery rates (80% - 96%) (CÂMARA, 


2019). It should be noted that the junction between reverse voltage, along with application of a 
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temperature higher than the ambient temperature, can recover the maximum power of a cell by 
up to 96% considering an irradiance of 1000 W/m? (HACKE et al, 2011). Although temperature 
is an important catalyst in panel regeneration, it is evident that very high temperatures can 
affect module materials, impairing their long-term performance (PINGEL, 2010). 

Regarding the panel, inthe manufacturing some presente characteristics of the materials 
that compose it can be modified in order to minimize the intensity of leakage currentis that can 
n case they lead the cells suffer the PID process. In order to prevent the entry of sodium ions 
into the cell and encapsulationççof the panel, a layer of phosphorus silicate (PSG) was used 
in the glass. This layer was elaborated in the solar cell emitter diffusion process (HYLSKY et 
al, 2018), and the goal was not to reduce the efficiency of the photovoltaic cell, parameters 
such as short circuit current and open circuit voltage, maintaining resistive properties against 
the PID, preventing the penetration of sodium on the surface of the silicon. The technique 
was achieved by making necessary modifications to the dissemination process (HYLSKY 
et al, 2018). The transport of sodium from the glass to the cell is one of the major factors for 
the formation of PID, however, it is not yet known for sure how much sodium is needed for 
the formation of PID (OH, 2015). Using a glass that does not have sodium in its composition 
effectively protects the cell against PID (KOCH et al, 2012). Two glasses with different sodium 
levels were tested, borosilicate Glass with 6.5% sodium by weight and soda-Lime Glass with 
16% sodium by weight. It was observed in the PID tester that the glass with a lower sodium 
composition has a higher resistance against PID after 100 hours of testing, while the soda- 
Lime Glass in 24 hours already presented accumulation of sodium deviation, indicating that 
the sodium in the glass is critical for PID (HACKE et al, 2011). 

Changes in composite materials in the separation of photovoltaic panel materials 
demonstrate improvements in the PID process. In (HYLSKY etal, 2018) it proves that lamination 
films intended to prevent thermal and electrical isolation of the cell that directly affecis PID 
degradation. Ethylene Vinyl Acetate (EVA), which is widely used today in the manufacture of 
panels, can be replaced by other materials that have different characteristics ffom EVA, such 
as conductivity, decreasing the migration of sodium ions, slowing the emergence of PID. The 
materials can be partially neutralised by ethylene methacrylic copolymer, EVA copolymer or 
polyolefin elastomer (HYLSKY et al, 2018). Power comparison between modules with ionomer 
and EVA films was diagnosed. Modules with ionomers retain, after 500 hours of exposure to 
PID, 99% power, while modules with standard EVA film lose almost all power in 24 hours of 
testing. A PID tester demonstrates in temperature conditions 85ºC and 85% relative humidity 
and no leaf performed at 96 and 192 hours, stating (Fig. 7) that the standard EVA film has 
low resistance falling to close to 11.5% of its power compared to the ionomer that maintained 
98,5% of its power. Through Laser Ablation-Inductively Coupled Plasma Mass Spectrometry, 
the accumulation of sodium ions can be obtained after the test. The reading demonstrates 
a high concentration of sodium in the EVA encapsulation and in the cell compared to the 
ionomer, reporting also that the deposition of sodium ions in the ionomer layer and cell occurs 
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slowly when compared to standard EVA (HOFFMANN, 2014). 
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Figure 7. Modules with ionomer and EVA demonstrating the resistance of both films against PID. 


Fonte: (KAPUR, 2015). 


At the cell level, modifications to the cell structure can decrease PID aggressiveness. 
In (OH, 2015), it describes the manufacture of two or three anti-reflective coatings to reduce 
the susceptibility of the PID (KOCH et al, 2012). However, the improvements achieved are too 
low to be practical PID solutions, and the disadvantages of this modification can considerably 
affect the efficiency of photovoltaic cells due to the refractive index of these layers. The 
refractive index of ARC was considered the only parameter capable of preventing the PID 
effect in PV modules. According to reports by (SCHUTZE et al, 2011) the PID does not occur 
If the refractive index of a silicon nitride ARC is 2,2% or more. 


41 CONCLUSION 


This article presented a bibliographic review of the Potential Induced Degradation 
(PID) mechanism that is becoming more relevant with the growth of the photovoltaic system, 
in which the system tensions are higher. It has been shown that the origin of the PID can be at 
the cellular, module and system level, as well as the reversal. The solution at the system level 
is to apply a reverse voltage in order to reduce the DDP between the ground and the system. 
At the cellular level many parameters influence the stability of the PID, in addition to the 
resistivity of the base and the resistance of the emitter sheet, the most important parameters 
were the properties of the anti-reflective coating. Since the changes of this layer can avoid 
the PID effect. Recently, most of the modules marketed are already manufactured with the 
PIiD-free description, showing that the industry has been trying to fix the problem, early in the 
manufacture of the modules (FIGUEIREDO, 2017). The high temperatures also proved to be 
an ally to the regeneration of the PV module detected with PID. 

It has also been reported on the detection techniques, in which the images obtained 
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through the thermographic technique and electroluminescence and the plotting of the |-V 
curve, among other methods, are used in several researches, Thermal images are a viable 
technique to be performed in the field and with great possibilities to infer the existence of PID, 
being the rapid possibility of detecting problems in one or more modules. 

Given these conclusions, it is expected that from this study, the understanding of 
the current state of research related to PID, other works could be developed to achieve a 
significant improvement in the effects of the PID and to provide stability of solar panels. 

We intend to continue this work by analyzing the methodologies and techniques 
presented here in photovoltaic modules installed in photovoltaic plants in the region of Campos 
dos Goytacazes (RJ) and Cabo Frio (RJ) in order to prove the analysis of field studies by other 
authors. 
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ABSTRACT: This chapter presents an article 
whose subject seeks to define the optimal 
allocation of capacitors for reactive power flow 
planning in distribution power system. The 
problem is based on a mathematical model that 
minimizes the operation and investment cost, 
where a sensitivity index and the lowest voltage 
are used to select the buses for allocation. A 
Particle Swarm Optimization (PSO) algorithm 
adapted to a discrete version is used to find 
the best quantity and location of the capacitor 
bank (CB). As a result, there was a reduction in 
power loss and an improvement in the voltage 
profile even with the addition of CB. The method 
is applied in the 10 and 34 bus system. The 
proposed algorithm presented a good solution 
compared to other optimization techniques 
presented in the literature. 

KEYWORDS: Power systems, Power distribution, 
Reactive Power Planning. 


Coleção desafios das engenharias: Engenharia elétrica 2 


ALOCAÇÃO DE BANCO DE 
CAPACITORES EM SISTEMAS DE 
DISTRIBUIÇÃO USANDO ALGORITMO 
DE PSO DISCRETO 


RESUMO: Este capítulo apresenta um artigo 
cujo assunto procura por definir a alocação ótima 
de capacitores para o planejamento do fluxo de 
potência reativa no sistema de distribuição de 
energia. O problema é baseado em um modelo 
matemático que minimiza o custo de operação 
e investimento, onde um índice de sensibilidade 
e as barras de menor tensão são utilizados 
para selecionar as barras a serem alocadas. É 
utilizado um algoritmo de Otimização de Enxame 
de Partículas (PSO) adaptado para uma versão 
discreta para encontrar a melhor quantidade 
e localização do banco de capacitores (CB). 
Como resultado, houve uma redução na perda 
de potência e uma melhora no perfil de tensão 
mesmo com a adição do CB. O método é 
aplicado no sistema de 10 e 34 barras. O 
algoritmo proposto apresentou boa solução em 
comparação com outras técnicas de otimização 
apresentadas na literatura. 

PALAVRAS-CHAVE: Sistemas de Potência, 
Sistemas de Distribuição, Planejamento de 
Potência Reativa. 


11 INTRODUCTION 


The power grid is a set of equipment 
responsible for transporting electricity from 
one point to another, to ensure that consumers 
have access to electricity as pleased, and it is 
managed by electrical utility distribution. Besides 
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continuous power supply, quality assurance is in the utility purview, meaning that the voltage 
levels shall be within limits established in technical standards for operating a distribution 
system (DS) (WILLIS, 2004). 

During the years of operation, it is not possible to accurately predict the growth of the 
consumed load by the system, which leads to a reduction in the voltage profile of the grid and 
an increase in the power losses of electrical system (PEREIRA JUNIOR, 2009). 

The DS is highly responsible for power losses in electrical power systems, and it has 
a low-voltage regulation. This issue is within the context of DS operation planning, which 
seeks to encounter a preferable operating scenario, aiming at minimizing costs. In this regard, 
several strategies are researched that reduces power loss in the grid, such as the use of 
distributed generation, the allocation of voltage regulators and the setup of capacitor banks 
(CB) (PRASAD REDDY P, 2014). The use of CB is a technique inserted within reactive power 
planning and usually employed to solve these problems (PRASAD REDDY P, 2014). 

The allocation of CB is justified by the fact that these devices improve the grid voltage 
profile, leading to a more horizontal profile and thus improving the quality of the energy 
delivered to the consumer. Besides that, as they are inserting reactive power into the grid, 
they reduce the apparent power flow in the system, and this leads to lower electrical losses 
(TABARES, 2016). However, mistaken the CB insertion location as well as its quantity may 
worsen the system voltage profile and increase active power losses (DIXIT, 2016). 

Thus, CB allocation can be handled as a combinatorial optimization problem, 
where there is a finite yet very large set of feasible solutions. For such problems, computer 
optimization algorithms are often employed, such as metaheuristics that do not guarantee a 
global solution to the problem, however, are efficient in finding high quality solutions without a 
high computational expense (DOS SANTOS PEREIRA, 2018). 

For this problem, the literature shows a variety of algorithms that can be applied, such 
as Plant Growth Simulation Algorithm (RAO, 2008), Particle Swarm Optimization (DIXIT, 
2016), Swarm robotics search & rescue algorithm, Genetic Algorithm (KALANTARI, 2011), 
Evolutionary Algorithms, Cuckoo Search Algorithm (PRASAD REDDY P, 2014), Gravitational 
Search Algorithm, Whale Optimization Algorithm (NEAGU, 2017), among others. The 
metaheuristic technique employed in this paper was the Particle Swarm Optimization (PSO), 
however, due to the nature of the problem, this algorithm has been adapted to work with 
discrete variables. 

On the basis thereof, this paper presents an algorithm capable of reactive power 
planning of a power distribution system, through CB allocation. The algorithm is guided 
through an objective function that seeks to reduce the active power losses of the system, 
combined with a lower investment cost of CB setup. 

The descriptions of mathematical model employed, the metaheuristic technique for 
objective function minimization, the adaptation of the technique to solve the problem, the 
results obtained, and the conclusion are described in the next sections. 
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21 MATERIALS AND METHODS 


This section describes in detail the mathematical model, aspects, formulation, 
adaptation of the meta- heuristic technique employed and a flowchart of the proposed 
algorithm. The mathematical formulation will be presented below. 

The objective function evaluation that guides the PSO algorithm in search of the best 
location to insert CB is based on the following mathematical model. 


min Z = KpPp + DiB KciQc; 


(1) 

—dp; =D) 1P;j + Bj En (2) 

—dq; + Nck *Qp = D;-10:/+0ji (ENKEN, (3) 
Pi; =V2gij— ViV;gi; *cos(0; — 0;) — V;V;b;; * sen(0; — 0; (4) 
Qi; = —VÊb;; — ViV;gi; * cos(9; — 0;) + V;V;b; ; * sen(0; — 0;) (3) 
Sij S Smax (6) 

Vuin < Vi <Vuax (7) 


The first equation represents the function whose goal is to minimize two cost terms: the 
first of these is related to the system operation, which seeks to minimize the costs with active 
power losses, the second refers to investments in equipment where it looks for to minimize 
the cost of setting-up CB. In the equation, K, ($/kKW) and K, ($/kvar) represents the cost for 
active power loss in the system and the cost of setting- up capacitor banks on bus i, P. (kW) 
represents the total value of active power losses, Q.. (kvar) represents the reactive power of 
CB installed on bus i, and NB represents the set of system buses that has a CB connected 
to them. 

The objective function is subject to a set of equality and inequality constraints. Within 
the equality constrainis, (2) and (3) represenis the system active and reactive power balance, 
where dpi and dgi represents the active and reactive power demands on bus |, respectively, 
P,; and Qu; represents the active and reactive power flow that flows from bus / to bus j, 
respectively, Nrefers to the set of system buses, N,, represents the number of CB installed 
on bus kand Q,indicates the base power each bank can inject into the system. 

Expressions (4) and (5) represents the equation of active and reactive power usage in 
line, respectively, V, and O, represents the magnitude and the angle of bus voltage í and, g,, 
and b ; fepresents the conductance and susceptance of line connecting bus ito bus j. 

Regarding inequality constraints, (6) represents the apparent power limit on the line, 
which cannot exceed a maximum value, and (7) indicates the minimum and maximum voltage 
level limits on buses. 

The optimization technique used is based on the particle cooperative movement within 
a swarm, like bird's flights, seeking the best solution for optimization problems (KENNEDY, 
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1995). The particles of PSO algorithm moves in three terms: the inertia term, which quantifies 
the particles tendency not to change direction, that is, keep its own movement, the cognitive 
term, which quantifies the particle tendency to follow the best path it has ever found, and the 
social learning term, which quantifies the particle tendency to follow the best path the entire 
swarm has already found, that is, the best position found. 

Each particle carries two pieces of information: the number of CB to be installed and 
the candidate bus for receiving the devices, thereby defining the solution. Thus, at each 
iteration, the particles swarm moves within the solution space, seeking the best possible 
position. The position of each particle and the speed with which it moves are represented 
through expression (8). 


Vit sã wvj +cl+ fpoe — pi) + C2 + (Ibest - pp) (8) 


In (8), the first part refers to the inertia term, where vi indicates the velocity of particle 
p in iteration / and wis the particle inertia index. The central expression indicates the cognitive 
learning term, where c1 is the cognitive learning index, pp“! is the best position where 
the particle p has ever been, ph is the position of particle p in iteration i. And the last part 
shows the global learning term, where c2 is the global learning index and Gpest is the best 
position, any particle has ever been. The position of each particle represents the number of 
CB connected to bus 1, where each CB is capable to provide a total of 150 kvar. 

The PSO algorithm was designed to solve problems with continuous variables, 
however, the goal of this problem is to establish the number of CB that will be connected in 
each bus. Thus, the algorithm required adaptation so that it could work with discrete variables. 
The authors (KENNEDY, 1997) who initially developed PSO, also presents a way to discretize 


variables in PSO algorithms, which is based on the formulation presented below. 


f=""A (9) 


1+e | 


Pp = (10) 


REÊ o f >rand 
0, f<rand 


Equation (9) represents the Sigmoid Function which is responsible for establishing 
the particles velocity between O and 1, thus, equation (10) performs a draw to determine the 
logical level that represents that speed. The relation of this coding to the position vector of 
each particle is illustrated below. 

10 1 
pel 2 34 s1=/o 1 ] (11) 
00 1 

Equation (11) shows that each particle position is now represented by a binary vector, 

to use integer variables in discrete optimization problems. In addition, the position speed of 


each particle is also represented by a vector, so expressions (9) and (10) are used for each 


element of position matrix. 
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The definition of which system buses are candidates for insertion of CB is guided by 
two indicators presented below. 


20ijTij 
LSF, = SE (12) 
Norm; = a (13 


Equation (12) indicates the Losses Sensitivity Factor on bus | (LSF.), where Fis 
represents the resistance of the line connecting bus to bus, and (13) represents a vector that 
normalizes voltage values on the buses. Both equations represent the criteria used to select 
the candidate buses for the CB installation. Firstly, the LSF is calculated and descending 
ordered, then are selected the buses whose corresponding normalization vector is greater 
than 1.01. The quantity definition and which are the candidate buses for CB allocation helps 
reduce the solving problem space (RAÃO, 2008). 

A schematic diagram with the running steps of the algorithm is illustrated in Figure 
1, which shows the execution process of the algorithm. The position of each particle is used 
to evaluate the objective function, so these values are compared to the best positions each 
particle has ever been, and with the best position the whole swarm has ever been. If these 
positions lead to a better value of the objective function, its stored. The best positions of each 
particle are used in equation (8) for updating the velocity of each particle. 


31 RESULTS AND DISCUSSIONS 


The proposed algorithm was evaluated in test systems presenis in literature, and the 
results compared with those obtained by other authors. For the algorithm execution, algebraic 
modeling software GAMS was used, with the CONOPT solver. In both cases, 10 particles, 50 
iterations, K, = 168 S/kW, w=cl =c2 =1 were used, and the substation bus was chosen as 
a system voltage reference with value of 12..0º [p.ul. 
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Figure 1 - Algorithm Flowchart. 


Table 1 shows the available CBs that were considered in this test. The algorithm was 
written such manner that CBs with best cost-benefit ratio were used first, that is, the smallest 
values from the third column of Table 1. 


Power (kvar) Cost ($) Ratio ($/kvar) 
1200 2040 ERRA 
900 1650 1.833 
600 1320 2.2 
450 1140 2.533 
300 975 3.25 
150 750 5 


Table 1 - Available Capacitor. 
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The 10-bus system operates at 23 kV and has a total load of 12368 kW and 4186 kvar. 
System data is available at (HAMADA, 2008). Figure 2 illustrates the system configuration. 
Table 2 indicates which system buses were considered in the optimization process, as well as 
the amount of CB to be installed. Table 3 shows a comparison between the base case where 
there is no CB installed, the results obtained by (DIXIT, 2016) and (HAMADA, 2008) and those 
found in this test. Figure 3 shows a comparison between system voltage levels before and 


after CB allocation. 


l 2 3 - 5 6 1 8 9 10 


Figure 2 - 10-Bus System Configuration. 


Results 
Candidate B04 BO5 BO6 BO9 B10 
Buses 
Installed Banks 10 12 (6) 0 2 
Table 2 - Installed Capacitors 10-Bus System. 
Base Case (DIXIT, 2016) (HAMADA, 2008) Discrete PSO 
Power Loss (kW) 783.78 703.97 684 081.07 
Operation Cost ($) 131674.9 118226.96 114912 114503.76 
FeRONtE OR 0 2850 5400 A800 
(kvar) 
died ROSA 0 588 199.4 939 
Total Cost ($) 131674.9 118854.96 116111.4 115441.98 
VMIN 0.838 0.873 0.9 0.878 
VMAX 0.993 0.995 - 0.997 


Table 3 Comparison of Obtained Results 10-Bus System. 


According to Table 3 it is possible to notice that there was a 13% reduction in active 
power losses decreasing as at 102.21 kW. In addition, the total reactive power shows that a 
large amount of CB are not required, but this amount must be allied to the smart choice of 
installation buses. 

The 34-bus system operates at 11 kV and has a total load of 4636.5 kW and 2873.5 kvar. 
System data is available at (HAMADA, 2008). Figure 4 illustrates the system configuration. 
Table 4 shows the CBs installed on the system buses, as well as the considered options. 
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Table 5 shows a comparison between the base case, where no CB is installed, the results 
obtained by (DIXIT, 2016) and (HAMADA, 2008) and those found in this paper. Figure 5 shows 
a comparison between system voltage levels before and after CB allocation. 


Voltage Profile - 10 Bus System 


Before Allocation 
After Allocation 





E 
E 
Figure 3 Voltage Profile 10-Bu s System. 
Results 
Candidate Buses B19 B20 B21 B22 B23 
Installed Banks 6 0 4 1 0 
Candidate Buses B24 B25 B26 B27 - 
Installed Banks 0 1 O 1 - 
Table 4 - Installed Capacitors 34-Bus System. 
Base Case (DIXIT, 2016) (HAMADA, 2008) Discrete PSO 
Power Loss (kW) 221.71 168.37 158 166.94 
Operation Cost ($) 37249.19 28286.16 26544 28045.36 
Reactive Tot. (kvar) 0 1950 3000 2100 
Installation Cost ($) 0 Se ls 1365 544.5 
Total Cost ($) 27249.19 28807.86 27909 28589.86 
VMIN 0.9417 0.9502 0.951 0.9504 
VMAX 0.9941 0.9950 - 0.9950 


Table 5 - Comparison of Obtained Results 34-Bus System. 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 9 





Figure 4 - 34-Bus System Configuration. 


Voltage Profile - 34 Bus System 


Before Allocation 


After Allocation 


Voltage [p.u.) 





Bus 


Figure 5 - Voltage Profile 34-Bus System. 


Analyzing Table 5 it was observed that there was a reduction of 24.7%, that is, 54.78 
kW in active power losses. In addition, the total determined reactive power shows, as well 
as the 10-bus system, that a large amount of CB is not required, but that amount must be 


consistent with the smart choice of installation buses. 
Figure 5, as well as Figure 3, shows that CB allocation, furthermore, to reducing active 


power losses, leads to an improvement to grid voltage profile, which means a better operating 


point for the system. 
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41 CONCLUSION 


This chapter presents an article that presents the solution of a cost minimization 
problem of a DS for the optimal allocation of CB using adapted PSO metaheuristic technique 
for discrete variable problems. The algorithm used was evaluated in tests systems presents 
in literature and compared with other optimization techniques. In general, the discrete PSO 
algorithm showed itself robust, presenting excellent performance. Furthermore, it was efficient 
in determining the quantity and location of installation of the CB on the electrical grid, with 


low total costs, reducing operation and investment costs and minimizing active power losses. 


5 | ACKNOWLEDGEMENTS 


The author would like to thank CAPES for the financial support for the development of 
this research. 


REFERÊNCIAS 


DIXIT, M., KUNDU, P., JARIWALA, H. R. Optimal allocation and sizing of shunt capacitor in 
distribution system for power loss minimization. /n: IEEE Student's Conference on Electrical, 
Electronics and Computer Science (SCEECS), IEEE, 2006, p. 1-6. 


DOS SANTOS PEREIRA, G. M., FERNANDES, T. S., AOKI, A. R. Allocation of Capacitors and 
Voltage Regulators in Three-Phase Distribution Network. Journal of Control, Automation and 
Electrical System, v. 29, n. 2, p.238-249, 2018. 


HAMADA, M. M., WAHAB, M. A., EL-SAYED, A. H. M., RAMADAN, H. A. A proposed strategy for 
capacitor allocation in radial distribution feeders. In: 12th International Middle-East Power System 
Conference, IEEE, 2008, p. 146-151. 


KALANTARI, M., KAZEMI, A. Placement of distributed generation unit and capacitor allocation in 
distribution systems using genetic algorithm. /n: 10th International Conference on Environment and 
Electrical Engineering (EEEIC), Conference Proceedings n. 2, IEEE, 2011, p. 1-5. 


KENNEDY, J., EBERHART, R. Particle Swarm Optimization. /n: International Conference on Neural 
Network (ICNN), Proceedings of ICNN'95, 1995, p. 1942-1848. 


KENNEDY, J., EBERHART, R. Discrete binary version of the particle swarm algorithm. /n: IEEE 
International Conference on System, Man, and Cybernetics, Computational Cybernetics and Simulation, 
v. 5, 1997, p. 4104-4108. 


NEAGU, B. C., IVANOV, O., GAVRILA, M. Voltage profile improvement in distribution networks 
using the whale optimization algorithm. /n: 9th International Conference on Electronics, Computers 
and Artificial Intelligence, Conference Proceedings, 2017, p.1-6. 


PEREIRA JUNIOR, Benvindo Rodrigues. Alocação Otimizada de Banco de Capacitores em 
Sistemas de Distribuição de Energia Elétrica através de Metaheurísticas Multiobjetivo. 2009. 
Master's thesis (Master in Electrical Engineering) — Faculdade de Engenharia de Ilha Solteira, 
Universidade Estadual Paulista, Ilha Solteira, 2009. 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 9 cus 


PRASAD REDDY P, D. Sensitivity Based Capacitor Placement Using Cuckoo Search Algorithm 
for Maximum Annual Savings. IOSR Journal of Engineering, v. 4, n 4, p.06-09, 2014. 


RAO, R., NARASIMHAM, S., RAMALIMGARAJU, M. Optimal capacitor placement in a radial 
distribution system using plant growth simulation algorithm. International Journal of Electrical 
Power & Energy Systems, v. 33, n. 2, p.133-139, 2008. 


TABARES, A., FRANCO, J. F., LAVORATO, M., RIDER, M. J. Multistage Long-Term Expansion 
Planning of Multiple Alternatives. IEEE Transactions on Power Systems, v. 31, n. 3, p.1900-1914, 
2016. 


WILLIS, H. L. Power Distribution Planning Reference Book. New York: Marcel Dekker Inc., 2004. 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 9 ELE 


CAPÍTULO 10 


DESIGN OF A TRANSMISSION-LINE METAMATERIAL 
WITH A NEGATIVE INDEX OF REFRACTION AT 


Data de aceite: 01/10/2021 
Data de submissão: 30/07/2021 


Lucas Douglas Ribeiro 

Pontifical Catholic University of Minas Gerais, 
Graduate Program in Electrical Engineering 
Belo Horizonte, Minas Gerais 
https://orcid.org/0000-0002-0230-8899 


Juscelino Júnior de Oliveira 

Pontifical Catholic University of Minas Gerais, 
Graduate Program in Electrical Engineering 
Belo Horizonte, Minas Gerais 
htitps://orcid.org/0000-0003-2364-4741 


Arthur Henrique de Lima Ferreira 

Pontifical Catholic University of Minas Gerais, 
Graduate Program in Electrical Engineering 
Belo Horizonte, Minas Gerais 
htips://orcid.org/0000-0002-8018-6823 


Rose Mary de Souza Batalha 

Pontifical Catholic University of Minas Gerais, 
Graduate Program in Electrical Engineering 
Belo Horizonte, Minas Gerais 
https://orcid.org/0000-0002-6157-1574 


ABSTRACT: A metamaterial based on a two- 
dimensional transmission-line network loaded 
with inductors and capacitors enabling to achieve 
negative-refractive-index (NRI) is developed. 
The dispersion characteristics are calculated by 
their equivalent circuit model and an operating 
frequency of 2.95 GHz in which there is 
impedance matching with free space is chosen 
in the S-Band. This NRI metamaterial supports 
transverse electric (TE) waves, thus it can be 
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S-BAND 


used in applications such as lensing. 
KEYWORDS: Composite  Right/Left-Handed 
(CRLH); left-handed media (LHM), metamaterials, 
negative permeability, negative permittivity, 
negative refractive index, periodic structures. 


PROJETO DE UM METAMATERIAL 
BASEADO EM PARÂMETROS DE LINHA 
DE TRANSMISSÃO COM ÍNDICE DE 
REFRACÇÃO NEGATIVO NA BANDA S 


RESUMO: Foi desenvolvido um metamaterial 
baseado em parâmetros de linha de transmissão 
bidimensional carregada com indutores e 
capacitores que permite obter índice de 
refração negativo (do inglês, negative refractive 
index, NRI). As características de dispersão 
são calculadas por meio do modelo de circuito 
equivalente e a frequência de operação de 2,95 
GHz na qual há casamento impedância com 
o espaço livre é escolhida na banda S. Este 
metamaterial NRI suporta ondas transversais 
eléctricas (TE), portanto, pode ser usado em 
aplicações como lentes. 


PALAVRAS-CHAVE: Composto Direito/ 
Esquerdo (CRLH); meio Esquerdo (LHM), 
metamateriais, permeabilidade negativa, 


permissividade negativa, índice de refração 
negativo, estruturas periódicas. 


11 INTRODUCTION 


In 1968, Veselago investigated materials 


with 'simultaneously negative | permittivity 
and permeability, due to the reversal of the 


Snell's law, negative phase velocity, and some 
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other electrodynamic properties, they are known as left-handed materials (LHM). These 
metamaterials (MMs) are artificial composites characterized by subwavelength features and 
effective negative value of the refractive index (NRI). LHMs can be applied to optical and 
microwave applications, including beam steerers, modulators, band-pass filters, microwave 
couplers, antenna radomes, and superlenses [1], [2]. As Pendry proposed in 2000 [3], the 
application of MM lensing, which avoids the diffraction limit by utilizing both periodic and 
evanescent electromagnetic waves, further increased the interest for NRI MMs. 

Whereas the split-ring resonators (SRRs) and thin wires were the first realizations 
of LHMs [4]based on a periodic array of interspared conducting nonmagnetic split ring 
resonators and continuous wires, that exhibits a frequency region in the microwave regime 
with simultaneously negative values of effective permeability mu(eff, the first MM lens in planar 
form using the NRI transmission-line (NRI-TL) approach was achieved by [5]. In that work, a 
host TL is loaded using lumped shunt inductors and series capacitors that, respectively, create 
the conditions for an effective NRI. 

The TL MMs can be used in many applications, which include electrically small 
antennas, subwavelength resonators, guided wave applications, phase compensators, LTE 
antennas, radar applications [6], superconductors [7] and perfect lens, and superlens [8] 

Some developmenis of volumetric arrays of NRI-TL metamaterials that interact 
efficiently with free space were achieved in [9], [10] and led to the experimental demonstration 
of free space NRI lens [8]. The free space metamaterial lens development in [11] allows 
electromagnetic waves to focus into a focal region. The designed lens employs a multilayer 
geometry of unit cells that are loaded using fully printed capacitors and inductors in a substrate. 

A similar development is done in this work. We propose to design a transmission-line 
(TL) network periodically loaded using discrete (surface-mounted) inductors and capacitors. 
The TLis modeled by a metamaterial unit cell using circuit theory, and an operating frequency 
of 2.95 GHz at S-Band, in which there is impedance matching with free space, is chosen. 
The unit cell size is approximately 11 times smaller than the wavelength, which is sufficient to 
consider it as an effective medium. The equation for the dispersion relation is compared with 
the simulation results to analyze the metamaterial behavior. 

The behavior of an NRI-TL metamaterial in the GHz frequency is studied in this paper. 
Section Il presents the metamaterials theory based on transmission line parameters. In 
Section Ill, the simulation of a unit cell is shown. The Conclusion is in Section VI. 


21 METAMATERIALS BASED ON TRANSMISSION LINE PARAMETERS 


Metamaterials can be defined as artificial structures of natural elements that when 
correctly arranged can exhibit unnatural behavior. When they have an average cell size much 
smaller than the wavelength, they become effectively homogeneous structures. Thus, the 
refractive phenomenon can overlap the diffraction and scattering in MMs when there is wave 
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propagation. 

It is possible to separate the MMs into two groups: resonant and non-resonant. The 
most common resonant is the Split Ring Resonator. In the group of non-resonant, we have the 
metamaterials modeled by transmission line theory, also called TL metamaterials [12]. 

The non-resonant metamaterial, also found in the literature as CRLH (Composite 
Right/Left-Handed), do not rely on resonances to provide negative parameters, thus obtaining 
lower ohmic losses and larger bandwidth when compared with the resonant metamaterials [5], 
[13]. The TL MMs can be used in two configurations: shunt and series [9], [14]. 

Any natural dielectric can be modeled by an electric circuit composed by a series 
inductance and a shunt capacitance, which are directly related to the material parameters 
permeability and permittivity. This theory, a circuit of unit cell size d with shunt inductance and 
a series capacitance, is used to model a TL metamaterial and obtain negative parameters. 

A traditional transmission line can be used as the host to implement a material with 
negative parameters, together with loading elements as inductors and capacitors. An example 
of this configuration is the microstrip line that can be made in shunt and series configurations 
[9]. There are two co-planar strip lines on the top of a substrate in the series topology, and 
there is no ground plane. In the shunt configuration [15], the stripline is on top of the substrate 
and the ground plane on the bottom. These two configurations are shown in Fig. 1. 





& 


” 
(a) fã 


Fig. 1 - Transmission line host topologies. (a) Series. (b) Shunt [9]. 





The series configuration is shown in Fig. 1 (a). It can be seen as the connection of four 
transfer-matrices, which relates the input and output terminal voltages and currents of a two- 
port network. The connected transfer-matrices are shown in Fig. 2. 
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Fig. 2 - Block diagram for the series unit cell. 


Every arm of Fig.2 is a block represented by its transmission matrix and presents the 
following relations: 


vo e Vie-iBxd To Ei Iie")Bxd 
E am i -JByd O -JByd (1) 
W = We lh = he 


Where B is the phase delay over the cell dimension d. Due to the symmetry of the 
circuit, the transmission matrix is inverted in the output side. 

The series impedance Z = j(w L,-1/(wC,) and the shunt admittance Y=j(wC,-1/ 
(WL,) can compose the circuit of each loaded block, where L, and C, represent the natural 
characteristic of the strip lines, also called host line elements, and L and C, are the loading 
elementis. 


The circuit of the loaded block can be seen in Fig. 3. 





Fig - 3. Circuit model of one block. 


Before obtaining the values of the terms L, e C,, it is necessary to calculate the 
effective permittivity e., of the medium. Considering a stripline with lands of width w 
separated by a distance s, on a substrate with relative permittivity =,, a height h, and 
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surrounded by air with permittivity e,. As discussed in [16], it can be considered that the 


same stripline is surrounded by a single dielectric with an effective permittivity c ,,., as shown 


NS 


o e o qe ap 44" e 





Fig. 4. Effective permittivity [16]. 


Then, effective permittivity can be calculated by: 





h 
Eeff = (ct o, 775 In (— ) +1 75] 


(2) 
Ee > [0,04 — 0,7k + 0,01(1 — 0,15,)(0,25 + o) 


where k=s/s+2W. 


The values of L, and C, can be obtained from: 

















120d 1+ 1 
in(25 dmastsi 
Vo vVk/ 2 
Lp = 377nd 1 
skS-= 
vom (2 ENE) 2 
0 1—k' 
Ferçá Ra] (3) 
1+vk) 2 
120voln | 2 
C 1=vk 
: £ am (285) 
eff — 7 
ia 1 ste 
377nvo 2 


where k' =v1—k? e v is the speed of light in a vacuum. 


The transmission matrix coefficients can be defined as: 


A=ZY+1 B=Z C=Y D=1 (4) 


Through the circuit theory in the diagram of Fig.2, it is possible to get the following 


equations: 
Vi = DV; — BI; = —CV +AL 
V = DV; — BI L = CV +AL 
V; = DV? — BIS kh =CVº+AIS (5) 
V, = DV? — BIS = CV? + Ig 


0O=V+V-—-Va—V, L=L=6b=L 
After solving (1) and (5): 
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vi 
E; 


ed 
1, 
Crie Aep-A 0 0 
0 0 Crie Hey=A 
[M] x Ce; Aex C , y (7) 


D-De, -B-Be, D-De, -B+ Be, 


Where e,= e -/Ê* and e = e -/?Yº. The dispersion relation shows the behavior of MM as a 
function of frequency. This relation can be found by solving: 


det(M) = 0 (8) 


For simplicity, only axial propagation in the x-direction will be considered, so that there is no 
delay in the y-direction and B, After solving (8) and simplifying it, the dispersion relation is: 


Bd = arccos(1 + 4BC) (9) 


31 SIMULATION OF A UNIT CELL 


The software CST Studio Suite? was used to simulate a unit cell, according to the theory 
developed in the previous section. A comparison between the dispersion relation generated 
by theory and simulation is made. The simulation also provides the material parameters that 
can be used in the development of the complete system. 

The values of the host line elements, L, and C,, are dependent on the unit cell 
dimensions and material. The loading elements L, and C can be freely selected. Itis important 
to stick to the matching condition with free space that is achieved when e, =-1 and ur=- 1, 
culminating in a negative refractive index n =-1 and a characteristic impedance Z = 3770. The 
cell size must be small when compared to the wavelength used, in this case, approximately 
11 times smaller. 

The free space matching condition can be achieved by many L and C parameter 
combinations. These values were defined by iterative calculations. In Fig. 5 (a), the top view 
of the simulated unit cell can be seen. The substrate material is Rogers ROS008 with e, = S 
and height h =1.54 mm. 
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Electric wall 

















Magnetic wall bm 


(b) 


9.000 mm 
(a) 


Fig. 5. (a) Top view of the simulated unit cell with main dimensions. (b)Boundaries conditions used in the 
simulation. 


The host line elements were obtained for each block using the main dimensions and 
referring to [16], [17]. Furthermore, there is no ground below the substrate. Then we have 
L,=2.78$nHand C, = 0.16 pF, and the loading elements selected are C = 0.6 pFandL,.To 
reduce reflection, the 90-degree bends were mitered. 

To simulate the system is necessary to impose boundary conditions, as shown in Fig. 
5 (b). These conditions are imposed to simulate the neighboring unit cells by keeping the 
behavior of the field as if the unit cell were in the full model. Through two waveguide poris, 
the unit cell of Fig. 5 was simulated with the propagating wave passing along the x-axis. 

The dispersion relation was drawn and the S-Parameters were retrieved, as shown 
in Fig. 6, along with the dispersion relation calculated with (9). Considering the frequency 
range of the left-handed (LH) region, between 1.85 and 3.39 GHz, the difference between 
the calculated and simulated curves is 11.9%. The LH region shaded in gray is where 
backward waves and negative refractive index occur. The matching point with free space at 
2.95 GHz is indicated by the intersection point with the light line. 
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Fig. 6. Dispersion relation found by theory and by simulation. 


Through the simulation, the unit cell parameters of permittivity =, permeability 4, 
refractive index n and the impedance relative to the characteristic impedance of free space 
Z were extracted, as shown in Fig. 7. The intersection point with the light line indicates the 
matching point frequency with free space at 2.95 GHz given by the dispersion relation of 
Fig. 6. 


Epsilon 1 (Zmax=4.5) : -0.97577318| 


Mu r(Zmax=4,5): -1,0182589 
n(Zmax=4,5): 099777165 
z (Zmax=d,5) : 0,99302029 





| 15 2 2.5 35 4 45 5 


Fig. 7. Unit cell parameters extracted from simulation. 


41 CONCLUSION 


In this work, a transmission-line (TL) network periodically loaded using inductors and capaci- 
tors was modeled by a metamaterial unit cell using circuit theory. The discrete capacitor and 
inductors enabled stronger loading than with fully printed ones. A unit cell size smaller than 
15/11 allowed the use of effective medium theory in this model. The equation for the disper- 
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sion relation was validated by simulations and used to analyze the metamaterial behavior. 
The simulations confirmed the NRI and impedance matching to free space at the operating 
frequency of 2.95 GHz at the S-Band. A structure using the unit cells in the shape of a lens 
with fixed curvature that can focus electromagnetic waves is being developed. The lens will 
have a ladder-shaped surface since it is made with unit cells. 
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RESUMO: Este artigo descreve um modelo de 
perda de propagação para o Sistema Long- 
Term Evolution (LTE), em uma ponte sobre o Rio 
Tocantins na Cidade de Marabá — PÁ. Os dados 
sobre aintensidade do sinalrecebidoforamobtidos 
em campanhas de medição sobre a Ponte. Um 
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sistema Neuro-Fuzzy é usado para desenvolver 
o modelo proposto. O cenário estudado foi em 
uma região com cidade e rio, ambiente bastante 
comum no Norte do Brasil e é pouco estudado. 
Um modelo adequado para a tecnologia LTE é 
de grande importância para o fornecimento de 
um serviço eficiente, especialmente em regiões 
que tem características únicas como no Norte do 
Brasil. Os resultados confirmam a aplicação do 
modelo proposto com erro RMS não superior a 
2 dB. 

PALAVRAS-CHAVE: Perda de propagação; LTE; 
Neuro-Fuzzy; Rádio Propagação sobre pontes. 


RADIO PROPAGATION AND MODELING 
FOR A BRIDGE OVER THE TOCANTINS 
RIVER FOR LTE 


ABSTRACT: This paper describes a model of 
propagation loss for the Long-Term Evolution 
(LTE) system on a bridge over the Tocantins River 
in the city of Marabá-PA. Data on the received 
signal strength is obtained in a measurement 
campaign. A Neuro-Fuzzy system is used to 
develop the proposed model. The studied 
scenario is a region with the city and river, a 
very common environment in northemn Brazil 
and, a few studied. An appropriate propagation 
model for LTE technology is of great importance 
to providing an efficient service, especially in a 
region of unique climate and characteristics. The 
results confirm the applicability of the proposed 
model with RMS error no greater than 2 dB. 
KEYWORDS: Propagation Loss; Long-Term 
Evolution (LTE); Neuro-Fuzzy, Radio Propagation 
Over Bridges, Recovery Effect. 
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11 INTRODUÇÃO 


Um modelo adequado de predição de propagação é essencial para a cobertura 
móvel em cidades localizadas no entorno dos rios amazônicos, no qual o deslocamento de 
pessoas, bens e serviços se dará por meio de pontes trafegáveis sobre regiões hidrográficas. 
A ponte de Marabá sobre o rio Tocantins tem a função de permitir o cruzamento do transporte 
ferroviário pela Estrada de Ferro Carajás e o cruzamento da BR-155, sendo responsável pela 
ligação de cidades do sudeste paraense, como Marabá, até a costa norte brasileira [1]. Logo, 
a mobilidade nesse trajeto se dará em ambiente misto de cidade e rio. Com isso, um modelo 
adequado de predição de propagação de dados móveis, nas estações bases, são cruciais no 
planejamento de rádio propagação para esse cenário [2]. 

O LTE é um dos sistemas de rádio frequência de celular mais importante baseado nos 
padrões desenvolvidos pelo 3GPP (3rd Generation Partnership Project). O sistema foi criado 
com a função de fornecer taxas de dados mais altas para aplicações de dispositivos móveis 
através da rede de rádio [3]. O LTE é também conhecido como a tecnologia de Acesso à 
Rádio Terrestre UMTS Evolved (E-UTRA), tendo como principal objetivo prover maior taxa de 
dados e escalabilidade na largura de banda [4]. 

Com o avanço da tecnologia, o aumento do número de usuários e novos serviços, 
faz-se necessário estudos de predição de propagação em redes sem fio. Tais estudos podem 
permitir o funcionamento adequado de sistemas de celular, sem a interferência entre as 
células que venham a utilizar as mesmas frequências em ambiente que envolva terrenos 
suburbanos e sobre rios. Com isso, as discussões sobre modelos de perda de propagação 
para várias faixas de frequências, tipos de trajetos, correções para perfis do terreno, alturas 
das antenas transmissora e receptora entre outros são frequentes para adequação de novas 
tecnologias de transmissão [5]-[8]. 

A seleção de um modelo de propagação de rádio adequado para a tecnologia LTE é 
de grande importância para os serviços de comunicação e de trafego de dados nas regiões 
do norte do Brasil. No entanto, quando se trata de perda de propagação em ambientes de 
terra e rios, poucos modelos têm sido desenvolvidos como em [6] e [9-10]. Adicionalmente, 
encontram-se poucos estudos de propagação sobre pontes, como o descrito em [11]. 

Um modelo de propagação descreve o comportamento do sinal enquanto ele é 
transmitido e proporciona uma relação entre a distância transmitida pelo transmissor em 
direção ao receptor e a potência recebida. A partir dessa relação, pode-se projetar a perda 
de potência no trajeto caminho e o máximo intervalo entre as células de transmissão. A perda 
de propagação depende da condição do ambiente, frequência de operação e condições 
atmosférica [4]. Poucos trabalhos têm demostrado que o uso de ANFIS (Redes Neuro Fuzzy) 
tem sido eficaz na predição da intensidade do sinal, nas faixas de UHF e VHF [12-13]. Nestes 
trabalhos foram mostraram ambientes outdoor, mas não incluía pontes sobre rio. 

Este trabalho propõe um modelo de propagação em ambientes misto de cidade e rio 
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adequado para o sistema LTE. Os dados são coletados através de uma campanha de medição 
e organizados através de ferramentas computacionais. São apresentadas simulações de 
rede feita com o sistema Neuro-Fuzzy para predição da perda de percurso. 

Este artigo está dividido em cinco seções. Na seção Il está descrita a campanha 
de medição realizada na ponte de Marabá sobre o rio Tocantins. Na seção Ill é exposto 
o desenvolvimento do modelo de propagação utilizando rede Neuro-Fuzzy. Na seção IV 
são apresentados os resultados da campanha de medição e, na seção V, as conclusões e 


perspectivas de trabalhos futuros. 


21 CAMPANHA DE MEDIÇÃO 


a. Cenário 


A Campanha de Medição foi realizada na cidade de Marabá, a qual está situada na 
região do sudeste paraense, a 554 km da capital do Estado [1], dentro da Região Amazônica. 
A cidade possui um clima tropical úmido e se encontra as margens do rio Tocantins. As 
medições se concentram ao longo da ponte que passa sobre o rio Tocantins [14], sendo este 
um trajeto misto composto por estrada sobre a terra com 600 m de comprimento e a ponte 
sobre o rio de 1.75 km de comprimento, conforme mostrado na Fig. 1. O cenário de medição, 
onde vários pontos de recepção foram coletados sobre a ponte é apresentado na Fig. 2. Não 
foram considerados os pontos de recepção perto do transmissor, localizado na cidade, pois 


eles não fazem parte da ponte. 





Figura 1 - Percurso da campanha de medições. 


O cenário de medição, onde vários pontos de recepção foram coletados sobre a ponte 
é apresentado na Fig. 2, em que: D é a distância do receptor e o transmissor; é a distância 


total horizontal e é a distância sobre a terra. 
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Hrx 





Figura 2 - Cenário de medições. 


b. Equipamento e coleta de dados 


A antena transmissora está localizada após a ponte sobre o rio Tocantins, com se 
observa na Fig.2. A altura da antena é de 40 metros, possui um ganho de 15.5 dBi, a potência 
de transmissão é de 40 dBm e a frequência de transmissão é de 783 MHz. Estes dados foram 
fornecidos pela operadora, sendo que o modelo da antena faz parte da confidencialidade de 
dados dela. Não entanto isto não influencia no desenvolvimento do modelo de propagação 
com redes neuro-fuzzy, devido a que os dados necessários para os cálculos foram os antes 
mencionados. 

Na recepção foi utilizado um celular Huawei PLITE 9, o mesmo que possui uma 
antena interna omnidirecional, com ganho de 2 dBi, com polarização circular. Para comprovar 
a polarização foram feitas medidas com o celular na posição horizontal e vertical, e não foram 
observadas mudanças significativas. A unidade de recepção foi colocada dentro do veículo, 
sendo assim a altura do receptor é de 1.28 m em relação do nível do solo. A velocidade média 
do veículo de 40 km/h. O software G-Net Track Pro [15] foi instalado no receptor, para fazer a 


coleta dos níveis da potência de recepção e das coordenadas de latitude e longitude. 


c. Tratamento dos Dados Medidos 


Após a campanha de medição, os dados coletados foram organizados para calcular a 
distância entre o transmissor e o receptor e o cálculo de perda de percurso [16], considerando 
que foi medida a potência de recepção. 

A distância horizontal, entre dois pontos, é calculada usando a fórmula de Haversine 
[17], na qual usa a latitude e longitude entre dois pontos em uma esfera, com raio R (neste caso 
o raio da terra em metros), latitudes Lat, e Lat,, (respectivamente, receptor e transmissor), 
longitudes Long, e Long,, (respectivamente, receptor e transmissor). Os ângulos são em 
radianos e a distância, DIST, entre dois pontos da esfera (neste caso o globo terrestre) 


HAVERSIN ? 
é relacionada as suas localizações. A fórmula (1) demonstra a equação utilizada. 
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Dist. Haversin = 6371000 - acos(cos(90 — Latrx) 
- cos(90 — Latax) 
+ seno(90 — Latx) (1) 
- seno(90 — Latpx) 


-cos(Longax — Longrx)) 


Com a distância horizontal (DIST, 


uaversin) entre a torre de transmissão e os pontos de 


coleta de dados obtidos, foi possível calcular a distância real entre a antena transmissora e a 
antena receptora através do teorema de Pitágoras, como mostrado na Fig. 3 a seguir. 






Distância Real 


Hx = Hx 


Distância Harvesin 


Fig. 3. Cálculo da distância real. 


Para o cálculo da distância real (DIST, vers)» adaptou-se o teorema de Pitágoras, 
através da diferença entre a altura H,, do transmissor, de 40 m, a altura H., do receptor, de 


1.28 me a distância de Haversive DIST, 


maversnr deixando o teorema da maneira representada 


abaixo na equação (2). 


(Hry — Hay)” + DISTuaversin” = DISTreaL” (2) 


Foram coletados dados de 234 pontos. Os valores da distância real e a potência 
recebida são apresentados na Fig. 4. 
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Fig.4. Gráfico da Distância e Potência Recebida, de todos os pontos coletados. 


Pode-se observar o Recovery Effect [8], verificado em trajetos mistos, a partir da 
distância de 1.6 km até 1.7 km porque é observado um reforço do sinal, isto é, quando o sinal 
passa da estrada para a ponte sobre o rio. Em seguida, em 1.8 km, o sinal decresce devido 
a proximidade de um arco de metal sobre a ponte, mostrado na Fig. 5. A proximidade com o 
arco claramente causa somatório em contra fase do sinal. Por outro lado, quando o receptor 
está posicionado de 1.9 km até 2 km, o sinal apresenta um acréscimo, justamente na frente 
do arco de metal. 





GoogleséFth 


Fig. 5. Arco de ferro presente na ponte. 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 11 





d. Redes Neuro Fuzzy 


Um Sistema Neuro-Fuzzy (SNF) é um sistema híbrido e constituído pela junção de 
duas técnicas de modelagem: as Redes Neurais (RN) e a Lógica Fuzzy (LF). Essas técnicas 
combinadas tornam-se de grande interesse para aplicações principalmente na área de 
engenharia [7]. 

Neste trabalho foi usada a SNF é do tipo Sugeno, o método de treinamento híbrido, 
que combina Least Squares e Back Propagation com Downward Gradient. Outros dados 
importantes são o valor do raio de influência da base da função gaussiana que foi 1, o fator 
multiplicador 1,5, o raio de aceitação no agrupamento 0,75 e raio de rejeição 0,25. 

Para o desenvolvimento do modelo de propagação utilizou-se um total de 234 
amostras. Estas amostras são a potência recebida coletadas na campanha de medição, que 
posteriormente foram transformadas a perda de percurso. As amostras foram misturados 
aletoriamente e divididos em 3 grupos: 70 % para treino, 15 % para validação e 15 % para 
teste. Para a criação do modelo baseado em uma rede Neuro Fuzzy foi utilizado o software 
Mailab. 

As entradas para a rede foram três: distância entre transmissor e receptor, perda de 
espaço livre [16] e a distância somente sobre a água. A saída da rede é a perda de percurso. 
Estes parâmetros são detalhados a seguir: 


Entradas: 


- Distância entre transmissor e receptor: Esta entrada permite que a rede pos- 
sa diferenciar que é o cenário é um trajeto misto, devido que considera a distân- 
cia total e posteriormente somente a distância sobre a água. Para calcular a 
distância entre transmissor e receptor, foi usada a equação (1). 


* | Perda por espaço livre: A maioria de modelos de propagação tem como base 
a perda pela distância sem considerar obstáculos. Isto porque as atenuaçãoes 
causadas pelo ambiente são adicionadas posteriormente, seja como fatores de 
correcão ou variaveis de atenuação, que dependem do tipo de modelagaem. 
Neste estudo a rede calculará os parâmatros adicionais de perda baseada nos 
dados medidos usados para o treinamento. Esta metologia poderia ser compa- 
rada a modelos como Close in [18] ou Floating Intercept [19]. A perda de espaçõ 
livre é calculada como a seguinte fórmula. 


Lp = 32.44 + 20 logd(km) + 20 logf(MHz) (3) 


onde: 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 11 EST 


L,: Perda de percurso no espaço livre; 

d: Distância entre transmissor e receptor em 
km; 

F: Frequência em MHz. 


* Distância somente sobre o rio: No caso de modelos de propagação de traje- 
tos mistos, estes consideram a distância somente sobre a água, como se veri- 
fica nos modelos clássicos como Okumura [20] e ITU-1546-5[21], assim como 
também na modelagem com técnicas de aprendizagem de maquina [22]. As 
distancias usadas na equação (4) podem ser observadas na Fig. 2. 


D, = Dnr — Dr (4) 
onde: 
D,.: Distância sobre o rio ou distância somente sobre a água; 
Dat: Distância horizontal total. 
De: Distância sobre a terra. 


Saída: 


* | Perda de percurso: Usando as entradas antes mencionadas e o treinamento 
da rede, é possível obter a predição da perda de percurso para um trajeto misto 
que considera cidade e rio. Para calcular a saída foram feitas várias simulações 
com diferentes números de iterações. A rede que obteve melhores resultados 
precisou de 1000 iterações, isto porque o parâmetro de avaliação para selecio- 
nar a rede neuro-fuzzy, foi o erro médio quadrático (ERMSE), sendo seleciona- 
da a rede que obteve o menor valor de ERMSE. 


31 RESULTADOS 


A precisão do modelo proposto é avaliada na comparação dos dados medidos no 
entorno da ponte e os dados preditos pela rede Neuro Fuzzy. É importante mencionar que 
os dados de teste são totalmente desconhecidos pela rede. Estes dados estão ilustrados 
na Fig.6, sendo que os vermelhos são dados da rede e os pontos da cor preta é o modelo 
proposto para a rede. Pode-se observar que há coerência entre os dados medidos e os dados 
preditos, mostrando que a técnica de Redes Neuro Fuzzy é satisfatória para a elaboração de 


modelos de propagação, no entanto, está limitado para o cenário no qual foi criado. 
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—— Dados Medidos 
Modelo Proposto 
Okumura Suburbano 
Okumura Misto 


Perda de percusso [dB] 





1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 
Distância [km] 


Figura 6 - Gráfico dos dados medidos, modelo proposto, modelo Okumura-Hata Suburbano e modelo 
Okumura-Hata para trajeto misto. 


O modelo Okumura-Hata é originalmente usado em ambientes urbanos, mas pode 
ser adaptado para cenários suburbanos. As equações (4) e (6) mostram os cálculos de perda 
de percurso utilizados para o modelo Okumura-Hata para ambientes urbanos e suburbanos, 


respectivamente. 
Lurb = 69.55 + 26.16 - log(f) — 13.82 - log(Hyx) 
—(1.1 -log(f) — 0.7) : Hay — 1.56 -log(f) — 0.8 (5) 
+(44.9 — 6.55 - log(Hyy)) : log(d)? 
onde: 
b=1parad<20km (6) 
f 2 
Lisitiio — Las — 2: log (55) + 5.4 + PCO (7) 
onde: 
PCO: Correção da ERP, porque o modelo é feito para 1 Kw de ERP (8) 


O erro RMS obtido é de 2 dB entre os dados medidos e os dados preditos pelo modelo 
proposto, o que indica que o modelo apresenta um bom desempenho para uma ponte sobre 
o rio. Outros trabalhos apresentaram erros de 3 a 3.7 dB para trajetos sobre pontes como 
detalha em [11], embora esta ponte não foi sobre o rio, porém permite reforçar a ideia de que 
o modelo de propagação proposto é válido e fiável para o cenário no qual foi desenvolvido. 

Adicionalmente, na Fig. 6 é mostrado que a partir da distância de 1.5 km até 1.7 km, 
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os dados medidos mostram o Recovery Effect, que se produz ao passar da estrada para a 
ponte. Tal efeito é acompanhado pelo modelo proposto. Uma segunda característica do sinal 
que é acompanhada pelo modelo é diminuição da perda em 1.9 km na frente do arco de metal 
sobre a ponte que se observa na Fig. 4. A partir da distância de 2 km o sinal diminui com a 
distância. 

Na fig. 6 são apresentados os resultados entre o modelo de Okumura misto, suburbano 
e o modelo proposto em cenário de cidade e rio. O modelo proposto se apresenta como uma 
importante alternativa para ambientes de cidade e rio, pois tem a menor perda de percurso 
se comparado aos demais modelos presentes na figura conforme é mostrado na Tabela 1. 


Modelo Erro RMS 
Okumura-Hata Suburbano 16.95 
Okumura-Hata Misto 13.14 
Modelo Proposto 1.98 


TABELA | — Erro RMS dos Modelos 


41 CONCLUSÃO 


Para avaliação dos resultados obtidos no método de predição de perda de propagação 
sobre a ponte de Marabá no Rio Tocantins, a métrica utilizada para avaliar o modelo de 
propagação foi o erro médio quadrático. Dessa forma, o modelo proposto conseguiu predefinir 
a perda de percurso com um erro de RMSE 1.98 dB. 

Os resultados da predição de perda de percurso no trajeto misto estrada sobre terra 
e ponte sobre o rio, foram obtidos através do processamento dos dados no sistema Neuro- 
Fuzzy. Com isso, os resultados foram consistentes para um sistema de comunicação móvel 
no cenário tipicamente amazônico, pois conseguiram acompanhar os dados medidos com 
menor perda de percurso e dessa forma, torna-se essencial para atender de forma satisfatório 
o sistema (LongTerm Evolution) LTE em ambientes não homogêneos. 

Adicionalmente, como forma de otimização do sistema de comunicação, faz- 
se necessário novas propostas de medição para o aprimoramento de estudos de perda 
de propagação para adequação de novas tecnologias, como a LTE-Advanced e o 5G. 
Adicionalmente, para trabalhos futuros serão consideradas as características do rio Tocantins, 


visto que o mesmo tem uma variação de nível do rio de até 10 m. 
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RESUMO - O Presente trabalho descreve 
um conjunto de testes de compatibilidade 
eletromagnética, realizados para avaliar a 
imunidade a surtos elétricos de alta tensão 
em dois tipos de eletrodomésticos: televisor 
e refrigerador. Os testes foram realizados em 
condições específicas definidas pelo padrão 
internacional IEC 61000-4-5 de 2015. Os 
resultados mostraram que o televisor suporta 
maiores níveis de tensão e número de surtos 
que o refrigerador, entretanto, considerando 
que os surtos elétricos induzidos por descargas 
atmosféricas podem atingir níveis da ordem 
de 8kV o uso de dispositivos de proteção 
contra surtos (DPS) é indicado. Os resultados 
mostram ainda que a utilização do DPS pode 
ser imprescindível para manter a vida útil do 
equipamento em regiões com grande incidência 
de descargas atmosféricas. 
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PALAVRAS-CHAVE - Surtos elétricos, proteção contra surtos, refrigerador, televisor e vida 
útil. 


ELECTRICAL SURGE IMMUNITY TESTS IN HOME APPLIANCES 


ABSTRACT - This paper describes a set of electromagnetic compatibility tests, carried out to 
evaluate immunity to high voltage electrical surges, in two types of home appliances: television 
and refrigerator. Tests were performed under specific conditions defined by the international 
standard IEC 61000-4-5 of 2015. The results showed that the TVs supports higher voltage 
levels and more surges than the refrigerator. However, considering that electrical surges 
induced by lightning discharges can be as high as 8kV, the use of surge protection devices 
(SPDs) is indicated. Results also shows that the use of SPDs can be essential to maximize 
useful life of equipment in regions with high incidence of atmospheric discharges. 
KEYWORDS - Surges, lightning protection, refrigerators, television and useful life. 


11 INTRODUÇÃO 


O desenvolvimento tecnológico propiciou o aumento do número de eletrodomésticos 
nas residências, sendo os mais comuns o fogão, a geladeira e a televisão que estão em 
quase todos os lares brasileiros - 98,8%, 97,6% e 97,1%, respectivamente [1]. Segundo a 
Norma Internacional IEC 62305-2 de 2006 [2], os surtos na rede elétrica ocorrem quando a 
carga elétrica é elevada em algum ponto da rede. O surto pode aumentar a corrente que flui 
na instalação e danificar equipamentos elétricos ligados à rede. Diversos eventos podem 
produzir um surto, sendo a descarga atmosférica um dos principais [3]. 

A utilização de DPS (Dispositivo de Proteção contra Surtos) nas instalações elétricas 
propicia vários benefícios aos equipamentos a eles conectados, não só na infraestrutura 
industrial e comercial, mas também no ambiente residencial, na proteção de eletrodomésticos 
[4]. A proteção oferecida pelos DPSs pode prevenir não só a inutilização definitiva e instantânea 
de um equipamento na ocorrência de surtos de alta intensidade, quanto prolongar sua vida útil 
nos casos de surtos de menor intensidade, porém, com maior frequência de incidência, que 
degradam os componentes gradativamente, causando perda de eficiência do equipamento 
até sua completa inutilização. 

Segundo a NBR-5410 [5] toda linha externa de sinal deve ser provida de proteção 
contra sobretensões transitórias, seja a instalação de telefone, de comunicação de dados, 
de vídeo ou qualquer outro sinal eletrônico. O impacto que a utilização do DPS gera na 
manutenção da vida útil dos equipamentos é uma informação importante para determinar 
a sua eficiência e, sobretudo, justificar o custo da sua aplicação nas instalações elétricas 
dos vários tipos de infraestrutura. Dessa forma, o presente trabalho apresenta testes de 
compatibilidade eletromagnética que mostram como o DPS pode ser empregado na 
manutenção da vida útil de televisores e refrigeradores. 
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21 APARATO EXPERIMENTAL 


Nos sistemas de distribuição de baixa tensão as descargas atmosféricas podem 
causar danos aos dispositivos elétricos, porém, apenas 10% desses danos são causados 
por descargas diretas. Os 90% restantes são provocados por tesões induzidas [6]. Segundo 
Silva [7], uma descarga atmosférica no primário de um transformador de distribuição (13.8kV 
— 220/127 V, 30kVA) gera uma tensão induzida no secundário que pode chegar a 8kV. Além 
disso, raios adjacentes podem induzir, por acoplamento, tensões transitórias na linha, as 
quais serão transferidas para os dispositivos que estiverem conectados a esta linha [8 — 9]. 
O DPS é utilizado para limitar a tensão transitória que atingiria os equipamentos a um nível 
mais baixo, atuando em um tempo muito curto e protegendo os equipamentos conectados 
ao circuito [10]. 

No Brasil os equipamentos eletrodomésticos são normatizados pela NBR 5410 [5]. 
Essa norma estabelece uma tensão impulsiva mínima de 1,5 kV que os eletrodomésticos 
devem ser capazes de suportar sem apresentar falhas. Porém, segundo o Institute of 
Electrical and Electronic Engineers — Surge Protective Committee (IEEE - SPDC), a tensão 
na entrada de uma edificação fica limitada em 6 kV [11], que é superior à proteção exigida 
pela NBR 5410 e pode levar a uma falha nos eletrodomésticos. Nesse contexto, o presente 
trabalho descreve a realização de testes de compatibilidade eletromagnética regidos pela 
Norma IEC 61000-4-5 de 2015 [12] em um modelo de televisor e um de refrigerador. Segundo 
esta norma, os testes devem ser realizados sob a configuração apresentada na Figura 1. 

O gerador de onda combinada produz os surtos elétricos padronizados que são 
injetados nos conectores Neutro e Fase de uma rede elétrica de testes com o auxílio da 
rede de acoplamento, conforme ilustrado na Figurai. O gerador utiliza um aterramento 
exclusivo, através do qual ele descarrega o surto que retorna do circuito de teste. A rede de 
desacoplamento, ilustrada na Figura 1, é utilizada para evitar que o surto injetado passe para 
as instalações elétricas locais e provoque danos em outros equipamentos que não fazem 
parte do teste. 

Uma tomada de 10 A é utilizada na saída da rede de acoplamento para que o 
equipamento sob teste esteja em condições reais de uso. A seguir descreve-se o gerador de 
onda combinada, a rede de acoplamento e desacoplamento utilizados nos experimentos bem 
como uma estimativa do impacto do uso do DPS na vida útil dos equipamentos. 
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Aterramento | Combinada 
Exclusivo | 


Rede de 
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Porta para 
Equipamento 
sob Teste 


Fig 1 - Configuração das conexões do equipamento sobre teste para os testes de compatibilidade. 


Fonte: Elaborada pelo autor. 
a. Gerador de Onda Combinada 


O gerador de onda combinada produz surtos elétricos com forma de onda 
padronizada pela IEC 61000-4-5 [12] e parâmetros específicos como tensão de circuito aberto 
e corrente de curto circuito. A Figura 2 apresenta a forma de onda da tensão normalizada, 
na saída do gerador, quando em circuito aberto, atendendo às exigências da IEC 61000- 
4-5. Essa tensão de saída foi aplicada às cargas durante os experimentos. Dois tipos de 
geradores de ondas combinadas são especificados em função do tempo de subida (T.) e 
tempo de duração (T,), ilustrados na Figura 2. Para a onda de tensão esses geradores são 
definidos por: 10/700 us (T;, = 10 use T, = 700 us) e 1,2/50 us. Cada um tem suas aplicações 
específicas dependendo do tipo de porta a ser testada. 

O gerador de ondas combinadas de 10/700 us é usado para testar portas de 
comunicação simétricas que se destinam a ser diretamente conectadas a redes de 
telecomunicações externas. O gerador de ondas combinadas de 1,2 / 50 us é usado em 
todos os outros casos e é exatamente o gerador aplicado para testar os eletrodomésticos 
neste trabalho. Esse gerador tem uma resistência interna de 2 O e os parâmetros das 
formas de onda na saída (de tensão em circuito aberto exibida na Figura 2 e da corrente 
em curto-circuito, que pode ser encontrada na norma) estão apresentados na Tabela 1. 
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Normalizada 














Fig. 2. Forma de onda da tensão normalizada aplicada aos equipamentos. 


Fonte: Adaptada da Norma IEC 61000-4-5 [12]. 


Tempo de Subida Tempo de Duração 
(TS) Us (TD) us 


Tensão em Circuito Aberto TS =1,2+30% TD = 50 + 20% 
Corrente de Curto Circuito TS=8 +20% TD = 20 + 20% 


TABELA |. Parâmetros da forma de Onda do Gerador 1,2/50 [12]. 





b. Rede de Acoplamento 


A rede de acoplamento é responsável por permitir a inserção do surto de tensão 
proveniente do gerador na rede elétrica do equipamento sob teste, conforme ilustrado 
na Figura 1, e impedir que a tensão da rede local entre no gerador. O circuito da rede de 
acoplamento é um varistor, que permite a passagem do surto de alta tensão e impede a 
passagem da tensão da rede local para o gerador, em série com um fusível térmico, que 
em caso de falha do varistor desconecta o gerador da rede de testes. A saída da caixa 
de acoplamento é uma tomada de 10A padrão, na qual os equipamentos sob teste são 
conectados e recebem a alimentação para operação junto com o surto elétrico na tensão de 
testes definida no gerador. 


c. Rede de Desacoplamento 


A rede de desacoplamento permite a conexão do equipamento sob teste à rede 
elétrica local, que no laboratório é de 220V,.. - 60Hz, e impede a passagem do surto 
injetado no circuito pelo gerador para a rede local. A rede de desacoplamento é formada 
por duas bobinas de indutância 1,5 mH, como definido pela norma IEC 61000-4-5 [12]. As 
bobinas foram colocadas dentro de uma caixa de proteção conforme ilustrado na Figura 3. 
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Um transformador de corrente foi instalado no condutor fase para medição da corrente e 
confirmação de que o surto elétrico injetado pelo gerador não atinge a rede elétrica local. 


31 RESULTADOS E CONSIDERAÇÕES 


O aparato experimental apresentado na Figura 1 foi montado no Instituto de Inovação 
Tecnológica da Universidade de Pernambuco IIT/UPE e um televisor e uma geladeira foram 
submetidos a surtos elétricos de alta tensão para avaliar a necessidade de utilização do 
DPS para proteção desses equipamentos. Também foram utilizados nos experimentos DPSs 
Classe Ill [13], aplicados em série à entrada de potência do equipamento sob teste. 

Os ensaios foram realizados sob duas hipóteses de resultados: quantidade de surtos/ 
pulsos até ocorrer a falha do equipamento ou não falhar. Ocorrendo a falha, a quantidade total 
de pulsos aplicados até então é registrada e, quando o aparelho não falhou, foi registrado 
700 que é o número máximo de surtos ao qual cada equipamento poderia ser submetido. O 
tempo entre surtos aplicados pelo gerador é de 30 segundos e, para que o teste em um único 
equipamento fosse realizado em um dia, definiu-se que o número máximo de surtos seria de 
700, que corresponde a aproximadamente 6 horas. As marcas e modelos dos equipamentos 
não serão mencionados, mas são equipamentos amplamente comercializados no Brasil. 







Saída para Rede 
de Acoplamento 


| 
| 


Bobina da Fase 


Transformador 
de corrente 


Fig. 3. Fotografia do circuito da rede de desacoplamento. 


Fonte: Fotografia do autor. 


a. Experimentos com Televisor 


Um televisor com fonte de alimentação externa foi adquirido e submetido a surtos 
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elétricos de até 6kV. A Tabela Il apresenta os resultados dos testes com esse equipamento. 
Observa-se que o televisor não apresentou falhas, mesmo sem a utilização dos DPSs Classe 
Il na sua proteção, para surtos de até 6kV. Esse resultado mostra que para surtos até esse 


nível de tensão o televisor escolhido não precisa da proteção de um DPS. 


Televisor sem DPS 





S = Não ocorreu Falha no Equipamento. 


TABELA Il. Dados Experimentais com o Televisor 


Entretanto, os surtos induzidos no secundário de um transformador podem chegar 
a 8kV [7], o que não foi testado pois o gerador disponível não atingia esse nível de tensão. 

Os resultados motivaram a abertura da fonte de alimentação do televisor para 
constatar se o equipamento possui proteção interna específica. A Figura 4 apresenta uma 
fotografia da fonte de alimentação do televisor, na qual destaca-se o varistor responsável 
pela proteção desse dispositivo. Assim, ficou constatado que o circuito interno de proteção 
do modelo de televisor escolhido possui proteção eficiente contra os níveis e quantidade 
de surtos aplicados neste trabalho. Vale salientar que apesar do modelo possuir proteção 
interna o uso do DPS seria recomendado, uma vez que evitaria que o televisor recebesse os 
surtos diretamente e garantiria um tempo de vida ainda maior para o equipamento. 


a os 


— 
o 
= 





Varistor 





q 





Fig. 4. Fotografia da fonte de alimentação do televisor utilizado nos testes. 


Fonte: Fotografia do autor. 
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b. Experimentos com Refrigerador 


Foram realizados dezoito (18) experimentos com um refrigerador comercial sendo 
nove com o refrigerador conectado diretamente ao circuito gerador de surtos e nove 
com o equipamento protegido por um DPS comercializado no Brasil. Por se tratar de um 
equipamento de custo elevado os experimentos foram realizados com um único refrigerador 
e as falhas apresentadas eram corrigidas pela equipe da UPE antes da próxima sequência 
de testes. A falha mais comumente apresentada foi a queima de um capacitor na placa de 
potência, apresentada na Figura 5. Essa placa foi substituída por uma nova em todos os 
testes. Realizava-se uma bateria de testes até a falha do refrigerador e, na sequência, as 
peças que apresentavam falha eram substituídas e o equipamento ficava em observação. 
Após comprovação de que o equipamento funcionava corretamente uma nova bateria de 


testes era realizada. Os resultados são apresentados na Tabela III. 





Refrigerador sem DPS 














Teste 3 kV 4 kV 5 kV 
1º 700 (S) 114 (FE) 1 (FE) 
2º 700 (S) 225 (FE) 1 (FE) 
3º 700 (S) 700 (S) 1 (FE) 














Refrigerador com DPS 














Teste 5 kV 6 kV 7kV 
1º 700 (S) 160 (FD) 74 (FD) 
2º 700 (S) 104 (FD) 55 (FD) 
3º 700 (S) 256 (FD) 53 (FD) 
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FE = Ocorreu Falha no Equipamento / FD = Ocorreu Falha no DPS / S = Não ocorreu Falha no 
Equipamento nem no DPS. 
TABELA III. Dados Experimentais com o Refrigerador. 


Os resultados experimentais apresentados na Tabela Ill mostram a importância do 
uso do DPS no modelo de refrigerador escolhido, uma vez que um único surto de 5kV produz 
falha no equipamento. Com a proteção do DPS o conjunto suportou em média 60 surtos de 
7kV. Vale salientar, que a falha apresentada foi apenas no DPS, ou seja, a placa de potência 
e todas as funcionalidades do refrigerador foram mantidos mesmo após o DPS apresentar 
falha com surtos de 7kV. Esses resultados mostram que o uso do DPS pode garantir que o 
equipamento não apresente falha precoce em decorrência de surtos elétricos de alta tensão. 
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E e . cá 


Fig. 5. Fotografia do componente que falhou na placa de potência do refrigerador após surto elétrico. 
Fonte: Fotografia do autor. 


c. Análise sobre a Vida Útil dos Equipamentos 


O estado brasileiro com a maior densidade de raios em 2019 foi Tocantins, com 19,8 
raios por km? por ano [14]. Com os experimentos realizados foram determinados o número de 
surtos que cada equipamento pode suportar. Dessa forma, determinando o número de surtos 
que pode atingir uma residência, obtém-se o tempo de vida do equipamento considerando 
que a única falha possível seja devido a surtos provocados por descargas atmosféricas. 

Um surto elétrico pode atingir um equipamento em uso na residência quando uma 
descarga atmosférica incide: a) na estrutura ou em uma área próxima à estrutura num raio 
de 500 m [15]; b) na linha externa de alimentação ou em uma área próxima a ela. A área 
considerada como próxima à linha consiste num comprimento de 1km de linha e distante 
até 2 km da linha (para um lado ou para o outro) [15], conforme ilustrado na Figura 6. 
Desconsiderando a sobreposição entre as duas áreas, a área total de influência de uma 
descarga na edificação é de 4,39 km?. Assim, no estado de Tocantins, uma residência 
estaria sujeita a aproximadamente 87 raios por ano. Diante do exposto, os resultados com 
o refrigerador mostram que o equipamento poderia apresentar falha logo no início de uso 
sem a proteção do DPS uma vez que ele falhou com apenas um surto de 5kV. Com o DPS 
o conjunto suportou os 700 surtos (podendo suportar mais já que 700 foi o limite de pulsos 
aplicados no teste), considerando 85 surtos por ano isso corresponde a aproximadamente 8 


anos. 
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Residência 


[| Área Próxima à Estrutura (Raio de 500m) 





[|] Área Próxima à Linha Externa de Alimentação (1km x 4km) 


Fig 6 - Área de influência de uma descarga atmosférica na edificação. 


Fonte: Elaborada pelo autor. 


41 CONCLUSÕES 


O presente trabalho mostra como os testes de compatibilidade eletromagnética, 
em especial os de imunidade contra surtos elétricos de alta tensão possuem relevância 
na determinação da qualidade dos produtos produzidos no país e na comprovação da 
necessidade do uso de DPSs para proteção de equipamentos instalados nas regiões 
de incidência de descargas atmosféricas. Os resultados experimentais comprovam a 
necessidade de utilização do DPS nessas regiões e mostram que seu uso pode evitar a 
falha em um equipamento como um refrigerador logo no início de sua vida útil. Vale salientar 
que nos testes realizados nenhum equipamento apresentou falha quando utilizado em 
conjunto com o DPS. A análise sobre a vida útil do equipamento, apresentada neste trabalho, 
mostra um resultado importante no tocante ao cálculo do número de surtos que podem ser 
produzidos pela ocorrência de descargas atmosféricas. Esse resultado é relevante para o 
aprimoramento de circuitos de proteção das edificações, dos eletrodomésticos e dos DPSs. 
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RESUMO: Alguns sistemas de 
existentes nas edificações são arcaicos 
e consomem muita eletricidade quando 
comparados às novas tecnologias. Este artigo 
relata ensaios realizados com o sistema Dali - 
Digital Addressable Lighting Interface e lâmpadas 
fluorescentes tubulares de alta eficiência para 
determinar a economia de energia gerada em 
relação à iluminação tradicional e avaliar a 
viabilidade econômica da solução. 
PALAVRAS-CHAVE: Eficiência energética, 
controle da iluminação, sistema DALI, automação 
de iluminação 


iluminação 


LIGHTING AUTOMATION AND ENERGY 
EFFICIENCY IN BUILDINGS - THE DALI 
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CASO 


LIGHTING CONTROL SYSTEM: A CASE 
STUDY 


ABSTRACT: Some existing lighting systems 
in buildings are archaic and consume a lot of 
electricity when compared to new technologies. 
This article reports tests carried out with the Dali 
- Digital Addressable Lighting Interface system 
and tubular fluorescent lamps to determine the 
energy savings generated in relation to traditional 
lighting and to assess the economic viability of 
the solution. 

KEYWORDS: Energy efficiency, lighting control, 
DALI system, lighting automation. 


11 INTRODUÇÃO 


Com o advento da revolução industrial e 
a produção de grandes volumes de bens através 
da construção de fábricas, onde os meios de 
produção eram aglutinados com as pessoas 
que iriam operá-los para depois embalá-los e 
distribuí-los, não houve uma preocupação com 
a degradação que ocorria no meio ambiente. 
Acreditava-se que a disponibilidade de recursos 
era infinita e que os rejeitos dos processos 
poderiam ser descartados sem problemas. 

As alternativas que surgiram para 
solucionar essa questão foram evoluindo ao 
passo em que o consumo e o crescimento 
urbano foram aumentando de forma exponencial. 
Processos de descarte e reciclagem, uso racional 
de materiais, novos produtos sintéticos que 
substituem matérias-primas de recursos naturais, 
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entre outros avanços, marcaram o fim do século XX e a entrada deste novo século. Entretanto, 
o foco das grandes organizações e potências econômicas está voltado para um problema: 
“como continuar atendendo à demanda por energias, sobretudo elétrica, considerando o 
acelerado crescimento populacional frente às limitações da produção energética”. 

Embasado por essa premissa, este trabalho conduz a uma análise sobre as vantagens 
de se usar um sistema de automatização eletrônico de fontes de energia, especificamente o 
controle de iluminação artificial, no intuito de satisfazer não apenas à demanda, mas melhorar 
a eficiência dos projetos e contribuir com as discussões para eficiência energética visando 
políticas de sustentabilidade para o país. 


2 | EFICIÊNCIA ENERGÉTICA EM SISTEMAS PREDIAIS 


A necessidade da utilização mais eficiente da energia elétrica em edifícios, levou ao 
desenvolvimento técnico em diversos sistemas prediais onde este insumo é utilizado. A figura 
2 indica, no gráfico, o uso final energia elétrica em edifício comercial e permite verificar que os 
maiores sistemas consumidores são: ar condicionado, 42%; iluminação, 24%; equipamentos 
de escritório, 15% e elevadores e bombas, 13%. 


Uso Final da Energia Elétrica em Edifício 
Comercial 


E iluminação 
E ar condicionado 
E equipamento escritório 


E elevadorese bombas 





Figura 2 - Uso final da energia elétrica em edifício comercial. 


Fonte: (LAMBERTS, 2013). 


2.1 Sistemas de Controle de Iluminação 


Os sistemas de controle de iluminação passaram de sistemas isolados, analógicos 
para digitais centralizados, conforme a figura 3. 
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Controle Analógico D 54 (padrão Europeu) 


AMX 192 (padrão americano) 


Padrão DSI 


Padrão DALI 


Controle Digital 


Padrão KNX 


Padrão DMX 512 






Wireless 


Bluetooth 


Figura 3 - Sistemas de Controle de Iluminação: Analógico, Digital e sem fio 


Fonte: Própria 


2.2 DALI — DIGITAL ADDRESSABLE LIGHTING INTERFACE 


2.2.1 Histórico 


Entre todos os protocolos digitais e analógico, o mais flexível e promissor para 
controles de iluminação e novas aplicações futuras é o DALI. 

O padrão DALI, que é um acrônimo de Digital Addressable Lighting Interface, é um 
protocolo aberto que surgiu no final do ano de 1990, com a iniciativa dos principais fabricantes 
e desenvolvedores de sistemas de controle digital e analógico de iluminação e luminárias da 
época, como OSRAM, Helvar, Trilux, Tridonic e outros líderes. Ele é especificado pela norma 
internacional IEC 62386, sendo irrestrito a qualquer fabricante. 

Até junho de 1996, apenas 34 marcas participavam do protocolo DALI e hoje já 
são mais de 120 empresas, organizações e desenvolvedores que participam do grupo de 
membros DALI. É possível conferir todos os membros, além de projetos já realizados com o 
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protocolo, no endereço eletrônico www.dali-ag.org. 
No início, intuito era desenvolver um padrão de protocolo digital intercambiável entre 
diversos fabricantes, que facilitasse a aplicação e permitisse: 


- Baixo custo de instalação, com redução dos cabos e fiação dos controles e 
comandos. 


- Controle individual de lâmpadas e luminárias dentro de circuitos complexos. 


“Capacidade de adicionar componentes, sensores e outros equipamentos 
propriedade. 


-Via de informação e controle bidirecional entre a fonte e o controlador 


Desde sua criação, foi submetido a vários testes que lhe possibilitaram muitos 
avanços e desde o 2001 assumiu importante papel na automatização de edifícios e 
indústrias, contribuindo diretamente para a eficiência energética das instalações. Um 
exemplo bastante divulgado do uso deste protocolo é o caso da Copa do Mundo de 2006, 
que teve o Estádio Olímpico de Berlim totalmente iluminado e controlado por centrais DALI. 

Para a indústria o DALI trouxe muitos benefícios e também mudou a maneira como 
projetistas, light designers, instaladores e o próprio conceito das aplicações luminotécnicas, 
que agora contam com mais flexibilidade para aplicações simples ou complexas. 


2.2.2 Áreas de aplicação do DALI 


Ao contrário dos protocolos digitais que surgiram logo em seguida dos analógicos, 
o DALI não é uma releitura simplificada, ele vai além disso. Em poucos anos de avanço 
tecnológicos as empresas-membro do DALI desenvolveram e aperfeiçoaram as ferramentas 
de análise e controle. 

Uma notável aplicação é marcada pela variação de tipos de luminárias, reatores e 
tipos de lâmpadas que hoje podem ser controladas por uma interface DALI. Inicialmente, o 
sistema era focado para luminárias de lâmpadas fluorescentes ligadas a reatores. Atualmente 
são agregados de maneira simples lâmpadas halógenas de 12 volts, lâmpadas de descarga 
e principalmente LED, não em sua forma pura, mas através de seus drivers de acionamento 
e controle. 

Em vez de controlar apenas o sistema de acendimento e tempo de funcionamento das 
lâmpadas, o DALI oferece um mecanismo simples de controlar a luminosidade, por meio da 
integração de sensores, retornando ao controlador uma mapeamento completo e instantâneo 
sobre a performance do sistema. Antes do seu desenvolvimento, esse tipo de automatização 
era muito oneroso e pendia de muita infraestrutura, além se ser restrito a marcas que não 
compartilhavam seus sistemas para intercambiar com outros fabricantes. 

Outra vantagem do DALI é a descentralização dos controles, permitindo mais 
flexibilidade ao processo como um todo. Desta forma, grandes edifícios e empresas podem 


ter polos de controle e operação distribuídos por toda a instalação. Isso não impede que o 
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projeto do DALI seja unificado a qualquer momento para uma única central de controle. 

Em termos de proporção e adequação do uso aos ambientes pode-se dizer que o 
DALI é altamente adaptável. Sua aplicação pode partir de uma única luminária com um único 
reator acoplados a um multisensor, que irá oferecer controle de luminosidade, detecção 
de movimento por sistema de infravermelho passivo (PIR — passive infra-red sensor) e até 
operação por controle remoto do tipo IR (infra-red). Esta configuração mínima requer um 
simples painel de controle ou o próprio controle remoto e se classifica como uma instalação 
autônoma de controle. 

Partindo para uma configuração mais ampla, o DALI oferece a opção de se combinar 
com outros sistemas de controle, como um subsistema e como um sistema agrupado, na 
qual a central de controle energético do edifício poderá monitorar e designar operações ao 
sistema, corrigindo demandas, acionamentos e até mesmo cenários de iluminação. 

A parte estrutural também ganha com o sistema DALI. Por se tratar de um sistema 
bidirecional, onde apenas um par de cabos envia e recebe informações protocoladas, a 
qualquer momento uma nova luminária pode ser inserida ou retirada do trecho sem ter que 
passar reformas ou alterações complexas. 

Em termos técnicos o DALI também suas limitações. Um dispositivo controlador pode 
coordenar até 64 dispositivos em suas saídas. Essa linha de controle dispensa o uso de relés 
e pode ser associada a vários sensores espalhados pelo percurso, que é composto por uma 
estrutura cabeada muito simples. As figuras 4 e 5 exemplificam a instalação básica e suas 


limitações. 


Neutral conduct 





DALI (DA) 


Protective 
earth 


DALI (DA) 


Phase 


Figura 4 - Esquema de cabos estruturados do DALI 


Fonte: OSRAM 
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The DALI principle 


— :"É pa 
= mu control unit 

Power supply ,3 

DALI 7 max. 

control cable 300 m 


Alternative group assignment 
without changing the wiring 





Maximum of 64 units can be connected 


Figura 5 — Princípio básico de ligação ao controlador 


Fonte: OSRAM 


Basicamente, um controlador pode ser comparado a um controlador lógico programável 
— CLP — com a diferença que o DALI é exclusivo para iluminação e seus mecanismos. 
Entretanto à sua extensiva funcionalidade e flexibilidade, projetos com esse protocolo podem 
ser facilmente integrados a outros sistemas de automação de um edifício via gateways. 

Como em todo sistema de controle o DALI também necessita de um certo treinamento 
na operação e desenvolvimento de toda a lógica que comandará a iluminação. Porém, as 
empresas já passam a oferecer cursos ou o próprio serviço de projeto, facilitando ainda mais 
para quem deseja instalar o DALI. 

Simplificadamente, um sistema DALI, composto por um controlador unitário de 64 
saídas, poderá ser programado via software, normalmente por uma plataforma Windows, e 
depois poderá ser controlado por algum mecanismo, seja por um computador ou software 
específico para algum mecanismo de interface homem-máquina. Entretanto, existem 
mecanismos controladores que acionam diretamente o controlador DALI, como por exemplo 
interruptores com dimerização digital. Alguns fabricantes também dão ao DALI um ar mais 
tecnológico, com sensores e acionadores que funcionam em paralelo e interoperáveis, 
permitindo ampliação das áreas de controle individual ou coletivas. 

Com esse conceito simplificado de instalação do DALI, muitos edifícios comerciais, 
indústrias e até mesmo residências ou pequenos comércios investiram na aplicação dessa 


tecnologia e passam a se enquadrar como “edifícios verdes”. 
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31 METODOLOGIA 


O sistema de controle Dali modelo ECO compreende um controlador, um sensor 
de presença e luminosidade e reatores. No projeto aqui descrito, ele foi instalado em uma 
residência, numa área com 2,75m de largura, 4,70 m de comprimento e 2,20 m de altura. 
Por meio do software Lumisoft, da Lumicenter, baseado no método dos lúmens. Foram 
especificadas duas luminárias de 2 x 40W para obter nível de iluminamento de 500 lux. Eles 
foram instalados no ambiente e a coleta de dados foi feita no monitor de uso de eletricidade 
P4400, da P3 International Corporation. As medições de luminosidade foram feitas com o 
Digital Lux Meter MLM-1011, da Minipa. 

O controlador “dimeriza” as luminárias do laço DALI-1 e no laço-2 dá um offset em 
relação à máxima capacidade de iluminação do sistema. 

O acendimento das lâmpadas ocorre quando uma pessoa entra no ambiente, 
conforme indica a figura 6, no primeiro símbolo no eixo das ordenadas. O desligamento 
ocorre 15 minutos após a saída do último ocupante da sala. 


EJINE: 





Figura 6 - Desenho esquemático do acendimento, “dimerização” e desligamento das lâmpadas do 
controlador DaliECO 


Fonte: Manual DaliECO application guide v01b 
Foram feitas medições de consumo com lâmpadas de 40 W e 36W com reator 


eletrônico tradicional sem dimerização e de 28W com reator Dali sem dimerização. Os 
valores obtidos são apresentados na tabela 1. 
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P Ss F Economia Economia 
Tipo de lâmpada V(V) I(A) (W) (VA) (Hz) FP (W)(2) (VA)(2) 


OsmmEDN atas [DM [65:84 RETOS | 60 | (0:66 
luz do dia especial 


Philips 40 W (1), TIO - 450 066 68 784 60 0,86 
branca luz do dia 


Slim Fluor ECP 
6400K, 36 W(1), T8- 116,9 0,51 58,6 60,6 59,9 0,96 -13,82 -22,70 


branca luz do dia 


Slim Fluor ECP 
6400 K,36W(1), T8 115,9 0,52 586 60,1 60 0,96 -13,82 -23,34 


branca luz do dia 


Coramad 19,6 049 565 585 60 096 -16,91  -2538 
T5 - cool white | 


(1) Considerando duas lâmpadas de cada potência (40, 36 e 28 W) 


(2) Economia comparada com a lâmpada Philips 


TABELA 1: Valores de consumo das lâmpadas fluorescentes tubulares de 40, 36 e 28W sem 
dimerização. 


Já a tabela 2 mostra os dados obtidos dos ensaios com a dimerização acionada 
em diversas condições: com as janelas interna e externa fechadas ou abertas e com a 


iluminação do corredor de acesso ligada ou desligada. 


Ensaio 1 Ensaio Ensaio Ensaio Ensaio Ensaio Ensaio Ensaio Ensaio Ensaio Ensaio Ensaio 


(1) 2 3 4 5 6 7 8 g 10 11 12 
sea Fechada Aberta Aberta 50% Fechada Aberta Aberta 50% Aberta Fechada (4) Aberta 
ia Fechada Aberta Fechada Fechada Fechada Aberta Fechada Fechada Aberta Fechada (4) Aberta 
Corredor Apagado Aceso Apagado Apagado Apagado Aceso Apagado Apagado Apagado Apagado (4) Apagado 
V(v) 117 1167 1167 1162 116,1 1164 1169 1169 1182 177 1162 1185 
(A) 0,98 0,82 0.86 0.94 0,98 0,94 0,91 0.94 0,69 0,97 1,04 0,52 
P(W) 111 90,7 96,5 105 111 105 102 106 69.4 109 118 51,7 
S (VA) 114 95,9 101 109 114 109 106 110 76,7 114 121 61,5 
f (Hz) 60 60 60 59,9 59,9 60 60 60 30 60 60 59,9 
FP 0,97 0.95 0.95 0.96 0,97 0,96 0,96 0,93 0,9 0,95 0,97 0.83 
(2) (lux) 830 850 820 905 850 850 750 750 350 890 610 360 
(3) (lux) 350 360 360 230 160 240 230 210 230 160 110 1240 


(1) Medido logo após término do setup 
(2) Lux, altura de trabalho 
(3) Luz externa a um metro da janela externa, sensor a 2,05 m 


(41 À noite. com todas as ianelas fechadas 


Tabela 2: Ensaio do sistema Dali com controlador, uma luminária no canal 1 e uma luminária no canal 2 
(luminária de duas lâmpadas T5/28W). 


Na tabela 3, é feita uma comparação entre as lâmpadas fluorescentes tubulares de 
40W e o sistema com controle dimerizado. Na melhor condição, considerando-se a potência 


aparente, tem-se uma economia de 51%. 
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Condição 
de 
trabalho 


Lâmpada 


Grandeza 
V(M) 

I(A) 
P(W) 

S (VA) 

f (Hz) 

FR 
Economia 
(W) (4) 


Economia 
(VA) (3) 


Sem dimerização 


Tubular Tubular 
Osram 40 W Philips 40 W 
(1), T1O— (1), TÃO — 
luz do dia branca luz 
especial do dia 
118,3 120 
1,24 1,32 
131,6 136 

152 156,8 
60 60 
0,86 0,86 


Sistema Dali com uma luminária 2 x 28 W na 
rede Dali 1 e uma luminária 2 x 28 W na rede 


Set 


up de 
fábrica 


TS, 28 W 


119,6 
0,98 
113 
117 
60 
0,96 


-16,91 


-29,38 


Após 
setup 
inicial 


116,1 
0,98 
111 
114 
60 
0,97 


-18,38 


-27,30 


Dali 2 


Pós 


setup - 


Caso 
médio 


116,7 
0,86 
96,5 
101 
60 
0,95 


-29,04 


-39,99 


Pós 
setup - 
Melhor 

caso 


116,7 
0,82 
90,7 
99,9 
60 
0,95 


-33,31 


-38,84 


Pós 
setup - 
Melhor 

caso 


(2) 


118,2 
0,69 
69,4 
76,7 
60 
0,9 


48,97 


-51,08 


(1) Para mesma base de comparação, foram consideradas quatro lâmpadas T10, 40 W, com quatro 
lâmpadas T5, 28 W. 


(2) Simulado um dia mais claro, pois as simulações anteriores foram realizadas em dias nublados. 


(3) Economia comparada com a lâmpada de 40 W, Philips. 


Tabela 3: Comparação entre o sistema dimerizado e uma lâmpada fluorescente tubular de 40W sem 


dimerização. 


Os resultados obtidos foram os esperados, visto que as lâmpadas T5, 28W têm 
eficiência luminosa superior às de 40W. Além disso, com a dimerização obteve-se uma 


economia extra, pois o consumo diminui quando existe iluminação natural aproveitável. 


Adicionalmente o sensor de presença desliga a iluminação quando não há pessoas na sala. 


Na tabela 4 tem-se a especificação das lâmpadas. 


Potência (W) 


Temperatura de cor (K) 
Fluxo luminoso (Im) 


Philips 40 W, 
Extra luz do dia 
TL 40 W/ 750 


40 
5000 
2700 


Eficiência luminosa 


(ImAN) 


Vida útil (horas) 


Diâmetro 


70 


20 000 
Tio 


Osram 40 W 
Luz do dia 
especial 


40 


Tio 


ECP 36 W  Osram 28 W 
Branca luz Cool white HE 


Tabela 4: Especificação das lâmpadas utilizadas. 
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do dia 28W/840L 
36 28 
6400 
2900 
103 
20 000 
T8 o 
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41 ANÁLISE FINANCEIRA: 

Para verificação da viabilidade financeira da solução, foi feito um estudo considerando 
a substituição de 484 lâmpadas no denominado bloco D do campus São Paulo da instituição 
de ensino IFSP, Campus São Paulo. 

A tabela 5 apresenta o consumo anual de energia, em kWh e R$, das lâmpadas de 
40W e das de 28W dimerizadas. Foi considerada a tarifa verde da Enel, de março de 2016, 
sem bandeira tarifária ou impostos. 

A tabela 5 mostra uma estimativa de custos para a substituição das 484 lâmpadas de 


40W por lâmpadas de 28W dimerizadas por um controlador Dali. 


Estudo considerando 484 lâmpadas 


aa horo-sazonal verde (R$/kWh ponta e R$/kWh fora de ponta), AES Eletropaulo, março 03624 022596 
484 lâmpadas 40 W (W) 16 456 

484 lâmpadas 28 W, com sistema Dali de controle (W), considerando economia de 45% 9051 

Utilização de horas fora de ponta em um dia 15 

Utilização de horas na ponta em um dia 3 

Utilização de horas fora de ponta (sábado) 10 

No mês 

Utilização de horas fora de ponta (considerando 22 dias) 330 

Utilização de horas na ponta (considerando 22 dias) 66 

Utilização de horas fora de ponta (quatro sábados) 40 

Total de horas fora de ponta 370 

Total de horas na ponta 66 

Comparação entre lâmpadas de 40 e 28 W: energia consumida e valor 

Consumo de energia fora de ponta (kWh) 6089 3349 
Consumo de energia na ponta (kWh) 1086 597 
Valor do consumo de energia fora de ponta, sem impostos (R$) 1375,81 756,69 
Valor do consumo de energia na ponta, sem impostos (R$) 393,60 216,48 
Valor mensal do consumo de energia (R$) 1769,41 973,17 
Total anual do consumo de energia elétrica (R$) 21232,90 11678,10 
Economia (R$/ano) 9554,81 


Tabela 5: Consumo de energia (em kWh e R$) de lâmpadas de 40W e de 28W dimerizadas. 


5| CONCLUSÕES 
De acordo com os ensaios realizados, o sistema de dimerização fornece economia de 


40%, na melhor condição. 
No caso de um edifício em que o sistema de iluminação representa 50% do total de 
consumo de energia elétrica, na melhor condição, tem-se economia de 20%. 


No caso de um retrofitting, o payback é de 6,3 anos. 
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RESUMO: A geração de energia eólica vem 
aumentando nas últimas décadas e com isso 
também aumentaram os estudos voltados 
para a caracterização da sua eficiência — que 
depende de seus parâmetros construtivos, assim 
como das variáveis ambientais envolvidas no 
processo. O objetivo desta pesquisa foi identificar 
a influência das variáveis no eixo de rotação de 
aerogeradores. A partir de valores fornecidos 
pela literatura e pelos fabricantes de turbinas 
eólicas foram combinadas as seguintes variáveis: 
diâmetro do rotor; velocidade do vento; e ângulo 
de passo da turbina para verificar sua relação 
com a rotação e potência disponível no cubo de 
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rotores eólicos. Tal combinação foi construída 
a partir da modelagem das equações que 
representam a geração eólica com os cenários 
sendo construídos sob delineamento composto 
central rotacional (DCCR). Os resultados 
apresentam dependência combinadas das 
grandezas construtivas e ambientais com grau 
de dependência apresentado por duas equações 
teóricas. A conclusão desta pesquisa mostra a 
necessidade de ampliar os estudos para validar 
a dependência matemática entre a rotação do 
eixo e o ângulo de passo a fim de aumentar a 
acurácia dos seus sistemas de controle. 
PALAVRAS-CHAVE: Velocidade do vento. 
Energia eólica. Rotação. 


INFLUENCE OF THE ENVIROMENTAL 
AND CONSTRUCTIVE VARIABLES ON 
THE AXIS OF WIND ROTOR 


ABSTRACT: The generation of wind energy has 
been increasing in recent decades and with this, 
studies aimed at characterizing its efficiency have 
also increased — which depends on its construction 
parameters, as well as on the environmental 
variables involved in the process. The objective 
of this research was to identify the influence of 
variables on the rotation axis of wind turbines. 
From values provided by literature and by wind 
turbine manufacturers, the following variables 
has been combined: rotor diameter; wind speed; 
and turbine pitch angle to verify its relationship 
with the rotation and power available in the hub of 
wind rotors. Such combination was built from the 
modeling of the equations that represent the wind 
generation with the scenarios being built under 
central composite rotational design (CCRD). The 
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results show a combined dependence of the constructive and environmental quantities with 
the degree of dependence presented by two theoretical equations. The conclusion of this 
research shows the need to expand studies to validate the mathematical dependence between 
shaft rotation and pitch angle in order to increase the accuracy of their control systems. 
KEYWORDS: Wind speed. Wind energy. Rotation. 


11 INTRODUÇÃO 


A demanda por eletricidade aumenta cada vez mais e a produção de energia a partir 
de fontes renováveis torna-se essencial, sendo a energia eólica uma delas. É a energia 
existente na movimentação dos ventos, ou seja, é a energia cinética contida nas massas 
de ar da atmosfera. A energia contida no vento é não poluente como outros tipos de fontes, 
que são dependentes de combustíveis fósseis, tais como gás natural ou carvão. Para 
Montezano (2007) é uma alternativa energética competitiva na qual seu grande incentivo 
acontece na forma de baixo impacto ambiental. 

O problema apresentado nesta pesquisa é representar matematicamente a rotação 
e a potência do rotor eólico de turbinas de eixo horizontal para controlar as variações e os 
efeitos das variáveis ambientais e das variáveis construtivas de um gerador sobre seu eixo, 


no ponto de acoplamento à máquina geradora. 


21 DESENVOLVIMENTO 


A pesquisa teve abordagem quantitativa descritiva apresentando as características 
da representação matemática de um gerador eólico na rotação do eixo acoplado ao cubo 
rotor eólico, estabelecendo correlações estatísticas de aderência entre variáveis analisadas, 
seguindo o propósito de uma análise descritiva (VERGARA, 1998). 

A técnica utilizada para definição da equação de regressão foi a otimização de 
processos apresentada por Rodrigues e lemma (2009). O delineamento das variáveis foi 
fatorial do tipo composto central rotacional caracterizado por quinze ensaios dos quais oito 
compuseram o modelo fatorial, seis representaram os pontos axiais e um retratou o ponto 
central. 

As variáveis respostas buscadas nessa pesquisa foram a rotação do eixo do cubo do 
rotor (R), cuja unidade é radianos por segundo (rad/s) e a potência mecânica no eixo do rotor, 
unidade watt (W). 

A construção das variáveis é dada por 2" cenários fatoriais, acrescida de uma análise 
de ponto central e oito cenários com pontos axiais. A forma de definir a expansão dos pontos 
axiais é relacioná-los ao valor fatorial nominal pelo valor de transformação da variável axial 


(a) para os pontos extremos, definida pela Eq. (1). 
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a = + (1) 


Os níveis de variação foram escolhidos com base nas informações disponíveis pelos 


fabricantes de turbinas eólicas, conforme é apresentado na Tabela 1. 


Variáveis -O -1 0 +1 +a 
D: Diâmetro do rotor (m) 70,00 100,40 145,00 189,60 220,00 
B: Ângulo de passo (º) 0,00 4,05 10,00 15,95 20,00 
V.: Velocidade vento (m/s) 3,00 7,46 14,00 20,54 25,00 


Tabela 1 — Níveis para o DCCR 
Fonte: Autoria própria (2020) 


A potência transferida do vento para o eixo da turbina eólica pode ser calculada pela 
Eq. (2). 


AxC,x Vê 
Re dis Topê [w) (2) 


Onde p representa a densidade do ar, A representa a área do rotor, C, o coeficiente de 
potência e V, a velocidade do vento. 
O coeficiente C, apresentado na Eq. (3) é adimensional. 


=C 


Cp = (5 — Cab — c4) ed (3) 


Onde À, é a relação entre a velocidade tangencial de ponta da pá com a velocidade do vento 
no ponto “” (adimensional); B é o ângulo de passo da turbina (*) e C, C,, C,C e C são 
parâmetros de ajuste de ganho. 

Neste estudo foram adotados os parâmetros de ajuste de ganho indicados por 
Montezano (2007), que apresentou os valores de ajuste Ótimos como sendo: C = 0,5176; 
C=98; C, =0,4; C= be C, = 16,5. 


31 RESULTADOS 


Os cenários mostrados na Tabela 1 são uma relação dos valores apresentados na 


Tabela 1 e ilustram a combinação para as variáveis respostas estudadas. 
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Cenário 


Oo No Aro Na 


ad do + A 
BB ONO 


15 


Tabela 1 — Cenários construídos para o DCCR para potência e rotação R. 


D(m) 
100,40 
189,60 
100,40 
189,60 
100,40 
189,60 
100,40 
189,60 
70,00 
220,00 
145,00 
145,00 
145,00 
145,00 
145,00 


B (9) 
4,05 
4,05 
15595 
15,95 
4,05 
4,05 
15,95 
15,95 
10,00 
10,00 
0,00 
20,00 
10,00 
10,00 
10,00 


V, (m/s) 
7,46 
7,46 
7,46 
7,46 
20,54 
20,54 
20,54 

20,54 
14,00 
14,00 
14,00 
14,00 
3,00 
25,00 
14,00 


Pm (MW) 
0,74 
2,63 
-0,18 
-0,63 
15,39 
54,89 
-3,68 
AB 
1,06 
10,45 
12,43 
7,92 
0,04 
25,84 
4,54 


Fonte: Autoria própria (2020). 


3.1 POTÊNCIA TRANSFERIDA PARA O ROTOR 


R (rad/s) 
119 
0,63 
1,19 
0,63 
RETA 
LRVAS: 
der 
IRTÉS: 
3,20 
1,02 
1,54 
1,54 
0,33 
2,16 
1,54 


Os resultados obtidos por meio da simulação apresentaram Pm com valores 


negativos nos cenários 3, 4, 7, 8 e 12 mostrando que altos valores aplicados ao ângulo 


de passo da turbina podem levar o rotor a ter seu giro sustentado pela energia elétrica da 


rede ao invés da energia cinética do vento, ou seja, vai consumir energia ao invés de gerar. 


Montezano (2007) explicou a relação de C, com B por meio da relação entre a 


velocidade da ponta da pá e a velocidade do vento, também definido pelo autor como 


lambda (A), com os valores apresentados na Figura 1. 


Coleção desafios das engenharias: Engenharia elétrica 2 


Capítulo 14 


0.6 


0.4 


cp 0.2 


0.0 





-0.2 
0.0 2.5 5.0 5 10.0 12.5 15.0 17.5 
lambda 


Figurai — Dependência do coeficiente de potência com o ângulo de passo da pá (B) e da velocidade 
de ponta de pá com a velocidade do vento 


Fonte: Montezano (2007) 


Para Montezano (2007), o coeficiente de potência difere de valor nulo com uma 
velocidade de ponta de pá aproximadamente 2,5 maior que a que a velocidade do vento 
para ângulos de passo testados entre 0º e 50º. A Figura 1 deixa claro que existe um valor 
ótimo de lambda para cada ângulo. Como por exemplo, para um valor de ângulo B de 0º o 
aproveitamento máximo é superior a 40% e quando esse ângulo aumenta o aproveitamento 
decresce. 

O autor também observou que À elevado pode tornar o C, negativo, fator esse que 
implica em uma potência negativa no aerogerador causando o fenômeno de consumo de 
energia elétrica para sustentação do movimento das pás. 

O cenário 6 apresentou valor de potência elevado e não constatado nos padrões 
construtivos atuais, denotado pela combinação de pequeno ângulo de passo e elevada 
velocidade do vento em uma turbina construída com pás com valor elevado de raio observa- 
se que a tecnologia atualmente existente não consegue atender tais condições construtivas 
e ambientais. 

Na simulação realizada, a potência (em MW) apresentou dependência diretamente 
proporcional à velocidade do vento V e inversamente proporcional ao ângulo de passo B, 
sendo representada pela Equação de regressão (4). 


Pm = (—9,19)B + (6,91)V, — (10,36)6 + V, + (6,83) (4) 


A Figura 2 apresentada a seguir mostra a superfície de resposta da potência máxima 
transferida para o aerogerador sob as variáveis independentes velocidade do vento e ângulo 
de passo. 
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Pm (MW) 
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Figura 2 — Superfície de resposta da potência com influências da velocidade do vento () e do ângulo 
de passo (f). 


Fonte: Autoria própria (2020). 


Esta figura apresenta a dependência diretamente proporcional da potência gerada 
pela turbina eólica com a velocidade do vento. Também caracteriza que elevado ângulo de 
passo reduz o rendimento da turbina, podendo fazer com que a usina deixe de produzir 
energia e, para manter suas hélices girando, passa a consumi-la. 

Nesse mesmo contexto Minnesota Municipal Power Agency [20--] também descreveu 
que para cada ângulo de passo diferente há um limite de relação À que torna o coeficiente 
negativo. Esse coeficiente negativo faz com que a turbina deixe de fornecer energia ativa 


para a rede. Por esse motivo os projetos de geração eólica devem prever o corte de geração 
para altas velocidades de vento. 


3.2 ROTAÇÃO NO EIXO DA TURBINA 


O diâmetro de um rotor varia de acordo com a potência desejada para o aerogerador. 


A Tabela 3 apresenta alguns valores construtivos de diâmetro para aerogeradores entre 3 e 
4,5 MW de potência de pico de geração. 
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Potência da turbina Diâmetro do Rotação do Área varrida 


Modo (MW) rotor (m) cubo (RPM) do rotor (m2) 
N155 A5 155 1 18.870 
N155 48 155 1 18.870 
N131 3,9 131 7,9 15,4 13.478 
N131 3,6 131 7,5 13,6 13.478 

AW148 323,6 148 6,2 12,4 17.203 

AW132 3,465 132 744 13.685 


Tabela 3 - Exemplificação do diâmetro de alguns modelos de turbinas eólicas. 


Fonte: Adaptado do site da fabricante Nordex (2020). 


Os dados da Tabela 3 não esgotam todas as possibilidades construtivas e servem 
apenas como exemplos trazidos de modelos utilizados por fabricantes de aerogeradores 
para facilitar as relações matemáticas aqui apresentadas. 

A literatura não apresenta uma equação para modelar a rotação do eixo, sendo 
que para esta simulação foi constatada uma dependência diretamente proporcional 
a velocidade do vento e inversamente proporcional ao diâmetro do rotor. A regressão 


estatística é apresentada na Eq. (5). 


R = (-0,58)D + (0,77)V, — (0,25)D * V, + (1,71) (5) 


A superfície de resposta da rotação do eixo do aerogerador sob as variáveis 
independentes velocidade do vento e diâmetro do rotor é apresentada na Figura 3, permitindo 
observar a dependência diretamente proporcional da rotação com a velocidade do vento e a 


dependência inversa da resposta (rotação) com o diâmetro do rotor. 
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Rotação (rad/s) 





Figura 3 - Superfície de resposta da Rotação (R) com influências da velocidade do vento (Vv) e do 
diâmetro do rotor. 


Fonte - Autoria própria (2020). 


41 CONCLUSÃO 


A pesquisa teve como propósito encontrar uma correlação das variáveis diâmetro do 
rotor, velocidade do vento e ângulo de passo, pois possuem atuação direta para a geração de 
energia eólica e assim influenciando a fabricação de modelos de turbinas eólicas. 

Foi possível verificar que as variáveis ambientais e construtivas possuem grande 
influência na rotação e potência disponível no eixo do rotor eólico pelos gráficos de superfícies, 
sendo necessária a realização de experimentos para validar as equações de regressão aqui 
apresentadas. 

Foi observado nesse estudo que o ângulo de passo é a variável imprescindível para 
o controle da potência transferida para turbina eólica quando da variação na velocidade 
do vento, sendo que essa hipótese constitui a recomendação para trabalhos futuros com 
objetivo de verificar existência de dependência matemática entre a rotação do eixo e o ângulo 


de passo da pá do aerogerador. 
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RESUMO: Uma das mais utilizadas energias 
renováveis no mundo é a do sistema fotovoltaico, 
por conta do seu desenvolvimento energético e 
sua economia financeira. Por suas vantagens, 
foi realizada implantação desse sistema em 
uma instituição de ensino superior situada no 
recôncavo baiano. Nesse contexto, o objetivo 
desse estudo consiste em avaliar a viabilidade 
do sistema fotovoltaico implantado na Faculdade 
Maria Milza localizado no centro do recôncavo 
baiano. Foram coletados dados, no qual foi 
realizado um levantamento da energia fornecida 
pela concessionária e o emprego de dados do 
sistema de energia fotovoltaica. O método de 
procedimento foi o estudo de caso. Foi utilizado 
o software R, um programa de estatística, 
para análise dos dados das contas de energia 
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recolhidos durante o período de Janeiro à 
Dezembro do ano de 2019, ainda na utilização da 
energia elétrica proveniente da concessionária 
baiana, COELBA, para se chegar a uma média 
de consumo e de custo de quilowatit-hora. Em 
seguida foram recolhidos dados de potencial 
solar do recôncavo com a utilização do programa 
Meteonorm 7.1 para o sistema fotovoltaico, 
obtendo assim uma comparação econômica 
e de eficiência entre ambas. Os resultados 
abordados apontam que com o investimento no 
sistema solar fotovoltaico, com vida útil de 10 a 
20 anos na instituição de ensino superior, teria- 
se um retorno financeiro do investimento em 
dois anos e um mês, uma produção superior a 
45.000 kWh, suficiente para toda a demanda e 
assim tendo uma grande economia significativa 
em comparação com o abastecimento da 
concessionária. Com isso, podemos concluir 
que está pesquisa, demonstrou que com um 
investimento no sistema solar fotovoltaico na 
Faculdade Maria Milza localizada no recôncavo 
da Bahia seria de máxima economia para o 
sistema de abastecimento elétrico da empresa, 
tanto quanto financeira, tendo em vista que a 
região é propicia para o desenvolvimento de 
energia solar, assim também contribuindo com a 
sustentabilidade ambiental. 
PALAVRAS-CHAVE: Energia solar. Eletricidade. 
Engenharia. Meio Ambiente. Sustentabilidade. 


FEASIBILITY OF THE PHOTOVOLTAIC 
SYSTEM IN THE RECONCAVO DA BAHIA 
REGION 


ABSTRACT: One of the most used renewable 
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energies in the world is the photovoltaic system, due to its energy development and financial 
savings. Due to its advantages, this system was implemented in a higher education institution 
located in the Bahian Recôncavo. In this context, the objective of this study is to evaluate the 
feasibility of the photovoltaic system implemented at Faculdade Maria Milza located in the 
center of Bahia's Recôncavo. Data were collected, in which a survey of the energy supplied by 
the concessionaire and the use of data from the photovoltaic energy system was carried out. 
The procedure method was the case study. The R software, a statistical program, was used 
to analyze the data from energy bills collected during the period from January to December 
of 2019, still in the use of electricity ffom the Bahia concessionaire, COELBA, to arrive at a 
average consumption and cost of kilowatt-hours. Then, solar potential data were collected 
from the Recôncavo using the Meteonorm 7.1 program for the photovoltaic system, thus 
obtaining an economic and efficiency comparison between them. The results discussed show 
that with the investment in the photovoltaic solar system, with a useful life of 10 to 20 years 
in the higher education institution, there would be a financial return on the investment in two 
years and one month, a production of more than 45,000 kWh, sufficient for all the demand and 
thus having a great significant savings compared to the supply of the concessionaire. Thus, 
we can conclude that this research has shown that with an investment in the photovoltaic 
solar system at Faculdade Maria Milza located in Bahia's Recôncavo, it would be of maximum 
savings for the company's electrical supply system, as well as financially, given that the region 
it is conducive to the development of solar energy, thus also contributing to environmental 
sustainability. 

KEYWORDS: Solar energy. Electricity. Engineering. Environment. Sustainability. 


11 INTRODUÇÃO 


No atual cenário mundial, com as diversas fontes de energia não renováveis, fez- 
se necessário que a engenharia buscasse propror novas fontes de energias mais limpas 
e renováveis para o consumo humano, desta forma foi-se abrangendo as diversas fontes 
providas da natureza para o abastecimento, uma delas é a energia fotovoltaica, oriunda de 
raios solares. 

O termo energia sustentável pode ser usado para denominar tecnologias e recursos 
energéticos, que além de serem capazes de suprir as necessidades da humanidade, mais 
que também o façam de forma adequada com a preservação e a integridade dos sistemas 
naturais essenciais, evitando as drásticas mudanças climáticas na terra (BORBA; GASPAR, 
2007). 

A energia solar é inesgotável, a escala de tempo terrestre para o ser humano, 
e deve ser aproveitado ao máximo como fonte de calor e de luz, sendo que hoje ela é 
uma das alternativas energéticas mais promissoras para fornecer energia necessária ao 
engrandecimento da humanidade (PINHO; GALDINO, 2014). 

Dentre todos os projetos, a visão central da energia elétrica e seu abastecimento 
caracterizam-se pela eficiência e o desenvolvimento do Brasil, que tem como objetivo 


principal a eficácia no processo produtivo, reduzindo o consumo de energia fornecida pelas 
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concessionárias em 40%, tudo isso graças ao setor industrial que trabalha para que haja 
maior segurança no atendimento no quesito da demanda de energia com destaque nas 
fontes renováveis, pois elas são de extrema importância para fortalecimento da indústria 
brasileira no mercado global (FGV, 2016). 

O sistema fotovoltaico, na qual se conhece atualmente, teve seu primeiro indício 
descoberto por Edmond Becquerel em 1839 através do aparecimento de uma diferença 
de potencial nos terminais de uma célula eletroquímica, acontecimento proveniente de 
absorção de luz. Posteriormente, em 1876 surgiu o primeiro equipamento fotovoltaico após 
estudos da física do estado sólido (SILVA; SOUZA, 2017). 

No ano 1956 houve um grande impulso para o crescimento desse sistema, pois foi 
quando ocorreu a evolução da área eletrônica e tecnologia, onde houve um aumento no 
setor de telecomunicações, de fontes de energia para sistemas instalados em localidades 
distantes e a corrida espacial que também contribuiu para esse feito, pois fornecia menor 
custo, peso e segurança, além de fornecer a energia necessária para as viagens espaciais, 
o que continua sendo até os dias atuais (PINHO; GALDINO, 2014). 

A energia solar fotovoltaica é obtida através dos materiais semicondutores que captam 
e convertem a incidência da luz solar transformando em eletricidade, conectando a rede de 
distribuição da edificação residencial, comercial ou industrial, fazendo com que a própria 
edificação produza sua energia e seja auto-sustentável (CARVALHO JÚNIOR, 2019). 

O funcionamento do sistema é bastante simples e explica que durante o dia os painéis 
absorvem a energia da luz solar e a mesma é convertida em energia elétrica, corrente 
contínua (CC) e o inversor a transforma em corrente alternada (CA), que pode ser consumida 
nas edificações, gerando energia para lâmpadas, eletrodomésticos e qualquer aparelho 
conectado à rede (CARVALHO JÚNIOR, 2019). 

Existem dois tipos de sistemas fotovoltaicos, de acordo com a Bluesol (2017, p.37) “a 
resposta simples e direta é: o gerador solar fotovoltaico on-grid é ligado à rede elétrica, e o 
gerador fotovoltaico off-grid não é ligado à rede elétrica.” 

A Instalação do sistema solar fotovoltaico on-grid interligado à rede elétrica pode se 
apresentar de duas formas diferentes, sendo que uma pode ser instalada de forma integrada a 
uma edificação ou de forma centralizada como em uma usina central geradora convencional, 
onde neste caso geralmente está um pouco distante do ponto de consumo (RUTHER, 2004). 

Conforme a Câmara Brasileira da Indústria da Construção (CBIC) de todas as fontes 
de energias renováveis, a energia solar fotovoltaica contém uma gama de vantagens, como 
por exemplo, exercer um papel fundamental em ajudar as hidrelétricas, que são as principais 
fontes de geração de energia do Brasil, diminuindo significativamente os aumentos de pico de 
demanda de energia durante o dia, não gera nenhuma emissão poluente ao produzir energia 
elétrica, sem a necessidade de uso de combustível o que minimiza drasticamente o custo. 
Contudo, a geração pelo sistema fotovoltaico pode ser gerada nos locais onde se deseja 


consumir a carga, reduzindo as linhas de transmissão e aumento da segurança energética 
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(BRASIL, 2016). 

Quando se produz mais energia pelo sistema fotovoltaico do que se consome essa 
energia, ela pode ser injetada na rede de distribuição da concessionária por meio de um 
relógio de luz, que tem por objetivo medir a energia da rua quando não há sol e também a 
energia solar gerada em excesso, a edificação assim ganha um crédito de energia (TRAPANI, 
et. al, 2013). 

Esse crédito é regulamentado pela Agência Nacional de Energia Elétrica (ANEEL), 
que enfatiza o fato de haver regras especificas que variam de acordo com a localização e a 
classe de consumo. A agência mede por cada KWh gerado em excesso pelo sistema solar 
fotovoltaico da edificação e o consumidor recebe então o seu crédito em KWh que poderá ser 
utilizado de noite ou nos meses seguintes (CARVALHO JÚNIOR, 2008). 

Há muitas vantagens para as edificações urbanas e o sistema de distribuição ser 
integradas ao sistema fotovoltaico, incluindo uma série de vantagens para o sistema elétrico, 
no qual se relaciona com a economia, que não vem sendo considerado e quantificado. 
(RUTHER, 2014). 

Entre as vantagens do sistema solar fotovoltaico podem ser citadas que as perdas por 
transmissão e distribuição de energia são minimizadas, redução de investimentos em linhas 
de transmissão e distribuição, edifícios solares fotovoltaicos não apresentam necessidade de 
àrea física dedicada, já que a mesma é ocupada pela edificação, os edifícios têm capacidade 
de oferecer suporte KVAR a pontos críticos da rede de distribuição melhorando a qualidade 
da energia, capacidade de oferecer um grande fator de capacidade e alimentadores da rede 
com picos diurnos, entre outros (KALOGIROU, 2016). 

Nesse contexto, a Faculdade Maria Milza (FAMAM), situada no reconcavo baiano 
no território de Governador Mangabeira, localizada na rodovia BR 101, procurou ter um 
papel de suma importância no requisito de energias sustentáveis, para suprir seu proprio 
guarnecimento e promover a sustentabilidade ambiental na região. 

A FAMAM sempre foi abastecida de energia elétrica pela concecionária do estado da 
Bahia, a Companhia de Eletricidade do Estado da Bahia (COELBA), porém no início do ano 
de 2020, a faculdade tomou a iniciativa de implantar um projeto de viabilidade energética 
através de energias renováveis, pretendendo instalar painéis fotovoltaicos do sistema on-grid 
para seu próprio abstecimento em um período de 20 anos. Nessa perspectiva, o objetivo 
desta pesquisa foi analisar a viabilidade economia e a eficiencia do sisitema de energia solar 
na região da Faculdade Maria Milza pela produção e o próprio abastecimento de energia 


elétrica provida do sistema fotovoltaico on-grid a ser instalado. 


2 | MATERIAL E MÉTODOS 


A metodologia que desenvolveu esta pesquisa teve como modalidade o estudo de 
caso de abordagem quantitativa. A localização onde a pesquisa ocorreu foi no Recôncavo 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 15 





Baiano, na cidade de Governador Mangabeira na BR 101 onde se encontra a instituição de 
ensino superior referida, a FAMAM. 

Toda instituição ou empresa tem um contrato de fornecimento de demanda de potência 
contratada com a concercionária na qual não se pode abater o fornecimento produzido no 
sistema solar, esse contrato de potência a instituição ou empresa ainda irá pagar normalmente, 
o sistema fotovoltaico apenas abaterá a utilização de energia de consumo ativo na ponta e 
consumo ativo fora da ponta da faculdade. 

A aplicação da viabilidade com a instalação do sistema fotovoltaico on-grid foi 
devidamente estudada no consumo ativo na ponta, que é o periodo onde o horário da 
procura de energia da concerssionaria é mais requisitado, tendo uma diaria consecutiva 
de 3 horas e consumo ativo fora da ponta, que é o período diário composto pelas horas 
consecutivas e complementares ao horário de ponta, abatendo assim todos os valores 
financeiramente ao da demanda que era fornecida pela concercionária e os valores de 
consumo em kWh. 

Foram recolhidos os dados das contas de energia durante o período de Janeiro 
a Dezembro do ano de 2019, ainda na utilzação da energia elétrica proveniente da 
concessionária COELBA, para se chegar a uma média de consumo e de custo de quilowatt- 
hora, foi utilizado o software R, um programa de estatística. 

Foram recolhidos dados de potencial solar da região do recôncavo com a utilização 
do programa Meteonorm 7.1 para o sistema fotovoltaico, obtendo assim uma comparação 
econômica e de eficiência entre ambas. 

Para a avaliação do potencial solar, foram escolhidos os dados climáticos e 
solarimétricos da base de dados Meteonorm 7.1, disponíveis no software de modelagem 
solar PVsyst. A base de dados meteorológicos do Meteonorm é privada e baseada em dados 
de estações em solo e dados de satélite. 

Para o caso específico deste estudo, os dados disponibilizados foram obtidos no ano 
de 2019 do Meteonorm, versão 7.1, e são constituídos de médias mensais de longo prazo. O 
software faz a interpolação de dados de satélite com dados de medições das três estações 
meteorológicas mais próximas ao local da instalação. 

Silva (2013) afirma que para estimar a produção média mensal de energia de um 


sistema fotovoltaico, pode ser encontrada utilizando a expressão da equação 1 a seguir: 


GMM (kWh) = Ptot (kWp) x RSM (h/dia) x PR x nº dias (1) 

Onde: 

GMM (kWh) é o valor médio de geração mensal. Diferentes épocas do ano geram 
diferentes valores mensais, uma vez que o nível de radiação solar varia bastante durante 
o ano, devido à mudança de posicionamento da terra em relação ao sol (estações do ano). 

Ptot (KkWp) é a potência total instalada em placas fotovoltaicas, dada pela multiplicação 


da potência de cada placa pelo número de placas. 
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RSM (h/dia) é a radiação solar média no período de um ano para o local onde o 
gerador será implantando. 

Neste caso, NMRS = É a região do recôncavo baiano, considerando o número de 
horas de sol pleno por dia. 

PR é a performance ratio ou eficiência global do sistema. Esse valor representa 
o percentual de eficiência do gerador como um todo, levando em conta todas as perdas 
inerentes aos equipamentos e a operação do sistema, como perda no inversor, perda nos 
condutores, queda de tensão, indisponibilidade do sistema etc., e às condições do local da 
instalação, como perdas por sombreamento, perdas por sujeiras nos módulos, variações de 


temperatura ambiente etc. 


31 RESULTADOS E DISCUSSÃO 


Diferentes épocas do ano geram diferentes valores mensais, uma vez que o nível 
de radiação solar varia bastante durante o ano, devido à mudança de posicionamento da 
terra em relação ao sol (estações do ano). No Recôncavo, devido a sua localização na 
zona tropical, com duas estações bem definidas (verão e inverno) a incidência solar direta 
ocorre durante seis meses (verão) e que no inverno não há uma diminuição significativa da 
incidência dos raios solares nessa área, o que não compromete significativamente o sistema 
fotovoltaico, onde se encontra localizada a Faculdade Maria Milza. 

O nordeste tem um grande potencial a respeito no desenvolvimento de produção 
de energia solar no Brasil, por exibir níveis elevados de insolação com baixa variabilidade 
durante o ano. (PEREIRA et al., 2006). 

Foi verificada uma radiação média ou NMRS de 5,04 h/dia ou 5,04 kWh/m?/dia, valor 
retirado do programa Meteonorm, referente ao ano de 2019, conforme a Figura 1. A partir do 
levantamento dos dados meteorológicos referentes à região do Recôncavo da Bahia, obtidos 
por meio do programa Meteonorm, referente ao número de irradiação solar, desta forma foi 
possível estimar os valores da fórmula de produção média mensal de energia do sistema 
fotovoltaico. 


Coleção desafios das engenharias: Engenharia elétrica 2 Capítulo 15 





Site Governador Mangabeira [Braz] 


[ha eteororm 7.2. Sat=B55 


Data SOUNce 


Horizontal Horizontal Temperature a tumocd Linke Relatire 
global diff se vw ebocity Turbidity Humidity 

radiation radiation 

EsithAmê. dam E innÊ. cas “4E ns: [-] pa 
amitary lsz0 [263 EFE ja30 [2086 [810 
February jsos [2610 266; pes [2876 [soo 
“Prer 52º 3 E 7 [520º | EC 
April [a PET ERR [240 ETR ES 
fito 555º — o Foco Er [725 Ecs 
ma EE 5 = EXTOR PEA ESC 
Rr FT EXE =o E5 Des ETR 
September [74 | P5 Bão Bam EXE EDER 
October EEE 7 EST [570º ETTA ETR 
November [E — EE BSS [Po [5220 57 
December [642 [240 EEE [34 ja3470 FE 
“Tear EA 2 34 25.0 1 2 986 g1.9 


Figura 1: Dados meteorológicos e irradiação solar diária média da região do Recôncavo da Bahia, 
obtidos por meio do programa Meteonorm, referente ao ano de 2019. 


Fonte: Dados da pesquisa, 2020. 


Segundo a Solenerg (2012), com os dados meteorológicos podemos assim encontrar 
os valores da equação e estimar a produção do sistema fotovoltaico de acordo com o autor 
supracitado Silva (2013). 

O Ptot (kWp) nesta instalação da instituição, o gerador fotovoltaico será produzido por 
1.200 módulos de silício poli cristalino de 335 Wp, 144 células, resultando em uma potência 
instalada de 402 kWp. Neste caso, Ptotal = 402 kWp. Também foi estimado que com os 1.200 
módulos, de medidas (1.992 mm x 992 mm / cada) ocuparam uma área total de 2.405 m?. 

O RSM (h/dia) como já foi citado, refere-se a um valor que é retirado de algum banco 
de dados solarimétricos e, neste projeto, foi considerado o banco do CRESESB. Já PR 
considerado foi de 75,8%, que é um valor razoável para considerar todas as perdas, já que 
não é possível mensurar todas elas fisicamente. 

Assim, tomando como base a fórmula e os dados anteriormente citados, a produção 
média mensal estimada para o gerador do projeto foi de GMM = 402 x 5,04 x 0,758 x 30 = 
46.073 kWh. 

No programa de estatística R, com a introdução dos dados de consumo mensal durante 
Janeiro à Dezembro do ano de 2019 da Instituição FAMAM ainda com o fornecimento da 
concessionária de energia COELBA, obteve-se uma média de consumo mensal de 41.010,35 
kWh. 

Comparando então o consumo médio mensal do período de Janeiro a Dezembro do 
ano de 2019 de consumo que era fornecido pela concessionária COELBA em que se verificou 
um resultado de 41.010,35 kWh, isso somando o consumo ativo na ponta e fora da ponta, à 
produção de energia fotovoltaica é de 46.073 kWh, sendo abatido no consumo ativo fora da 
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ponta, demonstrando que os 1.200 módulos foram suficientes para suprir todo o consumo 
energético da Faculdade Maria Milza, resultando ainda em um sobra de 5.062,65 kWh, na 


qual abateu no valor de consumo ativo na ponta, como mostra o Quadro 1. 


ESTIMATIVA DE PRODUÇÃO PELO SISTEMA FOTOVOLTAICO 
si E E 


[| Jamero | zaasão | ao | 622 [og | 3 | sara] 
Dam | anseass | ao | a33 Joss | 30 |  aossz6r | 
| mão | azssaas | ao | am | os | a | 3500539] 
Junho | aos | ao | 336 | 0758 | 30 |  aom537] 


| Juno | asas | ao | 365 [os | 3 |  aaazaeo 
| Setembro | azemso | ao | am | os | 30 |  azaaras 
| Novembro | anasas | ao | 625 | 0758 | 30 |  s7a3425] 
TOTAL (kWh) 

Meses | Do | Meses | 
Médiagewnjjmês| mooss 


Quadro 1: Diferença entre o consumo energético da concessionária e a estimativa da produção de 
energia pelo sistema fotovoltaico. 





Fonte: Dados da pesquisa, 2020. 


O investimento em um sistema fotovoltaico de 1.200 módulos, com todos os 
componentes como inversores, transformadores, dispositivos de proteção CC/CA, as 
estruturas metálicas, cabos de isolamento, conectores e sistema de aterramento para a 
instituição FAMAM, verificou-se um valor médio de R$ 1.330.000,00. Conforme o Quadro 2 


pode-se estimar o retorno financeiro investido, levando-se em consideração o aumento da 
tarifa de energia de 11% ao ano. 


Energia Gerada Valor Geração acumulada 
E ES 


538.577,14 3.802.913,52 2472 013,52 
535. 997,65 4.613.975,62 3.283.975,63 


E 


533.430,52 SA 4347,75 4.144 347,76 


518.283,85 11.540 508,10 0.310.508,10 
515.801 12.832.502,47 11.502.502 47 
1 


E 5 
528.333,07 7.341 061,90 011 061,90 
525.802,64 s.346. 432,72 016. 432,72 
523.284,34 9394 852,14 8. 064. 832,14 
513.331,14 14. 069 978,06 12. 759978,06 


.58 538,9 

520.778,09 1.330.000,00 10.494 461,37 9.164.461,37 
E 508. 1 
.55 S024 


| 
7 
| 
| 

| asióssio 


530.875 6.383.538,97 053.538,97 


510.872,57 15.352. 471,88 14.022 471,88 


EEE 
[6 | 
E 
DE: 
9 | 





Quadro 2: Estimativa de retorno financeiro do sistema fotovoltaico na Instituição FAMAM localizada no 
recôncavo baiano no período de 2020 com o aumento anual da tarifa de energia de 11%. 


Fonte: Dados da pesquisa, 2020. 
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A partir deste Quadro 2 pode-se ver que com um investimento de R$ 1.330.000,00, 
em aproximadamente dois anos e um mês a economia gerada é de R$1.364.288,02 onde se 
tem um retorno financeiro do investimento inicial por completo. 

Segundo Freitas e Miranda (2016) a tecnologia do sistema fotovoltaico vem ganhando 
espaço no mercado do Brasil, embora possua custos elevados e o retorno seja em longo 
prazo, ainda assim compensa o valor investido para obtenção do sistema. 

A utilização desse tipo de sistema traz grandes benefícios aos usuários, onde através 
do cálculo de valor unitário mostra que a energia produzida pelo sistema fotovoltaico tem 
menor custo do que a fornecida por concessionária, assim o usuário economizará no consumo 
da energia que o seu sistema produz (DANTAS; POMPERMAYER, 2018). 


41 CONCLUSÃO 


O sistema fotovoltaico é capaz de gerar energia elétrica através de seus módulos, que 
captam a radiação proveniente do sol. Esse é um sistema sustentável de grande importância, 
no qual, produz energia elétrica sem poluir ou prejudicar o meio ambiente, tem sido um dos 
grandes geradores de energia ao redor do planeta. 

Para uma boa produção de energia elétrica através do sistema fotovoltaico, 
obviamente é preciso que o sol incida nas placas ou módulos para sua captação durante o 
dia, nesse sentido existem regiões do Brasil no qual a incidência solar é mais propicia para 
esse sistema. 

O recôncavo baiano demonstrou ter um grande fator climático para o desenvolvimento, 
abastecimento e consumo de energias geradas pelo sistema fotovoltaico na região. Com isso 
verificou-se que na Faculdade Maria Milza, ao possui uma grande área para essa produção, 
levando em conta que a instituição possui uma estrutura que possibilitou todas as instalações 
dos 1.200 módulos, totalizando uma área de 2.405 m? no qual é necessária para suprir toda 
sua demanda energética. 

A instalação do sistema de energia fotovoltaico on-grid suprirã toda a demanda da 
Faculdade Maria Milza de consumo ativo fora de ponta e consumo ativo na ponta, uma vez 
que, deixa de pagar o consumo de energia para a concessionária, pagando apenas o valor 
de contratado de potência. 

Então com o investimento para a instalação do sistema fotovoltaico on-grid na 
instituição superior de ensino FAMAM, será produzido em média 46.073 kWh de demanda de 
consumo elétrico, que supriu toda a necessidade de consumo da instituição que em média 
era de 41.010,35 kWh, assim aliviando financeiramente em questão energética de consumo 
na ponta e fora da ponta na qual era contabilizada para a concessionária do estado da Bahia 
e que agora é produzida pelo sistema solar. 

O investimento financeiro de R$ 1.330.000,00 com a instalação do sistema fotovoltaico, 


obteve-se um retorno em dois anos e um mês, no qual os sistemas têm uma garantia de 
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fabricação de mais de dez anos. 

Concluindo assim que com a instalação do sistema fotovoltaico on-grid na instituição 
de ensino superior da FAMAM localizada no recôncavo baiano recolherá grandes benefícios 
financeiros e energéticos com a produção de energia solar, uma vez que possui grande 
espaço e estrutura para as placas geradoras, reduzindo o pagamento da conta de demanda 
energética da concessionária COELBA e produzirá sua própria fonte de consumo energética 
limpa e sustentável durante anos., por estar localizada em uma região onde a geração de 
energia por meio do sistema solara é bastante vantajoso. 
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